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SECTION  I 
INTRODUCTION 

Detailed  theoretical  investigation  of  the  principle 
of  Auto-Resonant  acceleration  has  proceeded  to  the  point 
where  a detailed  experimental  test  of  this  accelerator 
concept  can  be  formulated.  Such  an  experiment  would  have 
three  fundamental  goals: 

1.  Examine  and  verify  the  basic  concept  of  Auto- 
Resonant  acceleration. 

2.  Explore  the  nonlinear  limits  of  operation  of 
such  an  accelerator. 

3.  Establish  the  proper  scaling  laws  of  such  a 
device  to  higher  output  values  of  ion  current 
and  energy. 

A conceptual  experiment  to  undertake  these  tasks 
has  been  devised  under  the  following  basic  guidelines: 
First,  in  order  to  minimize  technological  problem  areas 
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not  directly  associated  with  the  physics  of  the  experiment, 
we  require  that  the  electron  beam  power  source  be  of  a 
standard,  state-of-the-art  type  not  requiring  new 
technological  innovations,  that  magnetic  guide  fields  be 
of  reasonable  values  easily  achievable  in  the  laboratory, 
that  total  accelerator  length  be  held  within  reasonable 
bounds,  and  that  frequencies  and  wavelengths  associated 
with  phenomena  of  interest  lie  in  a regime  of  easy  experi- 
mental access.  At  the  same  time,  the  ion  output  energy 
level  should  be  sufficiently  high  to  discriminate  against 
other  types  of  collective  acceleration  processes.  Finally, 
ion  output  current  should  be  sufficiently  flexible  to  allow 
examination  of  nonlinear  saturation  phenomena. 

Using  these  basic  guidelines  and  the  parameter 
survey  code  developed  earlier  [Appendix  Z,  Reference  1 
(Reference  written  as  Ref.  herafter)],  the  following 
accelerator  parameters  were  selected: 


1.  W.  E.  Drummond,  et  al . , "A  Theoretical  Investi- 
gation of  the  Principle  of  Auto-Resonant  Acceleration," 
AFWL-TR-75-296 , Air  Force  Weapons  Laboratory,  Air  Force 
Systems  Command,  Kirtland  Air  Force  Base,  N.M.,  Nov.  1975. 
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Electron  beam  energy 

3 MeV 

Electron  beam  current 

30  kA 

Minimum  e-beam  pulse  time 

100 

ns 

Nominal  ion  output  energy 

30 

MeV 

Nominal  ion  output  current 

30 

A 

Maximum  magnetic  field 

24 

kG 

Minimum  magnetic  field 

2 

kG 

Wave  Growth  Section  Length 

5 

m 

Accelerator  Length 

4 

m 

RF  liner  frequency 

240 

MHz 

In  this  accelerator,  the  relativistic  e-beam  is 
injected  into  the  vacuum  waveguide  from  a 3 cm  radius, 
conventional  foil  diode  configuration  immersed  in  a 2.4  kG 
magnetic  field.  Following  a roughly  30  cm  transition 
region,  where  beam  equilibrium  is  established,  the  beam 
is  compressed  in  a 1 m long  compression  section  to  a 1 cm 
radius  with  a 24  kG  to  2 kG  magnetic  field  at  the  output. 

In  this  technical  report,  the  details  of  the 
conceptual  design  criteria  and  specifications  for  such  an 
accelerator  are  presented.  The  requirements  for  the 
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electron  beam  power  source  and  diode  are  discussed  in 
Section  II,  with  specifications  and  diagnostics  for  the 
propagating  electron  beam  given  in  Section  III.  Section 
IV  contains  details  of  the  magnetic  field  system  and 
configuration.  In  Section  V,  various  wave  growth  schemes 
are  examined,  with  ion  loading  and  acceleration  discussed 
in  Section  VI.  Some  overall  system  requirements  are  given 


in  Section  VII. 


SECTION  II 


POWER  SOURCE  AND  DIODE 

The  nominal  electron  beam  parameters  required  for 
a prototype  Auto-Resonant  Accelerator  are  described  in  the 
preceding  section.  The  pulsed  power  system  needed  to 
produce  this  beam,  as  well  as  the  diode  itself  and  its 
associated  diagnostics,  is  described  below. 

A.  Pulsed  Power  System 

Specifications  for  the  pulsed  power  system  are 
detailed  in  Appendix  A and  listed  below  for  completeness. 
Electron  Injector  Specification 

1.  Load:  100  ft  (constant  impedance  dummy  load 

for  tests) 

2.  Peak  voltage  (V  , ) into  100  ft:  3 MeV 

pk 

3.  Duration:  200  ns  minimum  useful  pulse  between 

0.96  V points  (±  2%) 

4.  Reproducibility  of  V ^ shot  to  shot:  ± 2% 

5.  Rise  time:  s 100  ns  (5%  to  95%  of  V . ) 

pk 
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6.  Prepulse  suppression:  Maximum  prepulse  less 

than  3 % of  V . (measured  with  no  dummy  load) 
px 

7.  Ripple:  Less  than  ± 1/2%  high  frequency 

( 2 10  MHz)  ripple  during  useful  pulse 

(V  2 0.96  V ) 

pK 

8.  Triggering:  Command  triggered  with  less  than 

± 200  ns  jitter  between  tfiggering  and  firing 

9.  Pulse  clipping  crowbar  switch  to  reduce  energy 
into  100  (2  load  to  less  than  10%  of  initial 
stored  energy  after  200  ns  useful  pulse  (0  to 
+ 50  ns  jitter  allowed  in  crowbar  switch) 

10.  The  electrical  circuit  driving  the  load  should 
have  a low  internal  impedance. 

The  voltage  and  current  requirements  (3  MeV,  30  kA)  needed 
in  the  electron  beam  are  well  within  the  range  of  conven- 
tional technology,  although  the  duration  of  the  flat-top 
portion  of  the  voltage  pulse  (200  ns)  is  longer  than 
normally  used  in  flash  X-ray  machines.  However,  it  is 
certainly  attainable  with  conventional  techniques. 

These  specifications  were  submitted  to  two  major 
manufacturers  of  pulsed  power  systems.  Both  groups 
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submitted  machine  configurations  which  could  meet  the  pulsed 
power  requirements.  To  confirm  the  suitability  of  each 
design,  in  particular  the  ripple  specifications,  computer 
generated  digital  voltage  waveforms  were  submitted  to  Austin 
Research  Associates  by  both  manufacturers.  These  waveforms 
were  modelled  with  a numerical  code  to  study  the  phase 
stability  of  the  Auto-Resonant  Accelerator.  The  details  of 
this  study  can  be  found  in  Appendix  B.  The  results  of  this 
study  show  that  both  the  droop  and  ripple  from  either 
machine  were  "safe"  by  at  least  a factor  of  2.  Thus, 
it  appears  possible  to  obtain  a pulsed  power  system  suitable 
for  studying  a prototype  accelerator. 

B.  Diode  Considerations 

Realization  of  the  uniform  beam  required  for  the 
accelerator  is  most  easily  attained  using  a conventional 
cold  cathode/thin  foil  anode  diode.  Foil  survival  then 
becomes  an  important  consideration.  Detailed  calculations 
in  Appendix  C show  that  an  A-K  gap  spacing  of  about  3 cm, 
with  a corresponding  current  density  of  1 kA/cm2,  will 
ensure  long  anode  lifetime.  The  large  gap  spacing  also 
minimizes  the  diode  impedance  collapse  due  to  plasma  closure. 
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More  detailed  work  in  the  area  is  indicated  to  ensure 
transition  to  the  appropriate  equilibrium. 

After  an  equilibrium  is  established,  the  beam  must 
be  compressed  to  a cross-sectional  size  compatible  with 
keeping  the  cyclotron  wave  growth  section  length  to  a 
tractable  size.  The  compression  should  be  carried  out  over 
a long  enough  distance  to  minimize  beam  heating  (intro- 
ducing transverse  temperature) . Appendix  F develops  the 
theory  of  beam  compression,  and  gives  a formula  for 
numerically  calculating  the  amount  of  heating  introduced 
by  a given  compression  section.  Initial  estimates  assum- 
ing an  adiabatic  compression  indicate  the  amount  of 
transverse  energy  to  be  roughly  1%  of  the  total  energy. 

This  is  sufficiently  small  to  preclude  Landau  damping  of 
the  desired  cyclotron  wave. 

B . Electron  Beam  Characterization 

The  diagnostics  which  would  be  used  to  characterize 
the  electron  beam  at  the  anode  and  in  the  beam  transition 
region  are  described  in  Appendix  G.  The  types  of 
measurements — calorimetry,  Faraday  cups,  magnetic  loops, 
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etc.— are  standard.  Also,  X-ray  techniques  can  provide  a 
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great  deal  of  information  on  gross  angular  spreading  of 
the  electron  beam. 
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I,  the  results  of  a computer  code  which  integrates  the 
relativistic,  one  dimensional  ion  equation  of  motion  are 
given  for  acceleration  of  ions  to  the  25  to  30  MeV  energy 
range.  There  it  is  shown  that  stable  acceleration  results 
if  the  spatial  ripple  is  kept  below  5%. 

In  Appendix  J an  analytical  estimate  of  the  ripple 
to  be  expected  from  physically  realizable  coils  is  given. 
These  calculations  show  that  the  ripple  may  be  kept  below 
1%  so  that,  with  careful  field  coil  design,  ripple  should 
present  no  problem. 

An  unfortunate  property  of  the  guide  field  which 


produces  the  shortest  acceleration  length  for  a given 
final  energy  is  that  it  is  nonadiabatic  ( jpT  ~ 


for  the  first  few  centimeters.  The  effects  of  replacing 
this  nonadiabatic  section  with  an  adiabatic  section  have 
been  investigated.  It  is  shown  in  Appendix  I that  this 
patching  may  be  done  in  a phase  stable  way  which  adds  only 
about  30  cm  to  the  length  of  the  accelerator  section 
( ~ 4 m)  . 
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SECTION  V 


CYCLOTRON  WAVE  GROWTH 

The  cyclotron  wave  growth  section  is  one  of  the 
most  critical  components  of  the  accelerator  experiment. 

Here  the  correct  negative  energy  cyclotron  mode  must  be 
properly  excited  and  phase  locked,  and  then  grown  to  the 
appropriate  amplitude.  Care  must  be  taken  so  that  non- 
axisymmetric  modes  are  not  grown  to  disruptive  amplitudes 
and  that  nonlinear  saturation  phenomena  do  not  limit  the 
amplitude  of  the  axisymmetric  wave  below  that  desired. 

A . RF  Exciter 

The  method  of  exciting  and  phase  locking  the 
desired  cyclotron  wave  has  been  examined.  It  appears  that 
a current  loop  drive  will  efficiently  couple  to  the  desired 
cyclotron  mode  if  located  at  0.7  to  0.8  of  the  wall  radius, 
as  indicated  in  Appendix  K.  Location  of  the  loop  too 
close  to  the  wall  results  in  an  undesirable  drop  in  loop 
impedance  due  to  the  induced  shorting  currents  in  the  wall. 


T-.-g---  ".- 


This  method  of  excitation  should  result  in  a 
fairly  pure  axisymmetric  standing  cyclotron  wave  (equal 
mix  of  positive  and  negative  energy  contributions)  with 
small  nonaxisymmetric  components.  Subsequent  growth 
sections  will  amplify  the  negative  energy  component  and 
damp  the  positive  energy  component. 


B . Wave  Growth  Configuration 

A continuing  examination  of  passive  external  liner 
configurations  capable  of  growing  the  negative  energy 
cyclotron  wave  has  been  maintained. 

A generally  accepted  method  reported  previously 
(Ref.  1)  is  through  use  of  a resistive  liner  configuration 
with  azimuthal  shorting  rings.  This  edge  ring  liner  has 
the  disadvantage  of  fairly  long  growth  lengths  (which 
necessitates  the  5 meter  growth  section  in  the  conceptual 
experimental  design)  and  attendant  kink  instabilities, 
which  one  must  take  care  to  avoid. 

Continuing  work  on  the  use  of  helical  slow  wave 
structures  indicates  that  these  disadvantages  may  be  over- 
come. The  first  helical  structure  examined,  the  sheath 
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helix  configuration,  yields  very  short  growth  lengths 
(~15  cm  for  the  model  accelerator  parameters  given  in 
Section  I);  but,  unfortunately , a detailed  examination  of 
this  structure  over  a broad  range  of  operating  parameters 
shows  that  the  m = 1 kink  mode  always  grows  at  least 
slightly  faster  than  the  desired  axisymmetric  mode.  As 
such,  such  a structure  could  be  utilized  only  if  the 
axisymmetric  mode  were  sufficiently  preferentially  ex- 
cited. This  analysis  is  presented  in  Appendix  L. 

In  Appendix  M,  a sufficient  condition  is  derived 

for  the  suppression  of  the  fast  kinking  instability  of  the 

sheath  helix?  namely,  u)2  < n2,  where  u>  , n are  the 

P P 

relativistic  plasma  and  cyclotron  frequency  of  the  beam, 
respectively. 

A second  helical  structure,  the  so-called  tape 
helix,  appears  to  retain  the  good  growth  features  of  the 
sheath  helix  for  the  axisymmetric  mode,  while  allowing  the 
growth  rate  for  the  m = 1 mode  to  be  reduced.  This  arises 
from  the  fact  that  in  the  tape  helix,  the  finite  spacings 
betwean  the  windings  give  rise  to  stop  band  characteristics 
By  locating  the  frequency  and  wavelength  of  the 


; 


I 
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nonaxisymmetric  mode  near  these  stop  band  regions,  a 
reduction  in  growth  rate  appears  feasible.  Details  of  this 
analysis  are  given  in  Appendix  N. 

A fourth  interesting  wave  growth  configuration 
utilizes  discrete  azimuthal  current  loops  connected  to 
external  passive  RC  networks.  By  choosing  the  value  of  C 
(capacitance)  so  that  the  networks  are  resonant  with  the 
loop  inductance  at  the  driving  frequency  of  240  MHz,  very 
small  growth  lengths  for  the  axisymmetric  mode  can  be 
obtained.  Furthermore,  it  appears  (although  it  has  not 
been  conclusively  shown)  that  the  m = 1 kink  modes  can  be 
suppressed  by  bringing  the  walls  of  the  waveguide  suffi- 
ciently close  to  the  driving  loops.  (If  the  walls  are 
brought  too  close,  the  waveguide  must  be  slotted  longitu- 
dinally to  prevent  shorting  of  the  loops.)  A discussion  of 
this  method  of  wave  growth  and  some  of  the  attendant 
problems  is  given  in  Appendix  N. 

C . Detailed  Energy  Balance 

Of  considerable  importance  in  determining  total 
energetics  in  the  Auto-Rescnant  Accelerator  concept  is  an 
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understanding  of  the  detailed  exchange  of  energy  during  the 
growth  of  the  wave.  This  question  is  examined  in  detail 
in  Appendix  P.  Here  it  is  shown  that  energy  is  extracted 
from  the  beam  not  at  the  center  of  the  beam,  but  near  the 
outer  edges.  This  is  important  when  considering  the  total 
energy  source  available  for  growth  and  acceleration; 
particles  near  the  center  of  the  beam  have  relatively  less 
free  energy  available  than  particles  near  the  edge  of  the 
beam  due  to  the  virtual  cathode  they  must  traverse  when 
leaving  the  anode  foil.  If  the  energy  for  growth  of  the 
cyclotron  wave  came  predominately  from  these  particles, 
the  wave  might  of  necessity  be  limited  to  low  amplitude. 

D . Wave  Growth  Saturation  and 
Other  Nonlinear  Effects 

As  the  cyclotron  mode  amplitude  grows  in  magnitude, 
nonlinearities  increase  in  importance.  These  nonlinearities, 
in  addition  to  modifying  the  linear  dispersive  properties 
of  the  mode  (nonlinear  frequency  shifts) , will  influence 
mode  growth,  and  may  introduce  a variety  of  new  phenomena: 
resonant  decay  instabilities,  nonlinear  Landau  damping,  wave 
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breaking  (or  particle  crossing) , particle  heating,  and 
"resonance  broadening."  The  beam  cyclotron  modes  do  not 
appear  to  be  subject  to  resonant  decay  instabilities  (see 
Appendix  C,  Ref.  1);  and  for  a sufficiently  cold  electron 
beam,  nonlinear  Landau  effects  can  be  neglected. 

The  task  of  evaluating  the  remaining  effects  of 
nonlinearities,  understanding  their  relative  importance, 
and  estimating  their  effects,  has  been  undertaken,  with 
a major  emphasis  in  determining  the  exact  nonlinear  mechanism 
(or  mechanisms)  responsible  for  saturation  of  the  growth  of 
the  cyclotron,  wave.  This  question  is  not  yet  completely 
resolved. 

Several  analyses  have  been  undertaken  to  examine 
the  general  question  of  nonlinearities  and  to  determine  at 
what  field  amplitudes  saturation  should  occur.  One  such 
nonlinear  effect  is  discussed  in  Appendix  Q,  where  attention 
is  focussed  on  the  nonlinear  shift  Aw  in  the  wave  frequency. 
The  magnitude  of  Aid  is  a measure  of  the  nonlinear  correc- 
tions to  the  linear  dispersion  relation.  It  is  of  interest 
to  determine  the  magnitude  of  Aw  in  terms  of  the  mode 
amplitude,  since  such  changes  may  affect  optimum  growth  of 


beam  cyclotron  modes  by  resistive  dissipation,  or,  in  the 
case  of  growth  produced  by  resonant  coupling  to  other  modes 
of  the  system,  may  lead  to  mode  saturation  when  Au>  is  of 
the  order  of  the  growth  rate.  This  investigation  is  not 
yet  finished.  Final  conclusion  on  the  importance  of  Auj 
must  await  completion  of  the  analysis. 

Although  the  understanding  of  the  nonlinearities 
is  in  its  early  stages,  it  is  already  possible  to  obtain 
rough  estimates  of  the  wave  amplitude  at  which  nonlinear 
effects  on  the  electrons  first  become  significant.  A 
number  of  such  effects  are  evaluated  in  Appendix  R subject 
to  several  approximations  which  simplify  the  cyclotron 
eigenfunctions . 

This  analysis  primarily  examines  a number  of  ways 
in  which  linear  approximations  lose  their  validity,  such 
as  when  fluctuations  in  density  or  electric  field  amplitudes 
become  comparable  to  the  ambient  values,  when  finite  electron 
gyroradius  effects  become  significant,  or  when  the  radial 
electron  flow  becomes  nonlaminar  due  to  radial  electron 
crossings  or  overtaking.  Although  all  of  these  effects 
are  related  to  one  another,  the  most  sensitive  nonlinear 
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effect  of  this  variety  is  the  variation  in  the  electron 


cyclotron  resonance  |w  - k c 6 + ft|,  sometimes  referred 

z z 

to  as  "resonance  broadening,"  which  occurs  due  to  the  in- 
fluence of  the  waves  on  the  electron  quantities  6,  and 
y.  Appendix  R also  contains  some  discussion  of  physical 
consequences  which  could  arise  from  nonlinear  wave  ampli- 
tudes of  this  level,  such  as  wall  collisions  of  the  outer 
electrons  or  field  emission  of  electrons  from  the  conductor 
walls.  It  appears  that  such  effects  can  be  controlled 
through  the  presence  of  the  vacuum  region  between  the  beam 

I 

and  the  wall,  and  through  limiting  the  ion  density  to  a 

sufficiently  small  fraction  of  the  electron  density. 

. 

In  order  to  guide  the  analytical  effort  towards 
the  relevant  nonlinear  saturation  mechanism,  computer  sim- 
ulation  of  the  growth  and  saturation  of  the  axisymmetric 
cyclotron  wave  in  an  anisotropic  resistive  background  has 
been  undertaken.  Details  are  given  in  Appendix  S. 

Specifically,  the  cyclotron  wave  is  grown  from  noise  with 
the  beam  initialized  in  an  equilibrium  background  state. 

The  homogeneous,  anisotropic  resistance  is  oriented  in  the 
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streaming  direction  (zero  conductivity  in  the  transverse 
directions) . The  computational  results  show  a growth  rate 
of  approximately  twice  the  growth  rate  predicted  by  linear 
theory  for  this  configuration  (Appendix  T) . Because  some 
of  the  assumptions  made  in  the  theoretical  derivation  are 
not  well  satisfied  in  the  simulation,  this  constitutes 
satisfactory  agreement.  The  computer  simulation  shows  that 
the  instability  tends  to  saturate  when  the  perturbed  electric 
field  in  the  streaming  direction  reaches  10  to  20  percent  of 
the  transverse  DC  electric  field.  Although  this  is  the 
approximate  level  where  one  would  expect  electron  particle 
crossing  to  occur,  one  cannot  justifiably  infer  that  particle 
crossing  is  the  responsible  mechanism,  since  other  nonlinear 
effects  discussed  in  Appendix  R are  also  important  at  these 
field  amplitudes. 

E . Wave  Growth  Diagnostics 

Experimentally,  it  is  necessary  to  monitor  the 
wave  amplitude  and  phase  velocity  of  the  cyclotron  wave  as 
it  progresses  down  the  wave  growth  section.  The  basic 
diagnostics  for  detecting  the  wave  would  be  small  electric 


and  magnetic  probes.  Since  the  wave  electromagnetic  fields 
extend  beyond  the  beam  to  the  conducting  walls  of  the  vessel, 
the  probes  can  be  placed  outside  the  beam,  thus  preventing 
interference  with  the  beam  and  obviating  probe  damage. 
Further,  since  the  coupling  is  electromagnetic,  they  can 
be  absolutely  calibrated  in  a standard  transmission  line  for 
field  amplitude  as  a function  of  frequency.  The  electric 
probes  measure  the  radial  electric  field  at  the  wall  and 
are  effectively  small  patches  isolated  from  the  wall  upon 
which  induced  charge  can  be  measured.  They  are  quite  sen- 
sitive; their  frequency  response  can  extend  from  the  100  MHz 
of  the  lowest  waves  of  interest  to  several  GHz. 

The  magnetic  probes  are  simply  one  or  two  turns 
of  fine  wire  in  a loop  of  diameter  1 to  4 mm,  appropriately 
shielded  electrostatically  and  connected  to  cable.  They 
are  placed  just  inside  the  conducting  liner  and  aligned  to 
measure  either  B or  Bn . Exterior  to  the  beam,  B.  is 
small  for  the  axisymmetric  mode,  but  not  for  the  non- 
axisymmetric  kink  and  higher  order  modes.  Hence,  discrimina- 
tion between  the  desired  axisymmetric  and  the  undesired 
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nonaxisymmetric  cyclotron  modes  is  available  through  the 
B measurement. 

For  these  cyclotron  waves,  magnetic  probe 
sensitivity  is  comparable  with  that  of  electric  probes. 
However,  their  frequency  response  is  somewhat  restricted, 
with  the  most  sensitive  probes  peaking  at  200  MHz.  With 
some  sacrifice  in  sensitivity,  suitable  coils  can  have  a 
response  extending  to  1 GHz. 

The  probes  would  be  arrayed  in  clusters  at 
intervals  1 m)  along  the  drift  tube.  Each  cluster  should 
include  electric  and  magnetic  probes  for  both  components 
of  B.  The  azimuthal  array  would  give  unambiguous  deter- 
mination of  structure  through  m = 2,  and  the  longitudinal 
array  determines  kz  for  wavelengths  of  1 to  40  centimeters 
Each  cluster  should  provide  complete  data  on  the  mode 
structure  (real  part  of  the  dispersion  relation) . Com- 
parison of  successive  clusters  down  the  tube  would  determine 
the  growth  rates  (imaginary  part  of  the  dispersion  relation) 
This  would  provide  all  the  essential  information  on  wave 
growth. 
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SECTION  VI 

ION  LOADING  AND  ACCELERATION 

A.  Methods  of  Ion  Loading 

Several  analyses  of  the  ion -wave  interaction  have 
been  performed  to  insure  that  the  proof -of -principle 
experiment  on  Auto-Resonant  acceleration  would  not  flounder 
on  unanticipated  difficulties  with  loading  the  ions  into 
the  cyclotron  wave,  or  with  retaining  the  ions  in  the 
wave  troughs  during  the  process  of  wave  acceleration.  The 
experimental  design  contemplates  a puff -gas  injector 
section  (see  Ref.  1,  Appendix  V),  wherein  the  subject  ions 
will  be  created  by  primary  beam  ionization  of  a neutral 
gas,  then  trapped  or  picked  up  by  the  preformed  travelling 
cyclotron  wave,  and  finally  carried  out  of  the  puff -gas 
section  into  the  accelerator  section  where  the  wave  and 
ions  are  accelerated  above  the  initial  wave  velocity. 

In  Ref.  1,  Appendix  V,  conditions  were  presented  for 
the  choice  of  neutral  gas  density  and  the  spatial  volume 
of  the  puff-gas  section,  such  that  cascading  via  secondary 
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ionization  could  be  avoided  and  the  accumulation  of 
untrapped,  residual  ions  in  the  puff -gas  section  would 
not  disrupt  the  electron  beam  during  its  pulse  time.  In 
Appendix  U of  this  report,  this  question  of  ion  loading  is 
further  explored  through  theoretical  and  computer  analysis 
of  a model  of  ion  pickup  and  acceleration.  It  is  shown 
that  ions  which  are  untrapped  following  wave  amplitude 
growth,  even  if  injected  into  the  accelerator  section  with 
a finite  velocity  v , may  not  later  be  picked  up  solely 
by  accelerating  the  wave  velocity  to  near  the  injection 
velocity  v . Instead  it  is  found  that  the  primary  re- 
quirement for  a high  pickup  efficiency  of  the  ions  into 
the  travelling  wave  is  that  the  wave  must  be  grown  to  a 
sufficient  amplitude — roughly  ei~  j Mv^2 — and  that  such 
growth  and  trapping  must  occur  before  the  wave  is  accel- 
erated. Provided  that  there  is  sufficient  amplitude 
growth  of  the  cyclotron  wave  before  the  puff-gas  section, 
there  is  no  harm  done  in  applying  a small  axial  field  to 
the  puff-gas  section  to  help  eject  the  residual  untrapped 
ions  and  so  retard  their  charge  accumulation. 
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B-  Stability  of  Acceleration 


It  has  been  shown  in  a number  of  analyses  (see  for 
example  Figure  U3  and  discussion  thereof)  that  once  the 
ions  are  sufficiently  trapped  in  the  wave,  the  wave  may  be 
accelerated  without  spilling  them  if  this  acceleration  is 
sufficiently  gentle.  However,  there  has  been  concern 
that  the  ions  might  be  spilled  during  acceleration  due  to 
spatial  ripples  or  nonadiabatic  variations  in  the  magnetic 
guide  field,  or  else  due  to  voltage  droop  and  ripple  in 
time  produced  by  the  impedance  variation  during  diode 
closure  or  by  transients  in  the  Marx  circuitry. 

The  first  question  of  stable  ion  acceleration  in 
the  presence  of  spatial  magnetic  field  variations  is 
examined  in  Appendix  I . Computer  calculations  are  de- 
scribed there  for  a code  which  integrates  the  one  dimen- 
sional ion  equation  of  motion  for  an  arbitrary  guide 
magnetic  field  profile.  A number  of  runs  were  performed 
for  the  parameters  of  the  proof -of -principle  experiment 


(i.e.,  y = 7,  B . ...  , = 2 5 kG , e $ =0.5  MeV , 

e zinitial  o 

ei  initial  = 0*25  MeV,  ~ 20  MeV,  initial  ion 
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phase  = - 0.75  radian  w.r.t.  bottom  of  well) . For  the 
cases  considered,  spatial  ripple  at  or  below  5 percent  of 
the  guide  magnetic  field  amplitude  presents  no  acceleration 
stability  problem.  Furthermore,  replacing  the  non- 
adiabatic  portion  of  the  theoretically  ideal  magnetic 
field  profile  with  an  easily  physically  realizable, 
adiabatic  profile  is  found  to  be  completely  satisfactory 
from  the  viewpoint  of  acceleration  stability,  and  results 
in  a total  increase  of  the  accelerator  length  of  less  than 
8 percent. 

The  problem  of  stable  ion  acceleration,  or  phase 
stability,  in  the  presence  of  high  frequency  time  ripples 
in  the  voltage  has  been  examined  previously  (see  Ref.  1, 
Appendixes  Q,  R) . Although  potentially  more  serious  than 
spatial  magnetic  field  rippling,  due  to  a cumulative  en- 
hancement effect  by  the  number  of  wavelengths  transited 
by  electrons,  this  high  frequency  voltage  rippling  was 
found  to  be  acceptable  if  the  ripple  amplitude  could  be 
held  below  ± 1/2  percent  for  frequencies  above  ~ 10  MHz. 
Subsequent  discussions  with  qualified  vendors  for  such 
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equipment  indicate  that  this  requirement  can  indeed  be 
met  for  discharge  into  a constant  100  ohm  load. 

For  the  case  of  a changing  diode  impedance,  the 
effects  of  voltage  droop  and  ripple  on  ion  phase  stability 
have  been  reevaluated  in  Appendix  B.  There  it  is  shown 
that  phase  stability  is  satisfactory  unless  the  voltage 
droop  over  the  200  ns  pulse  time  exceeds  20  percent,  or  low 
frequency  ripple  exceeds  5 to  6 percent — figures  Which  are 
well  above  the  values  anticipated. 

C . Ion  Diagnostics 

A conceptual  design  of  the  experimental  diagnostics 
required  to  verify  the  ion  parameters  has  also  been  con- 
sidered, and  there  appear  to  be  no  serious  difficulties. 

The  initial  requirement  of  measuring  the  output  ion 
current  intensity  may  be  accomplished  by  use  of  a large 
aperture  Faraday  cup  which  is  fairly  deep.  Use  of  a 
positive  bias  voltage  of  several  kilovolts  is  desirable  to 
prevent  any  secondary  electron  emission  from  escaping 
the  cup  and  giving  an  anomalously  high  reading. 
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Next,  rough  proton  energy  analysis  may  be 
accomplished  with  time  of  flight  techniques.  The  radial 
profile  of  the  proton  beam  current  can  be  estimated  with 
either  segmented  Faraday  cups  or  with  time  integrating 
dosimeter  measurements.  Time  resolved  measurements  of  the 
instantaneous  ion  beam  energy  dispersion  may  be  carried 
out  using  a magnetic  spectrometer. 


SECTION  VII 


OVERALL  SYSTEMS  CONSIDERATIONS 

Besides  the  specific  system  subcomponents  which 
have  been  treated  in  previous  sections,  consideration  has 
also  been  given  to  three  aspects  of  the  total  system. 
Specifically,  these  consist  of 

(1)  An  operational  control  system  logic  for 
running  of  the  experimental  accelerator. 

(2)  The  vacuum  requirements  necessary  in  the 
electron  beam  drift  tube  (waveguide) . 

(3)  Selection  of  a "dump"  for  the  electron 
(proton)  beam  which  minimizes  radiation 
production,  and  in  conjunction  with  this  the 
amount  of  shielding  necessary  to  protect 
personnel  working  in  the  general  area 
surrounding  the  accelerator. 

Each  of  these  items  is  considered  separately  below. 


For  purposes  of  control  system  design,  the 


Auto-Resonant  Accelerator  may  be  divided  in  five  logical 
subsystems: 

Vacuum  System 

Magnetic  Field  System 

Marx/Diode  Electron  Beam  Injector 

R.F.  Exciter 

Ion  Injector 

One  possible  flow  chart  of  the  system  sequencing  necessary 
for  an  experimental  run  is  shown  in  Figure  1.  A 
brief  description  of  the  logic  follows. 

Prior  to  any  run,  it  is  necessary  to  insure  that 
the  vacuum  level  in  the  beam-forming  diode  and  the 
waveguide  are  in  the  proper  pressure  range.  Once  this 
is  satisfied,  a general  visual  check  is  made  to  insure 
that  signal  cables  are  connected,  major  subsystems  are 
ready,  etc.  Then  the  main  power  line  to  the  control 
system  is  energized,  monitors  are  checked,  and  safety 
interlocks  are  set  after  insuring  chat  all  personnel  are 
outside  of  the  unshielded  area  around  the  accelerator.  A 


Figure  1.  Accelerator  Operation  Flow  Chart 


warning  bell  is  then  sounded,  and  main  power  lines  energized 
to  the  various  subsystem  high  voltage  power  supplies.  A 
brief  waiting  period  has  been  allowed  at  this  point; 
shot  run  sheets,  etc.  can  be  prepared  at  this  time.  At 
the  start  of  the  experimental  shot,  flashing  warning  lights 
are  turned  on  and  safety  dumps  on  capacitor  banks  are  lifted. 
The  Marx  bank,  magnet  field  bank,  and  r.f.  exciter  supply 
are  then  charged  to  the  appropriate  voltage.  A brief 
holding  period  is  then  provided  to  arm  all  data  recorders 
(oscilloscopes,  etc.)  followed  by  initiation  of  the  main 
discharge  trigger  signal. 

The  discharge  sequence  is  primarily  involved  in 
controlling  the  relative  timing  of  the  various  system 
subcomponents;  the  primary  delay  is  in  the  roughly  5 

millisecond  ring-up  time  of  the  magnetic  field.  Follow- 
ing each  shot,  capacitor  bank  dumps  are  enabled  and  an 
arbitrary  length  of  time  allowed  before  the  next  shot  to 
analyze  the  data  obtained.  Sequencing  involved  with 
any  interlock  interrupt,  or  completion  of  a complete 
experimental  run,  is  also  indicated. 


The  control  system  should  be  automated  as  much  as 
possible  to  minimize  operational  requirements  imposed 
upon  the  machine  operator,  and  consequently  maximizing 
the  number  of  useful  shots. 


The  ambient  pressure  level  necessary  in  the  beam 
drift  tube  region  is  determined  by  the  requirement  that 
the  amount  of  space  charge  neutralization  of  the  electron 
beam  caused  by  ionization  of  the  ambient  gas  not  greatly 
perturb  the  rigid  rotor  equilibrium  of  the  beam.  The 
details  of  the  calculation  can  be  found  in  Appendix  V; 
the  result  is  that  a pressure  in  the  10"6  to  10~5  range 
is  necessary. 

To  achieve  this  level  with  a modest  sized 
diffusion  pump  (with  cold  trapping) , the  waveguide  should 
be  fabricated  so  that  it  is  "clean"  from  a vacuum  point 
of  view.  Thus,  the  waveguide  should  be  fabricated  from 
type  304  stainless  steel  and  electro-polished.  Use  of 
conflat  vacuum  seals  will  allow  waveguide  bake-out  to 
minimize  the  time  required  for  water  vapor  pump-out. 
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Given  these  general  specifications,  attainment  of  the 
necessary  vacuum  level  is  not  a problem. 

C.  Beam  Dump  and  Radiation  Shielding 

Selection  of  the  best  material  for  stopping  the 
electron  (and/or  proton)  beam  at  the  end  of  the  accelera- 
tor is  predicated  on  withstanding  the  high  energy  fluence 
in  the  electron  beam  and  minimizing  radiation  production. 
As  shown  in  Appendix  W,  graphite  is  the  best  selection. 

Based  on  obtaining  data  at  the  rate  of  one  shot 
every  5 minutes,  Appendix  W also  indicates  that  4 
feet  of  concrete  placed  around  the  beam  dump  is  necessary 
to  reduce  radiation  levels  to  national  safety  levels. 

In  addition,  a wall  of  concrete  2 feet-  thick  along 
both  sides  of  the  accelerator  is  required  to  prevent  high 
radiation  exposure  on  any  shot  in  which,  for  whatever 
reason,  the  beam  strikes  the  waveguide  wall. 

i 

To  insure  proper  monitoring  of  radiation  exposure, 
all  laboratory  personnel  would  be  required  to  wear  a 
radiation  film  badge  and  have  the  exposure  determined 
on  a monthly  basis.  Pocket  dosimeters  placed  on  the 
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ELECTRON  BEAM  INJECTOR  CONSIDERATIONS 


Introduction 

In  designing  an  experimental  program  to  demonstrate 
the  feasibility  of  the  Auto-Resonant  Accelerator  concept, 
the  following  parameters  have  been  chosen  for  the  out- 
put ion  (proton)  beam: 

Proton  Energy  =30  MeV 
Proton  Current  =30  amps 
Pulse  Duration  ;>  100  ns 


The  peak  power  in  the  ion  beam  is  then  9 x 103  watts, 
and  this  should  represent  a small  fraction  of  the  acceler- 
ating media  (electron  beam)  peak  power  level.  Thus,  we 
require 

9 x 10s  « V I 

e e- 

To  be  conservative,  the  following  values  were  chosen: 
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V =3.0  MeV 
e 

I = 30  kA 
e 

t = 200  ns  pulse  width 
e 

ion  beam  power  is  then  1%  of  that  available  in  the 
electron  beam. 

In  the  remainder  of  this  appendix  consideration  is 
given  to  specifying  the  electron  beam  injector,  based 
upon  these  nominal  characteristics . 

Time  Varying  Diode  Impedance 

For  real  cold  cathode  diode  systems,  the  electron 

source  is  a plasma  sheath  formed  on  the  cathode  surface. 

This  plasma  has  a high  thermal  temperature  and  expands 

ballistically  with  a velocity  of  2 to  3 cm/ys.  This 

produces  a continuously  decreasing  gap  spacing  in  the 

diode,  and  causes  the  effective  diode  impedance  (Z  ) to 

L 

decrease  with  time. 

If  one  assumes  that  the  electrical  network  feeding 
the  diode  can  be  represented  by  a battery  and  an  effective 
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internal  impedance  (Z^) , then  the  load  voltage  will 
scale  as 


V ~ " 

L ZL  + Zr 


The  change  in  load  voltage  introduced  by  the  time 


dependence  of  Z is  then 
* L 


avt  zt  az, 

— “ = ( — - — ) — - 
VL  UL  + ZI  ' ZL 


while  the  change  in  load  current  is  given  by 


UA  U 

— = - ( — — ) 

I,  \ z.  +z.  / 


L I 


Since  the  local  phase  velocity  of  the  accelerating  wave 


v = ~ V 

* 'ju  + n l 

o 


then  the  change  in  local  phase  velocity  is 
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To  insure  that  the  ion  energy  does  not  change  appreciably 


throughout  the  pulse  length.  Equation  (A5)  (Equation  written 
as  Eq.  hereafter)  should  be  as  small  as  possible.  From  Eq. 
(A2)  this  implies  that  the  electron  accelerator  internal 
impedance  be  small  compared  to  100  ohms.  Thus,  in  evalu- 
ating a machine  design  we  want 

Z < Z (=-  100  ohms)  (A 6) 

I L 


t 


High  Frequency  Ripple 

As  discussed  (Ref,  1,  Appendixes  Q and  R),  high 
frequency  voltage  ripple  (superimposed  on  a slow  voltage 
droop)  can  produce  phase  shifts  which  de-trap  ions  and 
terminate  the  collective  acceleration  process . Filter 
elements  are  therefore  probably  required  between  the 
high  voltage  source  and  the  electron  beam  diode  load. 
Also,  since  the  magnitude  of  voltage  ripple  generated 
by  inductive  components  in  the  system  scales  like  dl/dt, 
a slow  pulse  rise  time  is  desirable. 
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Both  analytic  and  computer  studies  of  the  effect 


of  high  frequency  ripple  indicate  that  in  the  frequency 
range  from  10  MHz  up  the  amplitude  should  be  kept  at 

0.5%  or  less.  This  has  been  incorporated  in  the  injector 
specifications . 


Electron  Injector  Specification 

The  following  specifications  have  tentatively  been 
chosen  for  the  electron  beam  injection  source: 

1.  Load:  100  Q (constant  impedance  dummy  load 

for  tests) 

2.  Peak  voltage  (V  ) into  100  Q:  3 MeV 

pk 

3.  Duration:  200  ns  minimum  useful  pulse  between 


0.96  V , points  (±  2%) 
pk 

4.  Reproducibility  of  V ^ shot  to  shot:  ± 2% 

5.  Rise  time:  s 100  ns  (5%  to  95%  of  V J 

pk 

6.  Prepulse  suppression:  Maximum  prepulse  less 

than  3%  of  V , (measured  with  no  dummy  load) 
pk 

7.  Ripple:  Less  than  ± 1/2%  high  frequency 

(2  10  MHz)  ripple  during  useful  pulse 


(V  2 0.96  V , ) 
pk 
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8.  Triggering:  Command  triggered  with  less  than 

± 200  ns  jitter  between  triggering  and  firing 

9.  Pulse  clipping  crowbar  switch  to  reduce  energy 
into  100  n load  to  less  than  10%  of  initial 
stored  energy  after  200  ns  useful  pulse  (0  to 
+ 50  ns  jitter  allowed  in  crowbar  switch) 

10.  The  electrical  circuit  driving  the  load  should 
have  a low  internal  impedance. 

The  prepulse  suppression  criterion  was  added  to 
eliminate  cathode  turn-on  prior  to  arrival  of  the  main 
voltage  pulse.  The  slow  rise  time  requirement  will  hope- 
fully minimize  ripple  generated  by  inductive  components . 
The  addition  of  a crowbar  to  eliminate  any  late -time 
energy  flow  in  the  diode  is  meant  to  protect  the  diode 
anode.  This  is  discussed  in  another  appendix. 

Finally,  it  is  worth  noting  that  the  testing  of  the 
injector  is  done  using  a 100  ohm  dummy  load.  This  was 
felt  to  be  necessary  since  at  this  time  the  actual  diode 
has  not  been  designed,  and  therefore  its  time  behavior 
cannot  be  predicted.  However,  it  is  anticipated  that 


Design  of  the  conventional  diode  which  will  be  used 
for  initial  experimentation,  as  well  as  foilless  diode 
design,  will  be  made  using  a computer  code  (Ref.  2 ) 
similar  to  that  developed  at  Sandia  (Ref.  3).  This 
self-consistent  code  is  relativistically  correct,  and 
includes  the  beams  self-magnetic  field  as  well  as  an 
externally  applied  field. 


2.  "Research  on  Beam  Stopping  Due  to  Self-Magnetic 
Compression,"  DASA  2167,  Ion  Physics  Corp.,  Nov.  1968. 

3.  Boers,  J.  E.,  "Digital  Computer  Simulation  of  High 
Current,  Relativistic,  and  Field  Emission  Electron  Tubes," 
Proc.  of  11th  Symp.  on  Electron,  Ion,  and  Laser  Beam 
Technology,  1972. 
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APPENDIX  B 


THE  EFFECTS  OF  VOLTAGE  DROOP  AND  LOW  FREQUENCY  RIPPLE 
ON  ACCELERATOR  PHASE  STABILITY 

Introduction 

Specifications  of  the  Marx  generator  and  diode 
assembly  electron  beam  source  detailed  in  Appendix  B were 
quoted  for  output  into  a constant  100  ohm  load.  Subse- 
quent discussions  with  qualified  vendors  of  such  assem- 
blies showed  that  the  ripple  suppression  circuits 
envisioned  are  tuned  to  this  constant  100  ohm  load 
impedance  to  provide  the  proper  low  frequency  voltage 
waveshape . 

An  actual  diode  will  of  course  undergo  closure 
during  the  200  ns  pulse  time  at  a rate  of  roughly  1 to  2 
cm/(js,  resulting  in  estimated  40  percent  droop  in  load 
impedance;  i.e.  from  100  ohms  to  60  ohms.  The  question 
of  the  effect  of  such  a changing  load  impedance  was 
subsequently  brought  to  the  attention  of  the  vendors. 

Two  general  conclusions  were  deduced: 

ii 
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1.  Such  a changing  load  impedance  will  have  a 


negligible  effect  on  the  ability  of  the  circuits  to 
effectively  suppress  high  frequency  ripple. 

2.  Such  a changing  impedance  will,  however, 
affect  significantly  the  overall  output  voltage  wave- 
form, resulting  in  a significantly  increased  voltage 
droop  and  the  introduction  of  low  frequency  ripple  com- 
ponents associated  with  the  pulse  forming  networks . 

In  order  to  ascertain  the  effect  of  such  droop  and 
low  frequency  ripple,  an  updated  version  of  the  phase  sta- 
bility code  was  run  with  a voltage  wave  form  modelling 
increasingly  larger  amounts  of  droop  and  low  frequency 
ripple.  Specifically,  the  low  frequency  ripple  was 
assumed  to  correspond  to  roughly  8 MHz  at  levels  of  1,  2, 

3,  4,  6,  8,  and  10  percent  with  voltage  droop  of  10,  20, 

30,  40,  and  50  percent  over  the  200  ns  pulse  time. 

Accelerator  parameters  were  chosen  for  an  output 
of  30  MeV  ions  from  an  initial  injection  energy  of  250  keV, 
using  a cyclotron  wave  amplitude  of  500  keV.  Two  specific 
accelerator  configurations  were  considered. 
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In  the  first,  a short  accelerator  with  a low  safety 
factor  of  only  1.5  (high  fraction  of  total  electric  field 
available  used  for  acceleration),  was  considered.  Results 
show  that  phase  stability  was  maintained  for  voltage 
droops  of  20%  or  less  or  for  ripple  level  at  6%  or 
lower . 

For  the  second,  more  conservatively  designed 
accelerator  with  a safety  factor  of  2.0  (i.e.  utilizing 
only  half  of  the  total  electric  field  for  acceleration) , 
phase  stability  was  again  maintained  at  voltage  droops 
of  20  percent  or  less  but  the  ripple  amplitude  was  here 
required  to  be  less  than  the  6 percent  figure  for  the 
first  case. 

A tabulation  of  the  results  of  this  series  of  runs 


is  given  in  Table  Bl. 


TABLE  Bl 


THE  EFFECTS  OF 
RIPPLE  ON 

VOLTAGE  DROOP  AND  LOW  FREQUENCY 
ACCELERATOR  PHASE  STABILITY 

Safety 

Factor 

Droop 

(Percent) 

Ripple  @ 8 MHz 
(Percent) 

Acceleration 
Time  (ns) 

1.5 

10% 

1% 

70 

1.5 

10 

2 

93 

1.5 

10 

3 

98 

1.5 

10 

4 

101 

1.5 

10 

6 

106 

1.5 

10 

8 

phase  stability  lost 

1.5 

10 

10 

phase  stability  lost 

1.5 

20 

1 

106 

1.5 

30 

1 

phase  stability  lost 

1.5 

40 

1 

phase  atability  lost 

1.5 

50 

1 

phase  stability  lost 

2 .0 

10 

1 

100+ 

2.0 

10 

2 

105 

2.0 

10 

3 

106 

2 .0 

10 

4 

107 

o 

• 

CM 

10 

6 

phase  stability  lost 

o 

• 

CM 

10 

8 

phase  stability  lost 

2.0 

10 

10 

phase  stability  lost 

2.0 

20 

1 

132 

2 .0 

30 

1 

phase  stability  lost 

2 .0 

40 

1 

phase  stability  lost 

2.0 

50 

1 

phase  stability  lost 

APPENDIX  C 

ANODE  SURVIVAL  SPECIFICATION 


Using  a real  diode,  that  is  a cold  plasma  cathode, 
the  current  density  as  a function  of  time  can  be  written 


j(t) 


f (V(t)) 
d*  (1  - t/T)2 


with  f (V)  the  functional  dependence  of  current  density 
on  diode  voltage  and  d^  the  initial  gap  spacing; 

d 

_ _o 
v 

o 

where 

vq  = plasma  closure  velocity 
« 2.5  x 10“ 3 cm/ns 

For  voltages  greater  than  200  keV,  and  assuming  a 0.001  inch 
titanium  anode,  the  dose  delivered  to  the  anode  is 


im  gat 


where 


dE/dx  = ionization  energy  loss  by  the  primary 
beam  in  volt.  cm2/gm 


If  time  is  measured  in  nanoseconds,  then  in  this  range  one 
can  use  the  approximation 


f V dE  , , ^ Q 0.693  a,A 

■ - — (v)  * lO-9  sec  = — -r — v1'274 
d2  dx  d2  MV 

o o 


where  V is  measured  in  units  of  10s  volts . 

For  energies  of  70  keV  and  lower,  the  electron  beam 
is  totally  stopped  in  the  anode  and  the  dose  delivered  is 
higher  than  that  given  by  the  above  expression.  Between 
70  keV  and  200  keV  there  is  a smooth  transition  to  the 


above  formula . 


At  time  t = t , the  diode  voltage  pulse  is  over 
and  no  more  energy  is  delivered  to  the  foil.  Thus, 
neglecting  the  enhanced  dose  below  200  keV,  the  foil  dose 
is  roughly 


0.693 

d2 


(1  - t/x)2 


For  titanium,  the  dose  required  to  raise  it  to  the  melt 
temperature  is  « 200  cal/gm.  Thus,  it  is  required  that 


:: 
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APPENDIX  D 
DIODE  DIAGNOSTICS 

A knowledge  of  the  injected  beam  parameters — energy 
and  current — is  essential  for  interpreting  experimental 
measurements  made  in  the  development  of  the  Auto-Resonant 
Accelerator.  While  the  manufacturer  of  the  pulsed  accelera- 
tor used  to  form  the  electron  beam  will  provide  diagnostics 
to  measure  the  capacitively  graded  tube  voltage,  inductive 
corrections  of  relatively  large  magnitude  must  be  made  to 
infer  the  actual  diode  gap  voltage  from  this  type  of  data. 

To  eliminate  this , the  diagnostic  configuration  shown  in 
Figure  Dl  is  deemed  necessary  for  diode  characterization. 

The  probes  consist  of  a capacitive  voltage  divider  to 
measure  the  shank  voltage  very  near  the  cathode,  and  a 
resistive  current  shunt  to  measure  the  wall  return  current 
driven  by  the  B.  field  generated  by  the  beam  current. 

The  resistive  current  shunt  consists  of  a dielectric 
break  in  the  wall  return  path  shunted  by  a low  impedance 
band  of  parallel  (carbon  composition)  resistors.  This 
provides  a low  impedance  voltage  source  for  driving  a 50  ohm 
signal  cable.  The  resistor  band  structure  has  an 


i 
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I 


intrinsically  low  self-inductance,  and  a subnanosecond 
rise  time  is  obtained.  This  is  more  than  adequate  for  the 
anticipated  slow  current  rise  time  of  the  pulsed  accelera- 
tor. Calibration  of  the  shunt  can  be  accomplished  using  the 
50  ohm  transmission  line  structure,  easily  fabricated  out  of 
sheet  metal,  shown  in  Figure  D2.  A transmission  line  pulser 
launches  a square  voltage  wave  onto  a 50  ohm  cable;  the 
cable  feeds  through  a GR  connector  onto  a 50  ohm  biconic 
line  which  matches  into  a 50  ohm  coaxial  line.  The  current 
shunt  is  located  in  the  middle  of  this  line.  A second  bi- 
conic line/GR  connector  feeds  the  output  cable  (50  ohm) 
which  goes  to  a matched  termination  oscilloscope  preamp 
input.  The  continuity  of  the  50  ohm  wave  impedance  along 
this  structure  ensures  that  no  reflections  will  occur,  and 
the  measured  voltage  amplitude  at  the  oscilloscope  (V2)  is 
actually  the  voltage  driving  current  in  the  shunt.  Thus, 
the  shunt  current  is  simply 


I 

s 


(Dl) 


A second  oscilloscope  measures  the  output  voltage  of  the 

current  shunt  (V^) . The  resistive  shunt  calibration  is 
then  obtained: 


' 

l 

I 


j 

j 
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amps 

volt 


(D2) 


;> 


Gi  = ^T 


= inverse  of  shunt  resistance 


It  should  also  be  noted  that  the  transmission  line  struc- 
ture can  be  used  to  calibrate  B and  B loops , as  well 
as  E-field  probes. 

The  cathode  voltage  monitor  depends  on  the  capaci- 
tive division  of  the  total  voltage  (V  (t) ) between  the 

o 

shank/monitor  band  capacitance  (C^)  and  the  monitor  band/ 

ground  capacitance  (C2) . A resistor  (R)  feeds  the  band 

signal  (V  (t) ) onto  a coaxial  transmission  line  of  wave 

impedance  Z . The  relationship  between  V and  V is 

o M 

given  by: 

, t 

Vo(t>  = Gv  {VM(t)  + 7 I VH(t')  dt'}  ,D31 

o 

where 


(R  +Z)  (C  +C  ) 

G = — 

v ZC. 


(D4) 


and 


(R  + Z)  (C1  +C2) 


(D5) 
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For  reference,  the  inverse  relation  is 


v«(t> 


t 


I 


dV  (t') 
O 

dt 


t '/t 


dt 


(D6) 


From  Eq.  (D3)  it  can  be  seen  that  over  most  of  the  pulse  the 
integral  correction  can  be  made  small  if  the  time  constant, 
t,  is  large  compared  to  the  pulse  length.  The  major 
correction  will  be  at  late  time,  low  voltage  levels  when 
the  correction  can  be  large.  This  will  not  be  a problem 
during  the  major  flat  top  portion  of  the  pulse  when  the 
main  Auto-Resonant  acceleration  is  occurring. 

Calibration  of  the  monitor  requires  a determina- 
tion of  G and  t-  This  can  be  done  in  situ  using  the 
v 

technique  shown  in  Figure  03.  A transmission  line  pulser 

generates  a square  pulse  which  is  split  by  an  unmatched 

tee.  One  leg  of  the  tee  goes  directly  to  the  cathode  shank, 

and  quickly  charges  up  the  tube/diode  stray  capacitance  to 

a voltage  V . The  other  leg  goes  to  an  oscilloscope 
max 

which  records  V . The  voltage  probe  output,  V,  is 
max  M 

recorded  on  another  oscilloscope.  From  Eq.  (06)  , 


assuming  V is  a step  function,  there  obtains 


SITU  CALIBRATION  OF 
)E  VOLTAGE  MONITOR 


The  slope  of  the  decay  of  V (t)  determines  t,  while 

M 

the  gain  calibration  is  given  by 
V 

, !S2_ 

v VM(t-0) 

A cross-check  on  the  gain  calibration  can  be  obtained  by 
noting  that 


G 

v 


(R  + Z)  (Cx  +C2) 


Since  t is  measured  independently  and  Z is  known,  it 
is  only  necessary  to  calculate  to  provide  another 

determination  of  . Since  the  shank/monitor  geometry 
is  that  of  a coaxial  line,  then  to  a good  approximation 


C 


1 


2rr  z 

o_ 

3m  r /r 

M s 


(D 


where 


r = monitor  radius 
M 


r = shank  radius 
s 

L = length  of  monitor  band 

I ' 

Previous  experience  has  shown  that  these  two  methods  of 
determining  give  good  agreement. 


' 


^1 


APPENDIX  E 


TRANSITION  TO  BEAM  EQUILIBRIUM  REGION 


Since  the  electron  beam  is  injected  into  the  drift 
tube  region  through  a conducting  anode  foil  which  is  at 
the  same  (ground)  potential  as  the  drift  tube  (waveguide) 
wall,  there  is  a transition  region  near  the  anode  where 
the  radial  electric  fields  are  small  and  the  axial  electric 
field  dominates.  In  what  follows  the  extent  of  the  transi- 
tion region  is  estimated,  and  comments  are  made  on  the 
beam  equilibrium. 

Assume  that  an  equilibrium  exists,  and  that  the 


beam  density  is  uniform  (=  nQ)  at  all  points.  The  poten- 


tial  distribution  is 

defined 

by 

R 1©  \ + 

52  a _ 

-ju2  a2 
( po  ') 

7 * ~ R SR  l 

K ) T 

dR  J 

3Z2 

\ c2  ) 

(El) 


where 


1 


= inner  radius  of  drift  tube  (also  assumed  equal 
to  the  electron  beam  radius) 


m c2/e  = 511  keV 
o 


Setting 


cp(R,Z)  = cp1  (R)  + cd2(R,Z) 


w2  a2 
po 


31  = ( 


7acp_  = 0 


with  the  boundary  conditions 


cp(R  = a,  Z)  = cp(R,  Z = 0)  = 0 


then  it  is  easy  to  show  that 


opL  (R)  = 


- 1/4  ( 


a)2  a2 
po  \ 

“V"  ) Ve  (1 


CD. 


(R,Z)  = 1/4  ^ 


)\ 


n=l 


JiV 


J (X  R) 
o n 


X 


-X  Z 
n 


where 


c 

n 


J (1-x2) 


J_  { X x) 
o n 


d(x2) 


X = 
n 


th  *1  - 

n value  of 


x for  which 


J (x)  = 
o 


Generally,  ©2  is  dominated  by  the  n = 1 term;  this  will 
be  negligibly  small  when 


lZo 


« 


1 


(E] 


This  will  be  true  when 


-X  Z 

(e  1 ° < 10-2) 


(E] 


or,  in  other  words,  the  transition  region  length  required 


to.  eliminate  any  z -dependence  in  the  potential  is  roughly 
one  drift  tube  diameter. 

The  cold  fluid  Lorentz  force  equation  yields  the 
following  equations  of  motion: 


_d_ 

dt 


(y  vr) 


o 


B 

z 


+ Jjl  / 

m \dr 
o 


+ v 

z 


(E12) 


_d_ 

dt 


(y  r2  0) 


r (v  B 
m \ z r 


(E13) 


_d_ 

dt 


(Y  v ) 
z 


(- 


o® 

3z 


+ v 

r 


(E14) 


where 


B = externally  generated  solenoidal  guide  field 
z 

Bg  = beam  generated  magnetic  field 

v = r 9 
9 


S 


V2  + V2  + Vs  ^ /c 
. r 0 z / 
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Assuming  axisymmetric  0/^0  = 0)  , time  independent 
0/9t  = 0)  motion,  the  relevant  Maxwell  equation  is 


V x B = 


•i 

uo  3 


-» 

n u v 
o o 


Consider  now  the  beam  equilibrium  beyond  the 
tion  region;  it  is  required  that 


v = 0 
r 


and 


_d_ 

dt 


(YVr) 


0 


We  assume 


B = B = constant 
z o 


moc2(Vo  - 1)  = injection  energy  of  electrons 
at  the  anode 


(Ell 


transi 


(Elf 


(El' 


(Eli 


(El! 


■I«WT’>TT 


From  Eq.  (E15) , 


Be  - - m no  % ; J Vr'>  r' dr' 

o 


and  if 


v (r ')  =-  v = constant 
z z 


then 


uu2  a2 


B„ 


r 

2a2 


Then  Eq.  (E12)  becomes 


Y Vq  , , uu2  a2 

—l  _ n v + M ( -E£ — \ v (1_92)_z_  = 
r o 0 m \ c2  / e U V 2a2 


where 


e B 


n = 


o m 


(E20) 
(E21 ) 


(E22) 


0 

(E23) 


while  energy  conservation  gives 


id  = u)  /y  = relativistic  beam  plasma  frequency 
P po  o 

and  in  terms  of  conventional  beam  terminology 


4c2  y 17  0 Y 

o 


where 


I = total  beam  current  in  kA) 

5 

The  two  solutions  to  Eq.  (E23)  are 


If 


= r 0 = 


n r 
o 


ur 


J2° 


2 Y 


Y(1  - P2) 


n2 

o 


1 ± 


1-2 


« 1 


(ju2 

po 

* y 


n- 


(i  - p2) 

Z 


(E2 


(E2 


then  Eq.  (E26)  bas  solutions  given  by: 


from  Eq.  (E26).  In  general,  it  is  assumed  that  one  operates 


with  the  slow  rotational  equilibrium 


6 


(E32) 


P° 

2 Qo  Y1 


The  slow  rotation  of  the  beam  generates  a diamagnetic  field, 

AB  , which  tends  to  cancel  the  main  quide  field  B . The 
z o 

maximum  cancellation  can  be  shown  to  be 


-AB^ 

B 


-E 2_ 


n2  ys 


'ju 2 as 

(^) 


« 1 


(E33) 


so  that  diamagnetic  effects  can  be  ignored. 

If  the  beam  is  immersed  in  a uniform  magnetic  field, 
then  Eq.  (E13)  can  be  re-written  as 


d / _ • lei  r2  \ 

— ( y ra  0 - -1 — 1 — B ) = 0 
dt  \ m 2 z / 


(E34) 


If 


r = initial  radius  at  the  cathode 
o 


and  assuming 


0=0  at  the  cathode 


APPENDIX  F 
BEAM  COMPRESSION 

The  magnetic  compression  of  the  initial  rigid  rotor 
equilibrium  of  the  beam  to  a smaller  radius  is  to  be 
carried  out  in  such  a way  that  miminal  heating  of  the  beam 
occurs.  The  following  analysis  assumes: 

(1)  The  compression  is  carried  out  over  a distance- 
Az — which  is  large  so  that  E^  force  effects  can  be  neglected. 
The  space  charge  electric  field  is  all  radial. 

(2)  Diamagnetic  effects  are  ignored. 

(3)  The  centrifugal  term  is  negligible. 

(4)  The  beam  orbits  follow  equipotential  surfaces. 
This  is  achieved  by  shaping  of  the  waveguide  wall.  Then 
the  v factor  remains  constant. 

(5)  Time  derivatives  are  replaced  with  axial 

derivatives,  that  is, 

d -*  ■+  d d 

——  = v • V =-  v — — =-  c — — 
at  z dz  dz 

(6)  The  axial  magnetic  field,  B (r=0,z)  = B (z) , 

z o 

is  assumed  given. 


T— 


To  lowest  order,  the  equilibrium  orbit  is  defined  by  the 


relation 


fir  2 ( 1 
o \ 


!fi2  ) 


U)  . 

Vi 


Y fi.  7 

i 


or,  given  assumption  (7)  , 


fir  2 = fi . r . 2 

O 11 


Now  set 


r = r +6r  = r + er 
o o o 


where 


e = 


6r 


« 1 


Then 


v = 


dr  ^ 

o d 

c — — + c — e r 
dz  dz  o 


dr  5v 

o r 

c — — + c 

dz  c 


and  to  first  order 


6v  . fi.r.2  ,2.  d2r 

d / r\  _ _1_  / i i ■,  d2 fi  o 

dz  V c J 16  \ c2  1 ro  dz2  dz2 


fi.r.' 
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d2  fi 


Y c / \ r 16  o dz“  / 
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(F6) 


(F7) 


(F8) 


(F9) 


(F10) 
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Ass \aming  that 


'o  daQ 


« 1 


16  a dz2 
and  using  the  following  normalized  units 


f = Q/fi.. 

i 


x = z/Az 


then 
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dx 
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r \ 
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-i  r. 
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d3f"l/S 
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dx3 
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) f e 

n-- r ^ a i ^3, 


A"  d f 

J f1/2  dx2 


(Fll) 


Assuming 


I <x*°>  - 0 


then  Eq . (f!1)  can  be  integrated  once  to  give 


6v  r . 

r _ _i  _d_  / e s 

c Az  dx  \f"' 2 ) 


fi.r. 


= _i  ;_df“1/2  . 1 , *ViN  df1/2 


Az  ! dx 


8 V c 1 dx 
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r /df ' \ , / 

+ — j I I dx 

2 \ c / J ^dx  1 

o 


' i 


n.  az  2 v 

— j f3/ 2 e dx'  ^ 

c J 


(F 12) 


This  equation  can  then  be  integrated  to  give  an  integral 
equation  solution  for  s. 

Alternatively,  one  can  note  that 


r . 

— « 1 
Az 


so  that  keeping  Sv^/c  small  requires  only  that 


_d_  ( £ N 
dx  \W2  J 


= g (x)  < 1 


(F13) 


Then  setting 


q (x)  = C f n ~~ , where  C is  a constant,  (F14) 

’ o dx  o 


gives  a unique  relation  between  s and  f.  inserting 
this  relation  into  Eq.  (F12) , a differential  equation  for 
f is  obtained  which  can  be  integrated  numerically  to 
give  the  field  profile.  Finally,  an  alternate  approach 
is  to  define 
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so  that 


n.Az  2 

d2h  / 1 ,a, 

77  + I fSh 

dx  \ c / 


d2f~1/2 

dx2 


1 ( niri\  rd2f1/2  ,/df1^8  , _ % 

8 \ c ) dxs  + 4 Vdx  ) ! (F16) 


which,  given  a specified  f (x) , can  be  used  to  determine 
e (x)  and  Sv^/c . 


APPENDIX  G 


ELECTRON  BEAM  DIAGNOSTICS 


The  types  of  diagnostics  which  will  be  used  for 
electron  beam  characterization  are  quite  standard.  Be- 
cause of  the  current  densities  and  pulse  length  involved, 
thin  film  dosimeters  are  excluded.  Specific  pieces  of 
apparatus  will  include 

A.  Total  Stopping  Graphite  Calorimeters 
The  incident  beam  energy. 


E = I VI  dt 


is  converted  to  thermal  energy  in  the  graphite . 


Its  temperature  rise 


-L 

J Cp(t) 


dT  = E 


is  monitored  with  a thermocouDle  and  read  out  on  a slow 


scope  trace. 


B.  Segmented  Calorimeters 

A segmented  graphite  array  is  used  to  collect  and 
measure  instantaneous  average  current  density  over  selected 
segments  of  the  beam. 

C . Current  Shunts 

Return  current  shunts  will  be  used  to  measure 
transmitted  current  efficiency.  These  consist  of  parallel 
bands  of  carbon  resistors  selected  to  give  a total  shunt 
resistance  of  about  10  milliohms.  By  tapping  off  the 
shunt  signal  at  various  azimuthal  positions  it  is  possible 
to  infer  the  centering  of  the  beam  in  the  wave  guide. 

D.  Diamagnetic  Loops 

These  single  turn  loops  encircling  the  beam  generate 
an  output  signal  which  is  a measure  of  the  average  trans- 
verse energy  (Ref.  4) . 

E.  X-ray  Diagnostics 

Using  a total  stopping  bremsstrahlung  target,  it 
is  possible  to  infer  anode  information  using  pinhole 

4.  Hammer,  D.A.,  "Pressure  Balance  in  Pulsed 
Plasmas  with  Large  Sel f -Fields , " Phys . Fluids  11 , #6, 
pp.  1260  (June  1974) . 


photography,  and  near-  and  far -field  radiation  profiles 
measured  with  dosimeters. 

F . Magnetic  Loops 

Both  B and  B loops  can  often  be  used  to  infer 
structure  in  the  electron  beam  current.  Care  must  be 

I 

exercised  in  shielding  these  probes  from  electrostatic 

j 

fields . 

G.  Segmented  Faraday  Cups 

These  devices  can  be  used  to  measure  the  radial 
current  density  profile,  and  are  described  in  Ref.  5. 


5.  Schuch,  R.L.,  et  al . , "A  Compact  Faraday  Cup 
Array  for  Measurement  of  Current  Distribution  from  Pulsed 
Electron  Beams,"  Rev.  Sci.  Instr.,  43_,  #8,  pp . 1097  (1972) 


APPENDIX  H 


t 

L 

i 

E 


MAGNET  DESIGN  CONSIDERATIONS 


The  magnetic  field  necessary  for  confining  the 
electron  beam  and  controlling  the  cyclotron  wave  phase 
velocity  will  be  generated  by  a solenoid -type  struc- 
ture. It  will  be  driven  by  a low  voltage  capacitor  bank, 
so  that  the  whole  system  can  be  treated  as  a series  RLC 
circuit.  Assuming  the  circuit  is  under -damped,  the 
winding  current  is 


I 


exp  (-  ou  t/Q)  Sin  <u  t / 1 - 1/QS 
o o V 

(HI) 


where 

L = solenoid  inductance 
C = capacitance  of  driver  bank 
Vq  = initial  capacitor  voltage 
R = total  series  resistance 


I 


■’ 
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0)  = l A LC 
o 


Q = 

R 


Peak  current  occurs  at 


t — tan"1  vQ2  - 1 
° “o 


and  has  amplitude 


max  t tan-1  V Q2  - 1 \ 

V\  ' Sm  > 


(d) 


The  conversion  efficiency  from  initially  stored  capacitor 
energy  to  magnetic  field  energy  is 


1/2  L I2 
_ max 

1/2  C V2 
o 


tan-1  VQ2  - 1 \ 

7q2  - i ^ 


The  ring-up  time  to  peak  current,  t , should  be  chosen  to 

o 

be  several  times  the  magnetic  diffusion  time,  t,  through 
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the  walls  of  the  vacuum  waveguide; 


r = — (At) 
D 


where 

d = electrical  resistivity  of  waveguide  material 
At  = wall  thickness  of  waveguide 

The  intrinsic  Q of  the  solenoid  can  be  estimated 
from  the  expression 


QI  - 2 "o  (R/i) 


where  the  solenoid  is  assumed  to  be  made  up  of  (n/ i) 
windings  per  unit  length;  each  winding  has  width  w and 
extends  from  r = a to  r = b.  The  gap  between  windings 
is  assumed  to  be  negligible.  Then 


/R\  TTC  ' (b  + a) 

( - ^ 

\ l ) w(b  - a) 

\ i ) 

(t)’3”  I1  + <3  + 1/3M} 


. .a ^ 


6 « 1 


(H14) 


then 


f -A 


Q =■  uj  — 7 a (b  - a) 

I op 


If  we  set  t = 5t,  then  from  Eq.  (H2)  we  have 


5 ~ (At)3 
o 


tan-1  V - 1 

V Q*  - 1 


and  using  Eq.  (H16)  one  obtains 


tan-1  V Q - 1 


I Q2Z  -1 


5 5- 


p a (b  - a) 


Having  chosen  a value  for  Q (for  example  from  Eq. 
(H4)),  then  Eq.  (H18)  can  be  used  to  size  the  solenoid. 
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APPENDIX  I 

NUMERICAL  INVESTIGATION  OF  MAGNETIC 
GUIDE  FIELD  PROFILES 

In  the  Auto-Resonant  Accelerator,  the  phase 
velocity  of  the  cyclotron  wave  is  determined  by  the  value 
of  the  local  guide  magnetic  field.  For  ions  trapped  in 
the  wave  it  is  the  spatial  change  in  the  guide  field  which 
thus  governs  the  acceleration  process.  In  this  appendix 
we  examine  by  means  of  a computer  code  the  effects  of 
various  guide  field  profiles  on  the  ion  acceleration 
process . 

Ion  Equations  of  Motion 

In  Ref.  1 the  one  dimensional,  relativistic  ion 
equations  of  motion  for  the  test  ion  (constant  wave 
potential)  case  have  been  derived  and  are  given  by 


1 du 
c dt 


Y-3  Q sin  >!) 


where 


¥ = >P. 


+ - f n(Z' )dz' 

C J-. 


OU  , 

+ — (z  - ct) 
c 


(II) 


(12) 


is  the  phase  of  the  wave  relative  to  the  ion.  In  these 
equations , 

1 

Y = —r  ' 

/ 1 - u3/ca 

^(z)  B - m c is  the  electron  gyrof requency 

e e 

and  is  the  angular  frequency  of  the  cyclotron  wave; 

is  the  (constant)  wave  potential,  e is  the  unit  of 

electronic  charge,  and  m and  m are  the  ion  and  electron 

e 

rest  masses  respectively. 

From  Sq.  (12)  it  follows  that  the  wave  phase  velocity 
is  given  by 


v 

¥ 


uu 

ou  + fi(z)  C 


(13) 


Setting  u = v so  that  the  ion  sits  at  a constant  phase 
point  of  the  wave  (and  thus  experiences  a constant 
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guide  field  given  in  Ref.  l.  This  field  is  again  given 


with 


For  v « l,  Eqs.  (14)  through  (17)  reduce  to 


where  v 


\ 

I 

i 


The  guide  field  given  by  Eqs.  (14)  through  (17),  for 

appropriate  values  of  cp  , i/i  , v , and  uu,  is  the  desired 

o o o 

profile  and  the  one  which  the  magnet  design  should  strive 
to  achieve.  There  are,  however,  several  reasons  why  the 
field  profile  actually  obtainable  may  be  expected  to  differ 
from  this  ideal  field.  The  first  is  that  since  the  guide 
field  magnets  are  made  from  discrete  current  loops,  we 
expect  the  field  to  exhibit  spatial  ripple.  This  effect 
has  been  investigated  both  theoretically  and  numerically 
for  nonrelativistic  ions  in  Ref.  1 and  was  found  to 
present  no  serious  problem  providing  the  ripple  is  kept 
below  5%.  A second  problem  is  just  that  of  the  actual 
coils  producing  the  "wrong"  field  profile.  While  this 
problem  could  be  considered  just  a ripple  problem  by 
fourier  analyzing  the  difference  guide  field,  the  previous 
analysis  has  considered  only  ripple  with  a single  fourier 
component;  also,  in  order  to  interact  with  the  magnet 
design  it  is  useful  to  be  able  to  integrate  Eq.  (II)  for  an 


i 


89 


arbitrarily  specified  guide  field  profile  without  reference 
to  the  "ideal"  field. 


Finally,  the  actual  guide  field  will  differ  from 
the  constant  acceleration  field  because  the  latter  is 
nonadiabatic  near  the  beginning  of  the  accelerator.  The 
theoretical  analysis  on  Auto-Resonant  acceleration  has 
assumed  that 


c dft  . 

—y  — <<  1 

n2  dz 


i.e.,  the  fractional  change  in  the  cyclotron  frequency  in 
a cyclotron  wavelength  is  small.  From  Eq.  (18)  we  see  that 


c dfl  ^ /•  c - 3 eCpo  . 

n2  dz  z=0  w ) me3  sin 


(110) 


Since  for  the  usually  contemplated  injection  scheme 

« c we  see  that  Eq.  (19)  is  violated  at  the  entrance. 
We  thus  expect  to  patch  the  ideal  guide  field  to  that  in 
the  linear  section  through  a short  adiabatic  section. 


du 
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L J ou  me- 


It 


* * *o 

y n „ , 

+ j - ay  + y - r 

O 

(112) 

II 

U 

(113) 

where 

u is 

the  ion  velocity  normalized  to  c and 

y s 

(JU 

z — 
c 

T = 

cut 

(114) 

Q(y) 

C bkg(z| 

U) 

cu  (1 . 707)  Y 

e 

where 

bkg  u) 

is  the  guide  field  in  kilogauss  to  be 

given 

as  a 

function 

of  z in  centimeters . 

Given  a guide  field  profile  (i.e.,  B„_  as  a func- 

KG 

tion  of  z)  we  numerically  integrate  Eqs.  (Ill)  through  (113) 

in  time  for  given  values  of  the  parameters  ecp  /me3 , i|r  , 

o o 

uu,  Yg  and  the  initial  ion  position  and  velocity. 

For  the  cases  run  we  kept  all  parameters  fixed  at 
the  same  value  and  changed  only  the  guide  field  profile 


runs  are 


I 

|! 


II 


ecp 

me 


5.33  X 10 '4 


*o  = • °‘75 

u)  = 1.48  x 109 


(115) 


The  initial  ion  velocity  was  taken  to  be  that  of  the 
cyclotron  wave  at  t = 0,  z = 0.  The  guide  field  at 
z = 0 was  taken  to  be  25  kilogauss  in  every  case.  For 
this  initial  guide  field  the  choice  of  parameters  used 
corresponds  to  a wave  potential  of  0.5  MeV,  and  an  ion 
injection  energy  of  0.25  MeV. 

A listing  of  the  basic  code  is  given  in  table  1-1. 


Discussion  of  Cases  Run 
Test  Case 

To  check  the  code,  a run  was  made  using  the  "ideal" 
guide  field  given  by  Sqs.  (14)  through  (17).  The  numerical 
results  of  this  run  are  shown  in  table  1-2  and  are  in 
agreement  with  analytic  theory.  For  future  reference  we 
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TABLE  1-1  (Continued) 
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TABLE  1-2  CONSTANT  PHASE 


TABLE  1*2  (Continued) 
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give  a plot  of  the  constant  phase  guide  field  as  a 
function  of  z in  Figure  II.. 

Ripple  Effects 

Several  runs  were  made  to  examine  the  effects  of 
spatial  ripple  on  the  guide  field.  For  these  runs  the 
guide  field  was  taken  to  be 

B = B (1  + f sin  z) 
o u 

where  Bq  is  the  ideal  guide  field  given  in  Figure  II  (or 
by  Eqs.  (14)  through  (17)),  f is  the  fractional  ripple  and 
?o  is  the  wave  number  of  the  spatial  ripple.  These  effects 
have  been  analyzed  previously  in  Ref.  1 for  the  nonrela- 
tivist ic  equations  of  motion  and  were  done  again  here  just 
to  check  on  the  relativistic  correction. 

The  results  of  these  runs  are  shown  in  Table  1-3. 
Although  differing  qualitatively  from  the  nonrelativistic 
runs  done  previously,  the  qualitative  features  of  stability 
are  the  same;  that  is,  for  ripple  below  5%  the  ions  remain 


nee  an 


m mug'.1 


s-rv-i*’rt  »•> 


trapped  in  the  wave,  although  the  time  required  to  reach 


30  MeV  is  somewhat  greater. 


Patching  Effects 


Since  the  constant  acceleration  guide  field  is 
nonadiabatic  near  the  accelerator  entrance,  we  considered 

a guide  field  which  starts  at  the  initial  value  of  25  kilo- 


gauss,  decreases  in  an  adiabatic  manner  and  then  patches 
onto  an  adiabatic  portion  of  the  constant  acceleration 


field  in  a smooth  manner.  Specifically,  we  considered  a 
g uide  field  of  the  form 


<1’V[1+T“7r^rp'J 


B u 
o o 


z ^ 2 


B U 
o o 


(1-u  ) Fl+ccz  +a  (z-z)  1 
o i_  1 J 


Z s z 


where  B =25  kilogauss,  u is  the  initial  ion  (and  wave) 
o o 

velocity  and  the  patching  point  is  given  by 


1 + az^  - (l  + az^)a/3 


c 

dz 

is  constant  with 

(1.707)  Ye  a 

£ li  (1  + az  ) 3/3 

o 1 

and  the  field  and  its  first  spatial  derivation  are  contin 
uous  at  the  patching  point  z.  For  z > z,  this  field  is 
just  the  (nonrelativistic)  constant  acceleration  field. 
The  parameters  a and  z^  are  varied  from  run  to  run.  The 
parameter  z^  is  just  the  length  in  centimeters  of  the 
truncated  portion  of  the  constant  acceleration  field.  A 
qualitative  plot  of  this  model  guide  field  is  given  in 
Figure  12  . 

The  results  of  the  patching  runs  are  given  in 
tables  1-4  and  1-5.  In  these  runs,  the  value  of  a = 3.35 

was  used  in  order  to  match  the  model  field  at  365  cm  from 
the  patching  point  to  the  constant  phase  field  given  in 
Figure  II. 


Figure  12.  Qualitative  Plot  of  Patched  Guide  Field 
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Conclusions 


We  now  have  a code  which  will  integrate  the  one 
dimensional,  relativistic  ion  equations  of  motion  for  the 
test  ion  case.  This  code  has  been  used  to  evaluate  the 
effects  of  spatial  ripple  (single  modes)  on  the  ion  accel- 
eration process  and  shows  (as  in  the  nonrelativistic  case) 
that  for  ripple  below  5%,  we  obtain  stable  acceleration  up 
to  the  desired  final  ion  energy  of  30  MeV. 

The  code  has  also  been  used  to  evaluate  the  effects 
of  inserting  an  adiabatic  patching  section  in  the  guide 
field  and  shows  that  for  the  models  tested  we  may  obtain 
stable  acceleration  by  patching  smoothly  to  a guide  field 
which  lies  close  to  the  constant  acceleration  field.  In 
all  of  these  cases  there  is  a small  increase  in  acceleration 
time  and  length. 

We  expect  to  make  further  use  of  this  code  in 
order  to  more  closely  evaluate  the  effects  of  the  actual 
magnetic  field  profile  supplied  by  the  magnet  vendors  and 
designs . 


APPENDIX  J 

ESTIMATE  OF  MAGNETIC  FIELD  RIPPLE 

Consider  a cylindrically  symmetric,  infinitely  long 
solenoidal  magnet  with  hoop  windings . The  windings 
extend  radially  from  r = a to  r - b;  the  windings 
have  thickness  w and  are  spaced  a distance  6 apart. 
There  are  three  regions,  defined  as  follows: 

Region  I r « a 

J = 0 

Region  II  a s;  r ^ b 


J 9 
o 


! “ 

, 1 

00 

/ sin  k 

n 

2 

1 1 r 17  \ 

| W + 6 

T* 

7T 

Lj 

n=l 

l » 

) 

" i 

where 


k = 2rrn/  (w  + 6) 
n 


Region  III  r ^ b 
J = 0 
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The  magnetic  field  produced  is  derivable  from  a vector 
potential  through  the  equation 


B = V x A 


with  only  the  0 -component  of  A being  nonzero.  The 
equation  relating  A (r,  z)  to  J is 


(v‘k\  ‘ Vr,z) 


1 J_  / 5A9  \ 

r Sr  (r  3r  ) 


- u J.  (r,z) 
o 0 


In  each  of  the  three  regions  the  solution  takes  the  form 


A (r,z)  = f(r)  + S a g (k  r)  my k z 
o , n n n n 

n=l 


with  specific  solutions  in  regions  I and  III  being 


A (r , z)  = A r + E A I,  (k  r)  cos  k z 
I o , n 1 n n 

n=l 
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D 00 

O . « 


AIII(r'z)  - T + S cn  VV>  « V 

n=l 


In  region  II  the  solution  is 


A (r  , z)  = Br+—  j — S£—  r2 

II  o r 3moow+5 


+ S l B I (k  r)  + c K (k  r) 
n-1  n 1 n n 1 n 


sin  kn  ~ ) 

Mo  Jo  n k * L1  (knr)  I **  \Z 


where  I1  and  K are  modified  Bessel  functions,  and 
the  modified  Struve  function,  all  of  order  one. 

Boundary  conditions  on  A arise  from  the  con- 

0 

tinuity  of  B and  B at  r = a and  r = b.  Thus, 

» Z 


9Att 

dA 

3A 

II 

and  11 

III 

dz 

ana 

dz 

dz 

r=a 

r=a 

r=b 

r=b 

(J9) 


r 


* ^5 


\ T<  «'  V 


and 


r to  <r  AII> 


r=a 


r=b 


IS, 

; s'1  An» 


= i f (r  A ) 
r dr  HE 


r=b 


r=b 


(J10) 


For  the  terms  with  no  2 -dependence  there  obtains 


2 A = 2 B - n J —7—  a 
o o o o w + 6 


(Jll) 


0 = 2 B - u J — 7— 

o o o w + 6 


(J12) 


or 


. 

& I 


u J w (b  - a) 

A = — _2 

o 2 w + 6 


Noting  that 

(!)- 


J w (b  - a ) 
o 

w + 6 


(J13) 


number  of  ampere  turns  (J14) 
per  unit  length 


then 


" T (7) 


(J15) 
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with  an  ideal  solenoid,  while  the  remaining  terms  corres- 
pond to  the  rippling  produced  by  the  finite  gap  between 


windings.  Also,  I (x)  increases  rapidly  as  x increases 

o 


(I  (x  = o)  =1).  Thus,  the  A terms  should  be  kept  small 
o n 


to  minimize  ripple  throughout  the  beam.  These  terms  are 
defined  by  the  boundary  conditions,  and  one  can  show  that 


k A = u ( ^ ) 
n n o \ i J 


(w  + 6) 


w 

sin  k — 

n 2 


2tt  w (b  - a) 


X (F(x  = k b)  - F (x  = k a)') 
\ n n / 


(J17) 


where 


L (x)  K (x)  - L.  (x)  K (x) 
O 1 1 o 


I (X)  K (X)  - I (x)  K (x) 


F (X) 


(J18) 


APPENDIX  K 


COUPLING  TO  THE  RELATIVISTIC  CYCLOTRON  EIGENMODE 
THROUGH  AN  EXTERNAL  AZIMUTHAL  CURRENT  LOOP 


Previous  analysis  has  indicated  the  requirement  for 
phase  locking  the  relativistic  cyclotron  mode  in  order  to 
prevent  phase  jitter  deleterious  to  the  Auto -Resonant 
acceleration  process.  More  conventional  methods  of  utiliz- 
ing RF  cavities  to  HF  modulate  the  electron  beam  in  the 
cyclotron  mode  appear  not  to  be  feasible  due  to  the  large 
mismatch  in  parallel  wavelength  of  the  free  space  RF  signal 
and  the  cyclotron  mode. 

In  this  appendix,  the  use  of  an  external  HF  azimuthal 
current  loop  to  generate  the  cyclotron  wave  is  examined  and 
the  electrical  characteristics  of  the  loop  derived.  The  use 
of  other  possible  coupling  elements  which  might  be  convenient 
for  specialized  wave-growth  regions--such  as  a helical 
coupler  for  a helical  wave-growth  section — was  not  examined. 


Derivation  of  the  Beam  Response 


We  consider  a relativistic  cyclotron  beam  propagating 
along  a guide  magnetic  field  interior  to  a cylindrical 


conducting  guide  of  radius  a.  Assume  interior  to  the  guide 

and  beam,  a current  loop  of  current  I and  radius  r of 

o 

the  form 


J = 8 I6(r-r  ) 6(z)  exp  (-iuut) 
o 


The  inclusion  of  this  source  term  modifies  the  usual  govern- 


ing equations  for  the  electron  beam-magnetic  field-waveguide 
configuration  (Ref.  1)  by  the  addition  of  an  inhomogeneous 
term: 


± u>  (k-uu)  ± 

□llT  = — r-  ip~  ± 2rri  I 6 (r  - r ) 

k — u>  ± W O 


Here  »p.  n are  the  relativistic  electron  plasma  and  cyclotron 
frequencies,  respectively;  k,  the  longitudinal  wave  vector. 


and  tb~ , the  appropriate  cyclotron  wave  eigenmodes . 


Appropriate  solutions  may  be  expressed  as 


- - Fj1  P r-Y*'  a1' 
2 L 1 J 1 W 

r 


fa*'  ar'-Y*  Jr'jf  a*' 


where  s J^fr);  Y*  ^ Y^fr)  and  <T  ■ * 2rri  I , 


with  the  constants  o.  and  0 determined  from  the  usual 
requirements  of  a conducting  boundary  at  r = a: 


E - Q =>  0 - 0 


= 0 


r=a 


Ez  = 0 =»  i ar[r(0+  + !/)  )] 


= 0 


(K4 ) 


r=a 


Use  of  appropriate  recursion  and  Wronskian  relations  for  the 
Bessel  functions  allow  a,  0 to  be  expressed  in  the  follow- 
ing form: 


Here  D is  the  usual  determinant  of  the  cyclotron  eigermodes: 


D = 


+ + 
? Jo 


? J 


(K7) 


r=a 


and 


P = i 2tt  I r J,  (?”  r ) 
o 1 o 


The  quantity  Ag  is  found  from  the  usual  Fourier  inversion 
of  0 , 0 : 

+ — 

Ag  = • Ml  + M (K8) 

In  Eq.  (K3) # the  first  two  terms  represent  contributions  to 

the  local  loop  response  and  hence  do  not  affect  the  structure 

of  the  cyclotron  eigenmode  propagated  down  the  wave  guide 

away  from  the  loop.  Furthermore,  for  r = r , the  first  two 

o 

terms  vanish  identically. 

For  our  purposes  we  need  evaluate  only  the  contribution 
from  the  last  term  of  Eq.  (K3)  in  order  to  determine: 

(1)  The  eigenmode  response  propagated  away  from 
the  loop,  and  also 
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(2)  The  value  of  A.  at  the  loop  (which  is  required  to 
determine  the  loop's  electrical  characteristics)  . 
In  the  following  analysis,  only  the  contributions  from  this 
last  term  are  considered — the  reader  must  bear  in  mind 
that  the  result  is  rigorously  valid  either  only  well  down- 
stream from  the  loop  or  at  radius  r~rQ. 

Examining  the  structure  of  a,  0 , we  see  that  the 
eigenvalue  poles  arise  solely  from  zeroes  of  D.  These 
modes  are  of  three  general  types : 


(1) 

? a ~ 

0 

(Lower 

cyclotron  mode) 

(2) 

? a ~ 

0 

(Upper 

cyclotron  mode) 

(3) 

+ 

% a ~ 

§~a  s §a 

(Light 

modes ) 

Making  use  of  these  approximate  modes  (which  clearly  separate 
2 2 2 

only  for  uu  < k c ) , we  obtain  the  following  prescription 
P 

for  performing  the  spatial  Fourier  inversion: 


J*(a) 


JL(a) 


J*(r) 

J>) 


JL(r) 

J~(a) 


+ iri 


Jx(5r)  PY1(5a) 

JjU'h) 


light  modes  of  J^^a) 


Performing  these  inversions  is  straightforward,  albeit 


tedious.  More  generally,  if  we  consider  N loops  at 


locations  z and  current  I upstream  of  z,  then  An 

n 0 

is  found  by  simple  superposition:  f 


Ae(r)  = S T rn  0 - ) f O'6  Jl<3'8  r/a) 

n = 0 


x J,  (3  .8  r /a)  exp  (-3.8  Iz-z  I /a) 

1 o n 

uu2  cS  r J.  (2 .9  r /a) 

+ P na  s _o  _1_ o 

„ a 2 s.2  8 c Va  0.357 

(k,  c + uj  ) a 
x P 


J (2 .9  r/a) 

O?? )»irO(*-zn)/c]  (KU) 


where  k^a  = 2.9.  Note  that  only  the  contribution  from  the 
lowest  Bessel  root  of  J, (?a)  has  been  retained.  For  a 
more  accurate .solution,  all  the  roots  should  be  retained  and 


the  resultant  series  summed.  Such  a solution  would  then  yield 


the  exact  calculation  of  vacuum  self  inductance  of  a loop 


interior  to  a conducting  guide,  and  the  response  would  depend 


l 


(K12) 


We  see  from  Eq.  (Kll)  that  the  solution  is  given  by  the 


vacuum  contribution 


(K13) 


Another  important  electrical  characteristic  is  the  value 


compared  to  the  driving  voltage  in  the  current  loop.  This 


may  be  calculated  in  a straightforward  manner  by  solving  for 
A^  and  then  evaluating  |cp|  from  the  relation: 

SA 

i (k  c - uu)  cp  - - (K14) 

2 Sr 

Proceeding  in  this  manner,  we  find: 


co  (r  = z = 0) 


r 

o 

a 


2 

C 


2 » 2 
+ uu  ) 
p 


2 

a 


J.  (2.9r  /a) 
1 o 

0.357 

(K15) 


so  that  this  ratio,  denoted  by  e,  is  given  by: 


e = 


JjlL 

V 


5 5x  lcr3  — 


e c 


Ba  v 


cp 


— - j,  (2.9  r /a)  /0.357 
a 1 o 

- 2 

(— ) J*  (3  .8  r /a) 

\ a / 1 o 


(K16) 


where  I is  the  beam  current  in  kiloamperes;  B,  the  guide 
e 

magnetic  field  in  kilogauss;  a,  the  beam  radius  in  centi- 
meters and  v /c  = uu  /Cl  is  the  ratio  of  phase  velocity  of  the 
wave  to  the  speed  of  light. 

In  Figure  Kl  the  value  of  s and  the  characteristic 
inductance  L of  the  loop  is  plotted  as  a function  of  loop 
radius.  It  will  be  noted  that  as  the  loop  is  taken  close  to 


2 


the  conducting  wall,  the  ratio  of  e improves,  but  at  the 

expense  of  the  loop  inductance.  (This  decrease  in  loop 

inductance  is  due  to  the  conducting  wall  acting  like  a 

shorter  secondary  winding.)  The  value  of  r /a  - 0.7 

o 

maximizes  total  loop  inductance  at  a value  of  L = 13  a nh . 

For  the  proposed  ARA.  experiment  utilizing  a 30  kA 

3 MeV,  1 cm  radius  electron  beam  machine  in  a 24  kG 

magnetic  field  with  — 1.5  x 109 , the  value  of  e for 

r /a  2:  0.7  can  be  calculated  to  be 
o 


s = 2.2 


(K17) 


APPENDIX  L 


SOLUTIONS  TO  SHEATH  HELIX  DISPERSION 
RELATION  FOR  m = 0,  1 

In  a previous  report  (Ref.  1)  the  arbitrary  mode 
number  dispersion  relation  has  been  derived  for  a rela- 
tivistic electron  beam  drifting  through  a sheath  helix 
cylindrical  conductor.  Analytic  solutions  to  the  rela- 
tion describing  the  growth  of  low  phase  velocity 
axisymmetric  cyclotron  waves  have  been  found,  and  pre- 
liminary computational  solutions  have  been  reported. 
This  appendix  presents  the  analytical  solutions  found 
for  both  the  axisymmetric  and  kink  modes  as  well  as  the 
confirming  computational  results.  A stabilizing  influ- 
ence on  the  m = 1 kink  mode  caused  by  increasing  the 
helix  pitch  angle  and  the  magnetic  field  strength  is 
demonstrated.  It  is  shown  that  the  sheath  helix  liner 
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configuration  allows  rapid  growth  of  the  axisymmetric 
cyclotron  wave  which  is  needed  for  Auto-Resonant  accelera- 
tion of  ions.  However,  this  growth  mechanism  is  suitable 
only  if  used  in  conjunction  with  an  additional  active 
mechanism  which  initially  excites  axisymmetric  oscilla- 
tions in  the  electron  beam  before  it  enters  the  growth 
section.  This  is  necessary  because  the  m = 1 "kink" 
mode  grows  at  a slightly  greater  rate  than  the  desired 
m = 0 mode.  Thus,  the  wave  exciter  must  give  the  m = 0 
mode  a sufficiently  larger  initial  amplitude  than  the 
m = 1 mode  so  that  after  the  m = 0 mode  has  e-folded 
several  times  to  the  amplitude  necessary  for  ion  trapping, 
the  m = 1 mode  will  not  yet  have  grown  large  enough  to 
be  disruptive. 

The  analytic  solutions  to  the  sheath  helix  dis- 
persion relation  for  modes  m = 0,  1 are  given  in 
Section  A.  The  relative  growth  rates  and  phase  velocities 
of  these  modes  are  examined.  Section  B gives  details  of 
the  computational  methods  used  to  solve  the  exact  dis- 
persion relation.  A comparison  of  the  numerical  and 
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analytic  solutions  is  presented.  A comparison  of  the 


sheath  helix  and  resistive  liners  is  made  in  Section  C. 
Possible  improvements  of  the  sheath  helix  liner  are  also 


discussed . 


A.  Analytic  Solutions  for  qa  > 1 


The  full  electromagnetic  dispersion  relation  for 
waves  in  a relativistic  electron  beam  propagating  inside 
a sheath  helix  liner  has  been  shown  to  be 


(i/)+  [K1- 3 ujs  - (mk/qa  + q tan  0)2] 


+ 2 qa  K+  [uos  + 'U0  (fi*  + >iT ) + 02  ip+  ip~  ] = 


where 


? a J (fa) 

4 = 


K*  = 


Vl  W 

2 K (qa) 
in 


+ 


tan  8 
a 


e = 


m (k  + m) 
q3a2 


uu2  (k  - uu) 

q2  - k2  - (u2  ; §“2  = - q2  + — — — 

uu  - k + Q 

Preliminary  analysis  of  Eq.  (Ll)  has  revealed 
solutions  for  low  phase  velocity  axisymmetric  waves. 

These  results  are  valid  for  small  pitch  angles.  It  was 
shown  that  as  the  pitch  angle  of  the  helix  was  increased 
the  phase  velocity  and  wavenumber  of  the  most  unstable 
wave  increased,  resulting  in  an  increase  in  the  value  of 
the  parameter  q,  the  argument  of  the  Bessel  functions. 
Thus,  to  extend  this  analysis  to  moderate  and  large  pitch 
angles  the  large  argument  expansions  of  the  Bessel 
functions  are  taken. 

The  small  uu2  case  is  considered  such  that  the 
P 

ju2  corrections  to  % are  retained  only  near  resonance 
with  the  negative  energy  cyclotron  wave,  >ju  = kc  -Q  + e, 
where  e « Q. 
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Using  the  Bessel  function  approximations  a simpli- 


fied equation  for  waves  on  the  helix  is  obtained  with  a 
correction  term  describing  the  effect  of  the  beam  on  the 


wave . 


Thus,  the  dispersion  relation  becomes 


cos 


^ [y  -k2  sin3  8 - SL^L  9 (1+  + 4k  tan  0 )] 


uu2  n 

p 


2 q e 


y 2 - (mk/q  + q tan  8 ) 2 ) + 2q  uu 3 + 2q2  0s  1 = 0 

(L2 ) 


The  free  modes  on  the  helix  are  described  by  the 
zero  beam  density  dispersion  relation 


„3  — 


= k2  sin2  0 + -9  (l  + k2/q2  + 4k  tan  e) 


(L3 ) 


These  modes  can  couple  with  the  beam  cyclotron  modes, 
uu  — k - Q,  only  in  those  regions  where  the  uu,  k for  the 
two  modes  coincide. 
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The  exact  values  of  uu,  k at  resonance  for  the 


axisymmetric  mode  are 


uu  = k sin  0 


1 - sin  0 


An  approximate  solution  for  higher  m modes  may 
be  obtained  for  m/2q2a3  « 1 and  ^ 0/2  « 1,  where 
one  finds 


_ Q sin  0 + m cqs  0 
^m  (1  - sin  0) 


0 + m rtf  0 

1 - sin  0 


It  is  noted  that  the  crossing  conditions  for  the 
axisymmetric  mode,  given  by  Eqs . (L4 ) and  (L5),  are 
identical  to  those  obtained  previously  for  low  phase 
velocity  waves. 


i 
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I 
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The  maximum  growth  rate  for  the  unstable  waves 
may  now  be  obtained  by  examining  the  imaginary  part  of  the 
correction  to  the  cyclotron  resonance  condition  for 
uu  = k - Q.  In  this  case  one  obtains 


r 

m 


uuf.  Cl  cos 2 6 [uu  + q (tan  9+mq/na)2]  \l/2 
p m m 1 


16  q“  tw 


if. 

m 


m 


(L8) 


where  uu  k are  given  by  Eqs . (L6)  and  (L7)  and  q is 
mm  m 

given  by 


+ sin  9 ) + 2m  cos  9 


] 


(L  9) 


For  the  axisymmetric  mode^  Eq.  (L8)  reduces  to 


1/3 

(1  - sin  0)  ) 


(L10) 


Thus,  the  growth  rate  for  the  axisymmetric  mode  is  indepen- 
dent of  magnetic  field  strength;  it  increases  as  the  square 
root  of  the  beam  density,  and  has  a maximum  for  sin  9 = 0.5, 


I 
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i.e.,  8 = 30°.  It  is  also  noted  that  this  mode  becomes 
stable  for  s!a  8 < 0.  It  can  be  shown  that  the  cyclotron 
mode  is  unstable  only  when  the  helicity  of  the  helix  and 
of  the  electron  trajectories  has  the  opposite  sense. 

It  is  now  desired  to  examine  the  m = 1 growth 


rate  for  comparison  with  that  of  m = 0.  Equation  (L3) 
may  be  rewritten  to  obtain 


““m  = tan2  0 + ^ (l  + + 4k  tan  8 'j 


In  order  to  effect  a simple  comparison  with  the  m = 0 

growth  rate  given  by  Eq.  (L10) , an  underestimation  of  the 

m = 1 growth  rate  is  made  by  replacing  ou  in  the 

m 

numerator  of  Eq.  (L3)  by  q tan  8. 


Then  the  m = 1 growth  rate  may  be  put  in  the 


form 


r “■  ( -f-  sin  8 (1  - sin  8)  (1  + n)  ) 


where 
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\ 
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t 


x . -,  « — ■ —*■*•-■* 


•*  -u  P*P 


1 

[4 ( Qa/ c) 2 tan  9 (tan  9 + c/Qa) ] 


(1*13 ) 


A comparison  of  Eqs.  (L10)  and  (L12)  shows  that 
m = 1 grows  faster  than  m = 0 for  all  values  of  0 and 
tan  9 - Indeed  it  is  seen  that 


1 1 1/2 
— = (1  + Tl)  > 1 

: 0 


(1*14) 


since  r]  is  always  greater  than  zero  for  positive  pitch 
angles,  it  is  evident  from  Eq.  (L13),  however,  that  the 
ratio  of  m = 0 to  m = 1 growth  rates  is  reduced  by 
increasing  either  the  gyrofrequency  or  the  pitch  angle. 
Although  the  desired  mode  never  grows  faster  than  the 
kink  mode,  the  two  growth  rates  can  be  brought  fairly 
close  at  reasonable  values  of  pitch  angle  and  magnetic 
field  strength.  A wave  initiator  can  be  used  to 
preferentially  excite  the  m = 0 mode  before  the  wave 
enters  the  growth  sector,  after  which  this  liner  could  be 
used  to  safely  grow  the  axisymmetric  wave.  Of  course. 
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the  m = 1 wave  would  be  growing  at  a slightly  greater 
rate  but  its  initial  amplitude  would  be  that  of  the 
thermal  noise  in  the  beam  while  the  initial  amplitude 
of  the  m = 0 would  be  that  generated  by  an  external 
driver.  Thus  the  m = 0 mode  could  probably  be  grown 
to  sufficient  amplitude  for  ion  trapping  before  the  kink 
mode  could  disrupt  the  beam. 

B . Numerical  Solutions  of  the  Full  Sheath  Helix 
Dispersion  Relation 

The  full  dispersion  relation  given  by  Eq. 
is  not  tractable  to  an  exact  analytic  solution  because 
of  the  transcendental  nature  of  the  Bessel  functions. 

As  a result  approximations  to  these  functions  have  been 
employed  in  obtaining  the  analytical  relations  for  phase 
velocities  and  growth  rates  of  cyclotron  waves  in  the 
system.  In  addition,  the  analytic  work  is  generally 
restricted  to  a specific  range  of  values  of  parameter 
space  in  order  for  these  approximations  to  remain  valid. 
Despite  these  limitations  the  analytic  approach  to  a 
solution  to  the  dispersion  relation  is  most  important 
because  it  gives  results  that  allow  a qualitative 
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understanding  of  the  processes  involved  and  directly  yields 
the  scaling  of  the  frequencies  and  growth  rates  with  the 
physical  parameters  of  the  problem.  This  approach  also 
yields  approximate  quantitative  results  which  are  accurate 
to  the  degree  that  the  approximations  used  are  correct 
representations  of  the  actual  functions. 

Once  the  analytic  results  are  obtained  it  is 
desirable  to  check  them  using  a computer  program  to  solve 
the  exact  dispersion  relation.  This  has  been  done  for 
the  sheath  helix  dispersion  relation  using  the  dispersion 
solver  described  in  Ref.  1.  This  solver  employs  the  full 
dispersion  relation,  Eq.  (LI),  and  calculates  the  required 
Bessel  functions  of  complex  argument  from  their  integral 
definitions.  A modified  Newton -Raphs on  iteration  proce- 
dure is  followed  with  the  complex  frequency  as  the  iteration 
variable.  Basically,  the  program  computes  the  value  of  the 
dispersion  function  €,  the  left  hand  side  of  Eq.  (Ll)  , 
for  a fixed  real  k and  an  initial  guess  for  uu.  Then 
is  computed  and  uu  is  stepped  to  bring  the  value 
of  € to  zero.  The  result  is  considered  to  be  converged 


when  the  incremental  change  in  uu  necessary  to  further 


reduce  the  value  of  € is  below  some  predetermined 
limit . 

An  additional  facility  of  the  computer  dispersion 

solver  is  that,  besides  allowing  confirmation  of  wave 

growth  due  to  the  mode  crossings  studied  analytically, 

it  allows  investigation  -»f  other  areas  of  uu,  k space 

for  possible  unstable  waves.  Two  such  unstable  m = 1 

modes  have  been  found  in  past  liner  studies . Both  have 

subsequently  been  examined  theoretically  in  Ref.  1 and  in 

Appendix  M of  this  report,  and  were  found  to  be 

stabilized  for  uu2/fi8  less  than  about  1/4.  This  was 

P 

then  verified  numerically. 

The  present  computational  study  has  involved  an 
extensive  examination  of  parameter  space  for  the  sheath 
helix-electron  beam  system.  The  m = 0,  1 wave  growth 
has  been  examined  for  parameters  in  the  following  ranges: 


tan  0 

0.01 

-♦  3.0 

Qa/c 

0.5 

-♦  7.0 

0.01 

0.4 

All 

the  qualitative 

features 

of  the 

analytic  solutions 

have 

been  confirmed. 

There 

is  also 

excellent  quantita 

tive 

agreement  between  the 

rju,  k position  of  maximum 

wave  growth  given  by  Eqs . (L6)  and  (L7)  and  the  computer 

results.  In  Figure  LI  the  phase  velocity  of  the  most 
unstable  wave  versus  tan  8 is  presented  for  constant 
Qa/c  = 1.  The  solid  lines  are  the  theoretical  predictions 
for  m = 0,  1,  and  the  data  points  represent  the  computa- 
tional results . The  growth  rates  of  the  waves  found  from 
the  computer  solutions  to  the  dispersion  relation  are  also 
in  good  agreement  with  the  theory,  Eq.  (L8)  . The  m = 0 
and  m = 1 growth  rates  are  shown  in  Figure  L2  for 
fia/c  = 1 and  plotted  versus  tan  0.  Again  the  solid  lines 
are  the  theoretical  predictions  from  Eq . (L8)  and  the 

data  points  are  the  computational  results  . 
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C.  Summary  and  Conclusions 


The  sheath  helix  liner  has  been  investigated  for 
wave  growth  of  the  axisymmetric,  m = 0,  and  the  "kink," 
m = 1,  modes  through  interaction  of  the  negative  energy 
beam  cyclotron  mode  with  positive  energy  modes  associated 
with  the  helical  structure.  The  axisymmetric  mode  has  a 
phase  velocity  given  by 


v /c  = sin  0 

P 


(L15) 


and  a growth  rate  given  by 


uu  «-*• 

= ( 4 c 2 ^ 9 t1  - & 9)  ) 


(L16) 


The  disruptive  "kink"  mode  has  been  shown  to  grow  faster 
than  the  m = 0 mode  for  all  values  of  pitch  angle  and 
magnetic  field  strength,  although  it  has  been  demonstrated 
that  increasing  the  pitch  angle  and/or  the  field  strength 
has  a stabilizing  influence  on  the  m = 1 mode.  in  fact 
for  Qa/c  tan  9 s 2.4  the  m = 0 growth  rate  is  at  least 
80  percent  of  the  m = 1 rate. 
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This  work  shows  that  in  this  parameter  regime 
such  a configuration  could  be  a useful  wave  growth 
structure  provided  a wave  initiator  can  be  used  prior  to 
the  growth  section  to  generate  axisymmetric  perturbations 
in  the  electron  beam.  Such  a "tickler"  would  also  serve 
to  select  a dominant  single  frequency  component  in  the 
wave.  This  frequency  component  would  then  have  a much 
larger  initial  amplitude  upon  entering  the  growth  region 
and  would  thus  be  the  major  component  in  the  large 
amplitude  wave  grown  there. 

A comparison  of  the  sheath  helix  liner  with  the 
previously  studied  resistive  liner  shows  that  the  maximum 
growth  rates  for  the  helix  liner  are  at  least  3.5  times 
greater  than  those  for  the  resistive  liner,  for  wave 
phase  velocities  corresponding  to  250  keV  ion  injection 
energies.  If  higher  injector  energies  are  allowed,  or 
if  the  wave  may  be  slowed  down  after  growth  but  prior  to 
ion  injection,  then  even  faster  growth  rates  and  thus 
shorter  liners  are  obtainable. 

Future  efforts  need  to  be  directed  toward  a study 
of  the  finite  geometrical  effects  associated  with  a 
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physical  helix,  including  finite  spacing  between  the 


C 


wires.  Special  consideration  should  be  given  to  a 
search  for  a structure  that  could  be  prejudiced  against 
the  ubiquitous  m = 1 disturbance  without  a significant 
reduction  in  the  growth  rate  of  the  desired  axisymmetric 
mode . 


APPENDIX  M 


SUPPRESSION  OF  THE  m = 1 HYDRODYNAMIC  INSTABILITY 
FOR  THE  SHEATH  HELIX  GROWTH  SECTION 


Methods  for  suppression  of  the  hydrodynamic  insta- 
bility for  the  edge  ring  resistive  liner  growth  configu- 
ration have  been  reported  previously  (Ref.  1).  In  this 
appendix,  the  requirement  on  magnetic  field  strength 
sufficient  to  quench  this  instability  for  a sheath  helix 
growth  section  is  derived. 

Computational  analysis  of  a helical  travelling  wave 
reactive  growth  structure  indicates  the  presence  of  a 
strong  nonaxisymmetric  (m=l)  instability  occurring  at 
very  long  parallel  wavelength.  The  governing  dispersion 
relation  has  been  derived  previously  and  is  here  re-stated 
for  completeness : 


(1  + |34)+)  (1  + giif) 


2 qa  u*  K4 


K+  (k2  - q2  ) - as  k2 


= 0 


(Mia) 


where 


a = 


m q 

— + , tan  0 
qa  k 


(Mlb) 
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Here  ;u,  k are  the  frequency  and  parallel  wave  number  of 
the  wave,  respectively;  m,  the  azimuthal  mode  number; 

9,  the  pitch  angle  of  the  helix  and  a,  the  radius  of 
the  helical  wave  guide.  Moreover 


i/T  = r J (Ta)/J  (fa) 
m m±l 


(Mif) 


where 


= uuc 


-k2  - -2- 


uu2  (k  - yu) 


1/2 


k - rju  ± n 


(Mlg) 


Here  ijj  is  the  relativistic  beam  plasma  frequency  and 
Q,  the  relativistic  cyclotron  gyrofrequency . 
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We  now  consider  the  m = 1 kinking  instability  for 
the  low  k limit. 

± ± 
For  k -»  0 and  £ a small,  the  functions  K*  , ij) 

assume  the  following  limits: 


ill*  -*  4 


i|T  -*  2 (ra)2 


and  with  the  further  assumption  that  > 1,  the  dis- 

persion relation  assumes  the  following  simple  form: 


L= 


uo  (uu  + Cl) 


+ 1 + h"1  (qa)  = 0 


where  h(qa)  is  given  by 


h (qa)  = 


2 K* (qa) 

qa (K+  (qa)  + tan2  6) 


The  function  h(qa)  is  monotonic  decreasing  in  qa,  from 
a maximum  value  of  2.  Utilizing  this  upper  bound  on 
h(qa)  in  Eq.(M4),  we  obtain  as  a sufficient  condition 
for  stability 
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APPENDIX  N 

TAPE  HELIX  DISPERSION  RELATION 

This  work  begins  a consideration  of  the  effects  of 
the  finite  geometry  associated  with  a helical  liner 
structure  made  with  wires  of  finite  size  and  with  finite 
spacing  between  the  wires.  Previous  analysis.  Appendix  L 
of  this  report,  has  dealt  with  the  sheath  helix  model,  a 
physically  smooth  conducting  sheet  at  the  outer  boundary 
of  the  electron  beam  with  perfect  conductivity  in  a 
direction  oblique  to  the  beam.  Such  a liner  could  be 
made  by  winding  many  fine  wires  to  form  a helical  solenoid. 
For  wavelengths  much  larger  than  either  the  wire  thickness 
or  the  spacing  between  the  wire,  this  should  compare 
favorably  to  the  mathematical  model  of  the  sheath  helix. 

The  present  model  is  the  tape  helix  wound  with  a single 
tape  of  width  6 and  turn  to  turn  spacing  d. 

The  dispersion  relation  for  waves  in  the  tape  helix- 
electron  beam  system  is  derived  and  the  normal  modes  of 
this  system  analyzed  and  compared  with  those  of  the  sheath 
helix.  The  dispersion  relation  for  these  normal  modes  in 
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the  absence  of  the  beam  is  examined.  Coupling  of  these 
helix  modes  with  the  beam  cyclotron  mode  is  considered 
under  certain  restrictions,  and  growth  rates  are  obtained. 
It  is  shown  that  for  mode  coupling  in  certain  regions  of 
uu,k  space  the  physics  of  the  interaction  is  only  slightly 
modified  from  the  sheath  helix  results  by  the  finite 
geometry  of  the  tape. 

Derivation  of  the  Dispersion  Relation 

In  Figure  N 1 a schematic  drawing  of  the  tape  helix 

model  is  shown.  The  tape  width  6,  the  turn  to  turn 

spacing  d,  and  the  radius  of  the  winding  a,  specify  the 

physical  system.  The  pitch  angle  9 is  related  to  these 

2nd 

parameters  by  cot  9 = — — - . The  mathematical  model 
for  the  tape  conductivity  is 


. a T d /2TT 

a (r,6  ,z)  = ; — H 


4a6R  \ 


^ (z+6/2)  - (6+2TTn)) 


x h - h(^~-  (z-6/2)  - (6+2TTn)yi  • 6 (r  - a) 

where  R is  the  resistivity  per  unit  length  and  the  H 
functions  are  the  Heaviside  or  step  functions  with  the 
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1 aa  (nrr6/d) 
8aR  (nTT6/d) 


ime 


k = g + 2Kjt/d 


are  now  found 


(N17) 


where  K = 


= - a • — k 

k,m  q 


K (q) 

m 


— K + p*  k 


K*  (q) 

m 


(N18) 


and 


qa  = k8  -uj: 


(N19) 


where  k (z)  and  I (z)  are  modified  Bessel  functions  of 
m m 


order  m and  argument  z . 

Applying  Maxwell's  equations  to  the  boundary  region 

provides  a set  of  matching  conditions  sufficient  to  specify 

the  four  sets  of  undetermined  coefficients,  a ,3 

k,m  k,m 

a'  , and  3 ' 
k,m  k,m 

Three  of  the  "jump"  conditions,  governing  the 
tangential  electric  field  components  and  the  components 
of  magnetic  field  parallel  to  the  tape  direction,  are 
satisfied  for  each  Fourier  component  separately: 


(1) 

k,m 

Es 

(beam) 

k,m 
’ E6 

(liner)  = 0 

(2) 

Bk,m 

z 

(beam) 

- Bk'm 
z 

(liner)  = 0 

(N20) 

(3) 

k,m 

B || 

(beam) 

k,m 
- B || 

(liner)  = 4TTj^*m  (beam) 

The  fourth  matching  condition  relates  the  magnetic  field 
components  perpendicular  to  the  tape  to  the  tape  current 


-r-*-'  - . • r 


and  the  parallel  beam  current.  This  condition  couples  all 

the  Fourier  modes.  It  is 

/ . \ #.  % 1c i m ( . . k / m ..  . 

(4)  Bj_  (beam)  - Bx  (liner)  = 4ttjx  (beam) 


. T k ',m ' . 

+ 4tt  / v' ' EH  (liner) 

Lu  f k k § rn  -m 
k 


Using  these  conditions  it  is  possible  to  solve  for 

the  coefficients  a/  , 0 ' , and  0,  in  terms  of  the 

k,m  k,m  k,m 

coefficient  ol  and,  after  considerable  algebra,  obtain 
k , m 

the  following  governing  equation  for  the  system: 


Vm  W 


7 — ??-(r!;/d)  a c (v<  , 

L>  2aR  (rrr/d)  k',m+r  m+r v n-r 


where 


k = 0 + 2nn/d  . 

n 


Here  £ (k  ) is  the  dispersion  function  for  Fourier  modes 
m n 

of  order  m and  wavenumber  k in  a beam  vacuum  system, 

n 

given  by 


„ ,,  . / +/mk  \ .+  \ 

W ‘ (*  (T*  + ®)  - 8 ’K) 


X ^"(^K'JU)  + ^ 


where  and  K are  as  previously  defined  and 

3 = (id-  - — 0X1 

\ (k-uu)  / 

The  function  € (k  ) is  the  sheath  helix  dispersion 
ran 


function  for  Fourier  modes  of  order  m and  wavenumber  k 


n 


(see  Appendix  L)  and  is  given  by 


^(k^)  = u (>p++  th  ) ^'jj3ot2S  - casQ  + qsin  9 j K2  J 


2qK(^0  +[^}  + sto  9 j ^(<»  - [^}  8]  0’) 


may® 


(N25) 


If 

region  at 

sin  ( rrr6/d) 
(rxr6/d) 


the  helical  conductor  fills  the  entire  boundary 
r = a,  i.e.,  6/d  = 1,  then  the  weight  function 

is  nonzero  only  for  r = 0,  and  Eq.  (N22) 


reduces  to  the  form 


a 


3,  A,m 


TT 

2aR 


0 


(N26) 


for  each  a coefficient. 

In  this  case  each  angular  mode  number  satisfies  a 
separate  dispersion  relation  independently  of  the  other 
angular  modes.  Eq.  (N26)  is  the  sheath  helix  dispersion 
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relation  for  finite  resistance  in  the  helical  conductor. 

For  R = 0 this  becomes  the  perfectly  conducting  sheath 
helix  liner  described  in  Appendix  L of  this  report. 

On  the  other  hand,  when  the  tape  width  is  reduced  so 
that  there  is  space  between  turns  of  the  tape  then 
combinations  of  angular  modes  are  necessary  to  satisfy 
the  boundary  conditions.  If  the  tape  width  6 is  reduced 
such  that  6/d  « 1 then  each  angular  mode  is  weighted 
approximately  equal  in  the  solution  to  Eq.  (N22)  and 
this  equation  becomes 


a C ( 3 + \ = n 

3,  l,m  m \ d / 


a„ 


2TT(j-r) 


2aR  L “0  , i-r  ,m+r^m+r  d 


) 

(N27 ) 


It  can  be  shown  that  the  right-hand  side  of  Eq.  (N27)  is 
independent  of  i or  m for  all  2 and  m such  that 
2+m  = constant.  A recursion  relation  for  the  a's  thus 
obtains 


a 


3, 2,m  m S, 


a / 2rr2\ 

'3+~)  = a 


3 * N ,0  o V 


a + 2^N  . 

3+  d J 


(N28) 
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For  the  case  where  the  wire  resistance  R - 0,  the 


following  dispersion  relation  for  modes  of  index  N is 
obtained,  using  the  recursion  relation  Eq.  (N28) : 

€n  ( P + 2n  (N-m)/d  ) 

£ (P+2tt  (N-m)/d)  = ° 


(N29) 


Normal  Modes  in  the  System 


Before  examining  the  narrow  tape  helix  dispersion 
relation,  Eq.  (N29)  , it  will  be  of  some  benefit  to  con- 
sider the  normal  mode  structure  in  this  system.  A normal 
mode  is  characterized  by  a parameter  N and  contains 
contributions  from  all  angular  modes  such  that  m+  l = N. 
Thus  a normal  mode  may  be  written  in  the  form 


_ V’  iz  ( 0+2tt  (N-m)  /d)  im<3 

E = ) a e e 

N L.  N-m,m 
m 


(N30) 


The  recursion  relation,  Eq.  (N28)  , allows  a 
calculation  of  the  relative  amplitudes  of  these  components 


to  obtain 


- , , 2n  v 

EN  = ^,o  L 6 ( 9«n  (H-m)  /d)  eXp  "L1”3  « )J 

m 

m 


(N31) 
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Normal  Mode  Dispersion  Relation  (uu^3  = 0) 

The  general  dispersion  relation  for  the  narrow  tape 
helix,  Eq.  (N29),  may  be  simplified  in  the  limit  of  zero 
beam  density,  where  the  expression  reduces  to 


o - £ K (q  ) I ' (q  ) l I ' (q  ) K ' (q  ) 

mmmm\  mmmm 


m 


+ ( tt  + 3tan<3 ) i (<l>  Km  ) 
\ q / m m / 


(N32) 


Now  if  none  of  the  q 1 s approach  zero  the  large  argument 

m 

expansion  of  the  Bessel  functions  may  be  used  to  obtain 
an  approximate  dispersion  relation 


0 « (w2  - q^  tan *9)  T (N33) 

m ra 

The  factor  jj2  -q^tan29  is  common  to  each  term  in  the 
sum,  and  thus  setting  this  factor  to  zero  provides  a 
solution  to  Eq.  (N33) « From  the  definition  of  q^  this 
may  be  written  in  a simplified  form 


jo  a 

c 


3a  sia  9 + N 


(N34 ) 


This  is  formally  the  same  dispersion  relation  found  for 
the  sheath  helix  liner  for  a mode  of  order  N and  wavenumber 
0 . However,  the  present  relation  characterizes  the  uu,P 
dependence  of  a normal  mode  of  parameter  N which  is  com- 
prised of  a combination  of  angular  modes  with  different 
wavenumbers  k = 3 + 2tt  (N-m)  /d . 

To  examine  this  normal  mode  structure,  we  note  that 

£ (3+2n  (N-m)  /d)  = uu2/  (K  12  (q)  I (q)  I (q)  ) (N35) 

m m m+l  m-l 

for  uu  2 - 0 ; and  for  q2a2  > 1 , 

P 

£ (k  ) ^ q 2 uu2/4  (N36) 

mm  m 

Thus  the  contribution  of  the  angular  components  to  the 

normal  mode  is  weighted  by  the  relative  value  of  the 

parameter  q = (k  2 - uu2)172  , so  that  the  normal  mode 
m m 

can  be  written 

E T = Z aM  _ (q  V qm2)  exp  i (mo  +k  z)  (N37) 

N N , o o m m 

m 

where  again  k = (3  + 2tt (N-m) /d)  . 

m 

Since  the  parameter  3 is  variable  only  over  the 
range  - rr/d  < 3 < rr/d,  it  is  seen  that  the  m = N 
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components  of  the  helix  modes  occurs  by  means  of  currents 
generated  in  the  tape. 

The  growth  rate  of  a normal  mode  may  be  found  by 

Taylor  expanding  the  general  dispersion  relation,  Eq.  (N29) , 

about  u)  3 = 0 and  keeping  the  uu  3 correction  only  for 
P P 

the  term  whose  uu,k  approximately  satisfy  the  beam 
cyclotron  dispersion  relation 


uu  » k - « P + 2tt  (N-m)  /d  - Q 


(N38) 


The  peak  growth  occurs  at  the  3 a such  that 

o 

3 a = (ft  + N cos  0)/(l-«y6)  - (N-m)  cot  0 (N3  9) 

o 


Since  the  parameter  3a  is  defined  only  in  the  range 

- < 3a  < the  above  condition  implies  a restriction 

2 2 

on  the  angular  components  of  a given  normal  mode  that  can 
couple  to  the  cyclotron  mode. 

The  peak  growth  rate  under  these  assumptions  is 


id  zClx?  e (uu  +q  (tanS+Tr))3 
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where  uu  is  found  from  Eq.  (N38)  for  g = 3 . 

n o 

This  growth  rate  has  been  written  in  form  that  allows 
simple  comparison  with  the  sheath  helix  results  in  Appendix 
L of  this  report-  In  fact,  the  ratio  of  the  growth  rates 
may  be  written 


rw  tape 
N , n 

T sheath 
n 


} 1/q  2 

u m 


(N41) 


Thus  the  tape  helix  growth  rate  is  always  less  than  that  of 


a sheath  helix  with  equal  pitch  angle.  Since 


q.s  s (3  + 2 tt  (N-l)/d)  i - uu: 


it  is  seen  that  the  growth  rate  will  be  a maximum  for  a 
normal  mode  if  the  cyclotron  mode  couples  with  the 
dominant,  angular  component  of  the  mode,  N = m.  In 


the  limit  3 « tt/ d the  fundamental  component  is 

much  larger  than  any  other  component  and  the  infinite 

sum  22  1/q  s is  approximately  equal  to  1/q  2 . Only  in 
mm  N 
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this  limit  does  the  tape  helix  growth  approach  that  of 
the  sheath.  The  lower  growth  rates  for  the  tape  helix  may 
be  thought  of  in  terms  of  an  increased  "inertia"  provided 
by  the  other  angular  components  of  the  normal  mode  which 
do  not  couple  to  the  cyclotron  mode. 

Thus  it  is  possible  to  arrange  the  crossing  of  the 
cyclotron  mode  with  the  N = 0 normal  mode  in  a region  of 
'D,  k space  where  the  normal  mode  is  predominantly 
axisymmetric . At  the  same  time,  the  N = 1,  m = 1 crossing 

can  be  avoided  entirely  and  the  N = 2 , m = 1 normal  mode 

crossing  occurs  in  a region  where  the  m = 2 omponent  is 
dominant,  Eq.  (N41),  thus  indicating  a reduced  growth  rate 
for  m — 1.  This  can  occur  if  fia/c  is  in  the  range 

C~-“  (1  - sin  6)  > “ > r?'  ^ (1-3  sin  9),  0 I largest  (N42) 

2 c l 2 J 

For  such  a case  the  axisymmetric  growth  rate  will  be 
approximately  that  obtained  for  the  sheath  helix  while 
the  nonaxisymmetric  growth  will  be  reduced  from  the  sheath 

helix  growth  by  a f actor  (( 1/q? ) /72  (1/q2)  ) . 

1 m 

m 
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Summary 


The  preceding  analysis  has  shorn  that  the  tape  helix 
couples  less  effectively  to  the  beam  cyclotron  mode  than 
does  the  sheath  helix.  In  addition,  the  wave  grown  in  the 
beam,  since  it  is  a normal  mode  of  the  system,  contains  all 
angular  modes.  However,  it  appears  possible  to  grow  an 
almost  pure  axisymmetric  mode  in  the  beam  by  arranging 
the  crossing  with  the  cyclotron  mode  to  occur  at  small 
values  of  the  parameter  3.  This  also  allows  the  tape 
helix  growth  rate  to  approach  that  of  the  sheath  helix. 

Thus  it  does  appear  possible  to  grow  the  desired  mode 
with  such  a liner  structure. 

Analysis  also  indicates  that  the  m = 1 "kink" 
instability  in  the  beam  will  have  a growth  rate  reduced 
from  that  for  a sheath  helix  liner  under  certain  condi- 
tions . In  this  case  the  relative  growth  rates  of  the 
m = 0 and  m = 1 waves  in  the  beam  have  a more  favorable 
value  than  for  the  sheath  helix. 


APPENDIX  0 


A MODIFIED  RESISTIVE  GROWTH  MECHANISM  USING 
DISCRETE  LOOP  DRIVES 

In  this  appendix,  the  use  of  discrete,  uncoupled 
loop  drives  connected  to  external  RC  lumped  circuits  is 
examined  for  growth  of  the  axisymmetric , low  frequency 
cyclotron  mode.  It  is  shown  that  a considerable  improve- 
ment in  growth  rate  over  the  usual  resistive  type  growth 
mechanisms  can  be  achieved  if  the  value  of  capacitance  is 
chosen  to  be  resonance  with  the  inductance  of  the  loops 
at  the  driving  frequency.  Operation  at  frequencies  well 
removed  from  resonance  reduces  the  growth  rate  (or, 
more  precisely,  increases  the  growth  length)  to  a value 
comparable  to  the  resistive  liner  configurations  dis- 
cussed in  Ref.  1. 

Arguments  eure  presented  for  achievement  of  a system 
which  is  relatively  stable  to  the  m = 1 kinking  type  mode. 


Voltage  handling  capability  is  also  investigated. 


It  may  be  required  to  either  detune  the  rings  for  the 
final  tens  of  centimeters  of  growth  to  avoid  power 
handling  problems  (thereby  resulting  in  an  increased 
growth  length) , or  to  include  more  rings . This  second 
possibility,  if  carried  to  the  limit,  would  result  in  a 
simple  modification  of  the  usual  edge  ring  liner  configu 
ration  in  that,  rather  than  utilizing  shorted  rings,  the 
rings  would  be  connected  to  a distributed  capacitive- 
resistive  network.  The  consequences  of  such  a configura' 
tion  have  not  yet  been  investigated.  Here  we  restrict 
our  attention  to  discrete  loop  drives  . 

Governing  Equations  for  Discrete,  Uncoupled  Loop  Drives 
In  Appendix  K,  the  amplitude  of  the 
cyclotron  eigenmode  established  by  a series  of  loop 
drives  is  derived.  For  loops  located  at  monotonically 

ordered  positions  z , of  current  I and  voltage  V , i 

n n n 

the  relevant  equations  may  be  put  in  the  following 
tractable  form: 
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For  each  loop,  the  loop  voltage  is  given  by 


N 

V = L i + a £ I sin  Q (z 
N N _ n N 

n=0 


z )/c 
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2rra  • 
c A0o 


(01) 


with  the  8 component  of  vector  potential  between  the 
NTh  and  N + l5 T loops  given  by 
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2rra 
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I sin  Q ( z - z ) /c 
n n 


+ A 


0o 


(02) 


The  principle  of  superposition  has  been  used  to  include 

the  effects  of  an  incident  wave  of  amplitude  A . Here 

9o 

L is  the  inductance  of  the  loop;  Q the  relativistic 
cyclotron  frequency?  a,  the  beam  radius;  and  a,  a 
coupling  coefficient  of  the  loop  drive  to  the  wave  given 
by 


a = 8TTsa 


(k2c2 


2 


0a 

c 


geometrical 

factors 


(03) 


Here  is  the  relativistic  electron  plasma  frequency 

and  k,  the  total  wave  vector  of  the  wave. 
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The  current  I is  strictly  evaluated  at  the 
n 

retarded  time:  t - (z  - z )/c.  However,  for  low 

n 

frequency  drives,  the  retardation  factor  can  be  neglected. 
This  requires  that  the  region  z - z^  be  limited  to 


(z  - z ) < 2ttcA) 
n o 


where  m is  the  driving  frequency. 

We  note  two  interesting  features  of  these  equations: 


(1)  All  circuits  upstream  of  the  N'M  loop  con- 
tribute to  V through  the  cyclotron  wave  coupling 

N 

coefficient  a,  while  the  current  in  the  NTH  loop  only 
contributes  through  the  loop  inductance  L. 

(2)  Each  loop  propagates  downstream  a standing 
cyclotron  wave  contribution. 

It  is  thus  obvious  that  at  least  two  loops  are 


required,  properly  phased,  in  order  to  generate  a travelling 
cyclotron  wave  by  an  appropriate  phasing  of  these  two 
standing  waves . 

We  thus  turn  our  attention  to  a specific  configura- 
tion of  two  loops  which  are  connected  to  external  RC 


circuits.  Let  these  loops  be  located  at 


V 


with  an  incident  wave  given  by 


A = exp  ( iuut  - inz) 

u o 


The  steady  state  equations  governing  V , V and 

o 1 

AQ (z  > d)  are  found  from  Eqs . (01)  and  (02) ; 


V = iouL  I + i2na  yu/c 
o o 


V = icuL  I + icua  sin  Q d/c  + l2r7a!U  e“lQ  d/c 
1 c 


2rra  A (z  > d)  — a 1 sin  Qz/c  + a I,  sin  — (z  - d) 
o o 1 c 


+ 2rra  exp  - iQz/c 

Let  Zo  and  Z1  t5ie  imPedance  of  the  external  RC 
circuit  elements  attached  across  loop  N = 0 and  N ■ 
respectively,  so  that: 


V + Z I =0 
o o o 


V1  + Z1  J1  = 0 
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Using  Eqs . (07)  through  (010b),  we  inquire  of  the 

choice  of  parameters  Z^,  Z^,  and  d to  achieve  both 
amplification  of  the  incident  wave  and  suppression  of  any 
standing  wave  components  downstream  of  the  two  loops. 

To  achieve  the  second  objective  we  demand  that 

A (z  > d)  ce  exp  - iQz/c  (Oil 

0 

which  requires 


Iq  + I exp(-  iQd/c)  = 0 


(012 


with  the  resultant  expression  for  A : 


a I 

A (z  > d)  = exp  (-  iQz/c)  1 - - — r2-  ( 1 - exp  (2iQd/c) 
o 4TTia  \ 


(013 


The  first  term  in  Eq.  (013)  represents  the  amplitude  of 
the  incident  cyclotron  wave;  the  second,  the  contributions 
from  the  two  loop  drives.  To  maximize  their  effect,  we 
take 


I 


exp  (2i  f.d/c ) 


1 


(014 


so  that  the  appropriate  spacing  of  the  loops  is  at  odd 


quater  wavelengths: 


d = — c/n,  N = 1,  3,  5,  7,  . . . 


(015) 


Using  Eqs.  (07)  through  (OlOb)  in  Eq.  (012)  we  find  that 


Z.  = Z + uua 
1 o 


(016) 


with  an  amplification  factor  A calculated  from  Eq.  (013): 


A = 1 + 


+ iuuL 


(017) 


Growth  Lengths  for  Discrete  Loop  Drives 


We  have  seen  that  a pair  of  passive  loops  located 


an  odd  quarter  cyclotron  wavelength  apart  operates  as  an 


amplifier  for  the  negative  energy  travelling  cyclotron 


wave  with  an  amplification  factor  given  by  Eq.  (017), 


provided  the  loop  impedances  are  chosen  consonant  with 


Eq.  (016). 


If,  for  the  first  pair  of  loop  drives  we  choose 


Zo  ifjuC 


(018) 


T-. 


where  the  capacitance  C is  chosen  so  that 


uu2  LC  = 1 
o 


(019) 


then  the  amplification  factor  is  infinite  (i.e.  we  obtain 
an  oscillator  or  wave  initiator) . Such  an  LC  circuit  could 
easily  be  pumped  up  to  high  power  levels  using  conventional 
RF  drives.  Furthermore,  the  total  power  dissipated  in  the 
first  loop  is  zero  (to  the  extent  that  the  loop  resistance 
is  zero);  i.e.,  only  reactive  power  is  required. 

The  second  loop  located  appropriately  downstream 
must  have  an  external  impedance 


Z = i ~ + uua 
1 l'juC 


(020) 


(i.e.  an  RC  series  network)  and  converts  the  initial 
standing  wave  set  up  by  the  first  loop  to  the  appropriate 
travelling  wave. 

Subsequent  pairs  of  loops  may  then  be  located 
downstream  of  the  first  pair  to  act  as  amplifiers . There 
is  no  requirement  on  the  spacing  of  the  pairs  of  loops 


relative  to  each  other,  other  than  the  fact  that  the 


spacing  must  be  sufficient  so  that  the  loops  are  not 
coupled  by  their  mutual  free  space  inductance. 

If  the  impedances  of  the  N™  loop  pair  are,  for 
example,  chosen  as 
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ituC 
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and  the  pairs  of  loops  are  spaced  roughly  every  cyclotron 
wavelength,  then  growth  lengths  £ of  the  order  of 

£ ~ 2rrc/n  (022; 


are  obtained. 

On  the  other  hand,  if  instead  of  operating  at 
resonance,  we  set  C = 0,  then  the  growth  length  increases 
to 


2ttc  L 
o a 


(023; 


Utilizing  the  value  of  a given  in  Eq.  (03)  , we  obtain 


and  note  that  this  growth  length  is  of  the  same  order  of 
magnitude  as  the  usual  type  resistive  liner  growth  length 
(Ref.  1) . 

Suppression  of  m = 1 Kink  Modes 

It  has  been  pointed  out*  that  the  high 
efficiency  and  strong  coupling  for  the  discrete,  resonant 
loop  drive  system  arise  from  the  fact  that  magnetic 
coupling  of  the  loops  to  the  electron  beam  is  being 
utilized,  as  opposed  to  electric  field  coupling  for  the 
usual  resistive  liner  configurations. 

Since  the  magnetic  field  of  a loop  drive  is  quite 
uniform  for  low  frequencies  (i.e.  the  field  configuration 
is  dominated  by  m = 0 contributions  for  jja/c  < 1)  , 
coupling  to  the  nonaxisymmetric  resistive  kink  mode 
should  not  occur  via  magnetic  coupling.  Coupling  can 

★ 

D.  Brower,  private  communication. 
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occur  through  the  electric  fields  as  in  the  usual  resistive 
liner  configurations;  however,  these  growth  lengths  would 
correspond  to  the  usual  resistive  liner  configurations 
and  would  be  correspondingly  much  longer  than  the  m = 0, 
magnetic  coupling  length. 

A second  attractive  feature  of  using  discrete  loop 
drives  is  the  stabilization  of  the  hydrodynamic  kink  mode. 
An  external  conductor  slotted  in  the  streaming  direction 
can  be  brought  arbitrarily  close  to  the  loop  drives 
without  affecting  their  operation.  This  configuration 
should  then  stabilize  the  hydrodynamic  kink  instability  by 
providing  a good  return  current  path  in  the  beam  direction 
(see  Appendix  I,  Ref.  1) . 


Some  Practical  Considerations  for  Resonant  Discrete  Loo 


is  straightforward . From  Eqs . (07)  and  (010a),  the  total 

current  in  loop  N = 0 is  calculated  to  be  given  by 


X _ i 2na  uj 

o (Z  + i;uL)c 
o 


with  total  gap  voltage  V given  by  V = Z I : 

G Goo 


V = - — . i 2rra  <ju/c 
G Z + ljjL 
■ o 


Relating  V to  the  total  potential  cp  of  the  wave 

through  the  relations  governing  the  cyclotron  wave 

(Qcp  ~ kxc  ~ 1kxc  A9  ~ for  1A^  |=1)#  we  obtain 


|VG/cp 


2rr  ijua  Q 


Z + iouL 
1 o 


where 


k,  = 2.9/a. 


Thus,  operation  at  resonance  and  low  value  of  resistance 
(e.g.  1/2  x-a)  in  order  to  obtain  high  amplification  factors 


per  pair  of  loops  may  result  in  a voltage  V,  exceeding 
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the  potential  of  the  cyclotron  wave.  For  large  amplitude 
cyclotron  wave  potentials,  it  may  then  become  necessary 
to  operate  considerably  off  resonance  to  avoid  high 
voltage  arc -over  problems.  Such  operation  will  of  course 
increase  the  total  growth  length  of  the  growth  section. 

A second  consideration  of  the  use  of  resonant 
discrete  loop  drives  is  the  required  ring -up  time  of  the 
loops.  Analysis  presented  in  this  appendix  is  for  the 
steady  state  behavior.  Care  must  be  taken  so  that  the  Q 
of  the  resonant  circuits  is  sufficiently  low  that  they 
can  reach  steady  state  early  in  the  electron  beam  pulse. 
The  ring  up  time  r is  given  by 


= L/R 


of  the  circuits.  Typically  L is  of  the  order  of  a few 
tens  of  nH;  and  for  amplication  factors  near  2,  we 
would  choose  R = 1/2  ijja,  a value  which  is  typically  a 
few  ohms.  Thus, 


t ~ 10-30  ns 


(029) 
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This  is  not  a severe  limitation  for  an  Auto-Resonant 
Accelerator  configuration,  since  e-beam  pulse  times  of 
100  to  200  ns  are  required  by  other  constraints  (i.e.  total 
ion  acceleration  time  requirements) . 
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APPENDIX  P 


DETAILS  OF  WAVE  ENERGY  BALANCE  IN  THE 
AUTO -RESONANT  ACCELERATOR 


: 

. 


Overall  energy  balance  in  the  Auto-Resonant 
Accelerator  is  governed  by  the  wave  energy  flux  relations 
and  the  conservation  equation  for  total  energy,  in  this 
appendix,  the  detailed,  radial  balance  of  wave  energy 
exchanges  between  the  electron  parallel  momentum,  electron 
perpendicular  momentum,  and  the  electromagnetic  field  is 
derived  for  a spatially  growing  negative  energy  cyclotron 
wave.  Integration  of  the  resulting  equations  over  the 
beam  results  in  the  usual  wave  energy  equation.  The  radial 
equations  demonstrate,  however,  that  the  energy  is 
extracted  not  from  the  center  of  the  beam,  but  rather 
near  the  outer  edge  of  the  beam  where  (k^r)  has  a 
maximum.  This  effect  is  important  when  considering  the 
virtual  cathode  formed  upon  injection  of  the  beam  into 
the  waveguide;  particles  near  the  center  of  the  beam  have 
relatively  little  free  energy,  so  that  if  energy  were 


required  to  be  extracted  from  these  particles,  total  field 
amplitudes  might  of  necessity  be  held  to  low  values. 

Before  proceeding,  it  is  perhaps  worthwhile  to 
recall  that  wave  energy  is  a quadratic  form  in  the  perturbed 
field  quantities  which  is  constructed  from  linear  theory 
and  which  can  be  shown  to  be  conserved  independent  of 
total  energy  conservation.  As  such,  wave  energy  is  calcu- 
lated using  only  products  of  first  order  quantities  (first 
order  in  the  electric  fields;  i.e.  those  quantities 
available  from  linear  theory)  and  therefore  does  not 
contain  zero  order  times  second  order  products.  It  is 
precisely  these  extra  terms  that  distinguish  total  energy 
conservation  from  that  of  wave  energy  and  renders  the  two 
equations  independent.  (It  is  also  true  that  in  some 
cases  wave  energy  and  total  energy  are  not  independent — 
e.g.f  nonrelativistic  plasmas  at  rest,  where  the  zero 
times  second  order  contributions  may  vanish  identically. 
However,  in  general,  wave  energy  and  total  energy  conser- 
vation do  yield  two  independent  constraint  equations.) 

With  this  point  in  mind,  we  now  proceed  to  a direct 
calculation  of  the  radial  wave  energy  balance  relations . 


In  the  following  treatment,  the  usual  high  gamma 
(ultra-relativistic)  limit  of  the  electron  equations  of 
motion  is  used,  so  that  we  consider  the  perturbed 
electron  parallel  velocity  - 0 to  lowest  order. 

We  consider  an  axisymmetric  negative  energy 
cyclotron  wave  with  the  basic  wave  structure: 

k± 

^ Az  k c - uj 
z 

iA  = A oe  j (k  r)  exp  i(k  z - uj)  ( 

9 r i x z 

with 

k * k - iX,  , 1 > 0 

z o 

Here  (k,  uj)  are  the  wave  vector  and  frequency?  with 
(A,  cp)  , the  usual  four  vector  potential  of  the  wave,  and 
c,  the  speed  of  light. 

For  relativistic  particles,  we  consider  the  expres 

sion 


<Ya>  = 1 + <pa>  + <pa> 


( 


where  p is  the  particle  momentum 


the  usual 


To  calculate  wave  energy,  we  must  pick  only  the  first  order 


The  calculation  of  (y)  can  be  effected  from  the 


equation 


iL/\  e/~  \ er-*-+*  -»  3A*  1 

^<Y>  --<E-V>  =-[-v7„  -V  • — ] 


Here,  again,  v,  cp,  and  A are  calculated  only  to  first 


order . 


Using  the  relation 


£<**> 


- v • Vcp*  + iuucp* 


we  obtain: 


^ <y  + e/m  cp*> 


e I'ju  * . 

v,  • A.  ) 

m c i J- 


Here  the  quantity 


<**>  - •' 


(P10) 


where  ftr  is  the  first  order  radial  variation  and  can  be 


related  to  the  perturbed  density  6n  by  the  equation  of 
continuity: 


7MWN 


(Pll) 
(P12) 

from  Poisson's  equation. 

Finally,  using  Eq.  (p4)  (noting  v^  = p±/y) , we 
obtain  the  following  expression  for  <y>: 


1- 

acp 

H 

3r 

n me2  I id  -k  c - Qla 
o z 1 


(P13) 


Here  we  have  assumed  tua  « kacs . 

We  see  that  indeed  the  second  order  changes  in  (y) 

are  proportional  to  J,  (k,r)  (i.e.  A,  , 3 m)  . Hence, 

lx  1 r 

energy  trams fer  occurs  near  the  edge  of  the  beam  and  not 
at  the  origin. 


Wave  Energy  Equation 


To  obtain  the  usual  conservation  equation  of  wave 


energy,  we  need  the  longitudinal  component  of  Poynting 


flux  S correct  to  second  order: 


c * * 

S = — (E  B„  - E B ) 
z 4rr  r 0 0 r 


Defining  the  quantity  W ; 

z 


W 3 n me2  <y>  + S 
* o z 


algebraic  manipulation  allows  W to  be  expressed  as 


- U)  v_ 


(kac3  + ou3)  |a  I3 

P r 1 


where 


2ou3  Q k c 

v < JSL  , c Tl  - 2 1 

g 3k  Lx  (k3c3  + tu3)3  J 

z p 


Using  the  relation 


(k3c3  + uu3)  3 k3 
P x 


E 3 as  k3  | cp  1 3 


we  then  obtain 
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Detailed  Energy  Balance 

It  is  of  interest  to  determine  the  amount  of  energy 
extracted  from  the  streaming  direction  and  utilized  to 
grow  and  sustain  the  cyclotron  wave.  Utilizing  Eq.  (P2) , 
we  can  compute  to  second  order: 


(p  > - (Vy2  - 1 - p2  (y  > - 


<P*> 


Utilizing  Eqs . (p6)  and  (pl3)  for  <p2>  and  <y), 
respectively,  we  obtain 


<PZ>  = 


e2  I A Is  (k2c2  + uj2)2  k2 
r p x 


uu2  ft2y  k3 

P 


X f-  i - *2  SI 

l U)2  k2 


- & £ tero)  } 


The  three  terms  in  the  braces  all  represent  energy  drains 
on  the  parallel  energy.  The  first  term  may  be  interpreted 
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NONLINEAR  FREQUENCY  SHIFT  OF  BEAM  CYCLOTRON  MODE 

In  the  linear  theory  of  the  electron  beam  cyclotron 
mode  (which  has  been  discussed  in  previous  reports) , the 
mode  amplitude  is  assumed  small.  For  finite  mode  ampli- 
tude, the  nonlinearities  may  no  longer  be  negligible^  and 
they  introduce  the  possibility  of  additional  effects 
associated  with  resonant  decay  instabilities,  nonlinear 
Landau  damping,  and  frequency  shifts.  The  mode  does  not 
appear  to  be  subject  to  resonant  decay  instabilities; 
and  for  a cold  electron  beam,  there  is  no  nonlinear  Landau 
damping.  In  this  report,  we  focus  our  attention  on  the 
nonlinear  modification  of  the  beam  cyclotron  mode  which 
is  manifested  in  the  appearance  of  a shift  in  the  wave 
frequency.  Such  nonlinear  frequency  shifts  may  influence 
the  time  evolution  of  the  beam  cyclotron  mode  especially 
in  situations  where  growth  of  the  mode  is  determined  by 
resonant  coupling  to  another  mode  of  the  system. 


i 


It  will  be  assumed  throughout  that  a fluid  descrip- 


tion of  the  electron  beam  is  adequate . In  evaluating  the 
nonlinear  frequency  shift,  the  procedure  is  to  calculate 
the  nonlinear  generation  of  zeroth  and  second  harmonic 
components  of  the  fundamental . The  nonlinear  interaction 
of  the  fundamental  with  these  harmonic  components  then 
leads  to  corrections  to  the  linear  dispersion  relation. 

A.  Equilibrium 

We  consider  a slab  model  equilibrium  in  which  a 
cold  relativistic  electron  beam  of  density  Nq  propa- 
gates with  velocity  uq  -♦  c along  a uniform  guide  field 

B in  the  z -direction.  The  beam  is  infinite  in  the 
oz 

y-direction  but  bounded  by  conducting  plates  at  x = ± a. 

The  electron  beam  dynamics  is  described  by  the 
cold  fluid  equations : 

3n 

— + V • n v = 0 
ot 

d „ e 

— m y v = - e E - — v x B 

at  — — c — — 

Y - (1  - v2/c2)-1/2  (Q  1) 
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f 


r 


perturbations  of  the  form  exp(-  i uot  + i k z)  , then 

o 


linearizing  Eq.  (Ql)  and  Eq.  (Q4),  we  obtain 


± le 


my  (uj  - k u ± Cl) 
o o o 


_af 

ax 


1 3 ± 


a A + H A 


= 0 


where 


±2 


- W ♦ - 3* 


3 ou_  (to  - k_u_) 


o o 


o c2  c2  (uj  - k u ± fi) 


o o 


4tt  N ez 


2 _ 

U)  = 

p my 


e B 


n = 


my  c 
o 


The  boundary  conditions  at  x = ± a are: 


(1)  A = a' 


(2>  + 1 r * - <Vc’>  v 


1 
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the  dispersion  relation 


If  k3  c2  > (ju3,  | K~a | < 1, 

o 

for  the  mode  symmetric  in  E is  given  by 


U)  = 


k u 
o o 


n k2c2 


(k2ca  + uu3) 


where 


k2  = k+2  + k3 
o 


+ + 
tan  h a = - h a 


(qs: 


The  corresponding  eigenfunctions  are; 


A+  (x)  = A sin  *+x 
o 


. + 

_ - , , . sin  h a , 

A (x)  = A sin  k x 

o sin  k a 


. . + x 

» A sin  k a • — 
o a 

In  order  to  calculate  the  nonlinear  frequency  shift, 
we  therefore  introduce  the  following  expansion  in  terms 
of  an  ordering  parameter  z: 


ini liv«  ivii1*' 


A+(r,t)  =ACo)±(r,t)  + e A(  1 5 ± (r,  t)  + 


where  A<0)±(r,t)  are  the  linear  solutions: 


ACo)±(r,t)  = A*(x,kQ)  exp(i  k^z) 


+ A~  (x,k  ) exp(-ik  z) 
o o o 


A+(x,k  ) = a e10  sin  k+x  s A sin  k+x 
o o o 


- , . . 18  x . + , X J + 

A (x,k  ) = a e sin  h.  a = A — sin  h a 

o o a o a 


where  the  amplitude  a and  phase  angle  9 are  given  by 


— = e D1(a)  + e2  D2(a)  + . - - 


Uy  _ 

=-u)  + s ©^a)  + s ©2  (a)  + . 


uu  is  determined  by  Eq.  (Q5),  and  k a » 2 for  the  mode 


with  the  smallest  perpendicular  wave  number. 
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n 


( o) 


e N k3  c3  2 ou3  . . 

o r/,  . P \ a + 1 a 


2m  y f23  (i)3 
o P 


Ul  + rr4)f  A+  - if  A-l 

L\  k3c3  /5x  o 2 Sx  o J 


,(0)  . . 


e k3x  + e a , + 


2m  na  o 2m  Qa  ax  o 


o) 


B(°) 

z 


2m 

na 

ic 

_a_ 

2Q 

ax 

c 

_a_ 

2n 

ax 

A - 


(A.  + A ) 


.(  o ) ± _ 


Ca  a2 


2Qa 


f*  K 

3xa  l o 


+ A~  l 
o I 


B 


< o)± 


ica  3s 


= ± 


. ± 


2na  ax3  ' o 


(A  + A ) T iA' 


where  we  have  neglected  terms  of  higher  order  in 

oua/kaca  < 1 and  1/y8  « 1.  For  convenience,  we  will  not 
P o 

always  denote  explicitly  the  arguments  x,  k , a,  8 of 

o 

the  dependent  variables . 
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where  H (x,k  ) is  determined  by  the  substitution  of  the 
o 

zero  order  and  first  order  solutions  in  the  nonlinear 
terms  of  Eq.  (Ql)  : 


+ * + 2 / 
H (X'V  " "in"  ( 


,|A  |*  A 


my  Q / uu4 
o p 


o 1 o 1 o 


h (x,k  ) = sin 
o 


H*X  {2  - 


sin  2 h*  a -i 
4k+  a I 


- an  3h+x 


K"*"  X 

+ ~p~  (flJJ'3x+X  - rQS  h+x) 


H (x,k  ) » 
o 


vC‘  / 

i6n  \ 


x+  2 / e \3  k6  c6 


n ) — 5 — 11  (x#k  ) |a  I2  A 

my  n / uu6  o 1 o 1 o 


h (x,k  ) = 3 an  3k*x  - sin  k+x 
o 


we  have  neglected  terms  of  higher  order  in  uu2/k2c2  < 1 

P 

and  1/y2  « 1 . 
o 


...  . . ; 


4ireN 


4-rTeN 


r 


■ * 


• — jppgnn 


We  solve  Eq.  (Q12)  for  v{3)±(x,k  ) and  substitute  in 

o 


Eq.  (Q13)  and  Eq.  (Q14).  we  multiply  Eq.  (Q13)  by 


A (x,k  ) and  Eq.  (Q14)  by  A (x,k  ) , and  after  inte- 
o o o o 


grating  from  x = - a to  x = + a,  we  sum  the  resulting 
equations : 


[-  ®2  + 


^2  ■ 

la  + -L  i 

e 

a . 

J gi  4n  ^ 

^ mv  n 

o 

(Q15) 


where 


gx  = 1 + ca?  H+a 


+a 


gQ  = h+  J dx  st  h*x  |~h+  (x)  - 


(sin  3it*x  + sin  x+x)  i 
2 


x*  sin  K*  a 


dx  ^ 1 h (x)  + 2^3  sin  3x+x  - sin  x*x^ 


-a 


For  x+a  « 2 


g ~ 3.73  , g,  « 1.17 
o 1 


1 
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Equating  real  and  imaginary  parts  of  Eq.  (Q15) 


D 

2 


= 0 


0. 


6(ju  = 


2-2 


o ksc 

4H+ag.  uu2 
1 P 


ex.+A 
c 

my  Q! 
o 


n 


Eq.  (Q16)  yields  the  shift  6uj  in  the  frequency  due  to 
nonlinearities . 

The  conditions  for  validity  of  the  expansion 
utilized  in  deriving  Eq.  (Q16)  are 


1 

2 


k3c: 


(JO 


en+A 


my 

o 


< 1 


and 


im  < 


k2c3 


a 
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This  calculation  does  not  include  changes  in  the 
time  averaged  variables  which  are  produced  by  a time 
variation  in  the  amplitude  of  perturbation.  It  is  clear 
that  such  changes  cannot  be  neglected  if  we  are  to  properly 
discuss  the  effects  of  nonlinear  frequency  shifts  on  the 
evolution  of  an  unstable  wave  growing  from  small  amplitudes. 
We  are  continuing  the  analysis  to  determine  the  role 
of  these  time  averaged  changes  induced  by  growing  amplitude 


perturbations . 
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NONLINEAR  EFFECTS  FOR  THE  ELECTRON 
BEAM  CYCLOTRON  MODE 

Much  of  the  previous  analysis  of  the  Auto-Resonant 
Accelerator  has  been  based  upon  the  linear  theory  of  the 
electron  cyclotron  mode,  which  is  valid  when  the  mode 
amplitude  is  sufficiently  small.  For  example,  the  detailed 
investigation  of  various  means  of  growing  the  cyclotron 
eigenmode  has  been  conducted  within  linearized  theory, 
as  has  the  analysis  of  ion  loading  and  acceleration.  How- 
ever, it  is  of  interest  to  determine  the  mode  amplitude 
level  at  which  nonlinear  effects  first  become  important 
for  a number  of  reasons : 

(1)  The  successful  functioning  of  the  accelerator 
for  nonlinear  wave  amplitudes  has  not  yet  been 
theoretically  established  (or  disproven) . 

(2)  Nonlinear  effects  may  well  cause  saturation 
of  the  linear  growth  mechanisms,  establishing 

an  upper  bound  to  the  achievable  mode  amplitude. 


(3)  Theoretical  knowledge  of  such  nonlinear 
saturation  amplitudes  and  mechanisms  is 
needed  for  comparison  with  computational 
and  experimental  investigations* 

(4)  Knowledge  of  the  magnitude  and  parameter 
scaling  of  nonlinear  mode  amplitudes  is  a 
valuable  design  tool,  since  this  information 
bears  heavily  on  the  required  length  of  both 
the  growth  and  acceleration  sections  of  the 
accelerator . 

In  this  appendix,  a brief  survey  of  a number  of 
nonlinear  effects  is  presented,  and  rough  estimates  are 
given  of  the  magnitude  and  scaling  of  the  mode  amplitude 
level  at  which  these  effects  become  important.  Many  of 
these  effects  may  be  seen  to  be  related  to  one  another. 
Although  care  has  been  taken  to  preserve  the  numerical 


coefficients  for  both  cartesian  and  cylindrical  geometry, 
the  very  nature  of  these  formulas  is  merely  to  indicate 
the  level  at  which  linear  approximations  lose  their  validity, 
and  this  occurs  gradually  rather  than  sharply  as  a function 
of  amplitude. 


215 


Although  it  is  possible  in  some  cases  to  construct 
nonlinear  theory  which  retains  its  validity  to  the  point 
of  growth  saturation,  so  that  the  saturation  amplitudes 
are  known  precisely,  such  investigations  are  not  described 
herein.  Consequently,  the  correlation  between  growth 
saturation  and  the  breakdown  of  linear  approximations 
must  remain  inferential,  pending  the  accumulation  of 
additional  theoretical  or  computational  evidence. 

Finally,  the  effects  discussed  herein  are  confined 
to  the  electron-wave  interaction;  nonlinearities  associ- 
ated with  ion  dynamics  or  with  a wave -growth  inducing 
medium  or  additional  driving  fields  are  not  considered. 

The  symbols  which  are  used  are  defined  below; 

e,m  = electron  charge,  mass 

0qc  = DC  electron  axial  flow  velocity 

y = (1  - 02)-1/a  = DC  electron  v-f actor 

O O 

y = perturbed  electron  y-factor 

nQ,fi  * DC,  perturbed  components  of  electron  density 
a,b  = electron  beam,  wall  radius 

x^  = perturbed  radial  (perpendicular)  electron 

coordinate 
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UN*1-  W*.  J.P«  i KWU>. J — 


•]  1 

9,-= 

= perturbed  radial  (perpendicular)  electron 

! 

velocity 

1 1 • 

8,c 

» perturbed  axial  electron  velocity 

B 

oz 

= external  magnetic  guide  field 

kl'kZ 

= perpendicular,  axial  components  of  wave 

vector 

k3 

= k3  + k3 
■>•  z 

“p 

= (4tt  n e3/y  m)l/3  = relativistic  electron 

O 0 

plasma  frequency 

n 

= (e  B /y  me)  = relativistic  electron 

OZ  o 

cyclotron  frequency 

6(ju 

= f)  uu3/k2c2 
P 

1 

U) 

= frequency  of  electron  cyclotron  wave 

r 

= linear  growth  rate  of  electron  cyclotron 

wave 

iJ;NL 

= nonlinear  frequency  shift  of  electron 

s 

H 

' 

cyclotron  wave 

EOi(xi> 

= perpendicular  component  of  DC  electric 

* 

field  at  radius  x^ 

. 

i , (xj 

= perpendicular  component  of  wave  electric 

field  at  radius  xx 

• 
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Sz  (XJ.> 


*.  (xl> 


$(xx) 


axial  component  of  wave  electric  field 
at  radius  x 

x 

DC  electric  potential  function 

amplitude  of  perturbed  electric  potential 

1,  cartesian 

geometrical  parameter  = (slab)  geometry 

2,  cylindrical 

geometry 


In  addition,  the  above  perturbed  or  wave  quantities  refer 


to  the  magnitude  of  the  amplitude  only;  phase  information 
is  supressed. 


The  analysis  which  underlies  these  conditions  of 
nonlinearity  is  based  in  part  upon  the  following 
approximations : 


b - a 
a 


« 1 


(Rl) 


Y » 1.  0 - 1,  0 « 1 

o z 


(R2) 


axisymmetry,  with 


= xx>  (R3) 


uua  « k2c2 

z 


(R4) 
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L .A. 


OB 


k3  + T~a 
2 2 c3 


« k3,  k a 
± z 


2**1 

C 


(R5) 


Condition  (Rl)  allows  use  of  the  analytically  simple, 
one-region  eigenmode.  Condition  (R2)  implies  strong 
decoupling  of  the  longitudinal  plasma  mode  from  the 
cyclotron  mode.  Condition  (R3)  is  a reminder  that  we 
have  considered  only  the  symmetric  eigenfunctions.  Con- 
dition (R4)  for  nonrelativistic  wave  speed  is  well 
satisfied  for  the  proposed  proof  of  principle  Auto- 
Resonant  Accelerator  experiment,  and  allows  simplifica- 
tion of  the  cyclotron  eigenfunctions  by  decoupling  the 
light  modes.  Condition  (R5)  implies  strong  decoupling 
of  the  positive  and  negative  energy  cyclotron  modes,  and 
also  allows  simplification  of  the  cyclotron  eigenfunctions, 
Subject  to  these  approximations,  one  finds  the 
following  radial  mode  structures: 


9>(xx)  = 5(o)  Tjo  (kxr)  , ajyfkjXj] 


(R6) 


Ex(xa)  * kx5(o)  [ji  (kxr)  , sin  (kxxx)] 


(R7) 
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k^a  =»  1 + G (for  the  broadest  eigenmode) 


* (xx>  “■  2lTe  no  XI/G  = 2 xx  Eoj.(xj.) 


(RIO) 


(Rll) 


where  the  Bessel  functions  apply  for  cylindrical  geometry 
and  the  trigonometric  functions  for  slab  geometry. 

The  various  conditions  wherein  nonlinear  effects 
become  significant  are  listed  below  for  ready  comparison 
with  one  another;  a discussion  of  the  different  nonlinear 


effects  follows. 
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Condition  (R12)  prescribes  the  wave  amplitude  level 

at  which  the  electron  density  fluctuations  due  to  the  wave 

become  comparable  with  the  mean  electron  density.  These 

density  fluctuations  will  be  strongest  on  axis,  and  at  this 

wave  amplitude  the  electron  density  will  vary  from  near 

zero  to  around  twice  the  ambient  density. 

Condition  (R12)  also  indicates  that  these  density 

nonlinearities  will  become  important  at  an  amplitude  level 

smaller  by  yu2/k2c2  from  that  for  which  longitudinal 
P 

trapping  of  the  electrons  by  the  wave  potential  would  occur. 
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Condition  (R13)  indicates  three  nonlinear  effects 
where  each  occurs  at  a wave  amplitude  higher  than  that  for 
density  nonlinearities  by  only  the  geometrical  factor. 

The  first  of  these  is  the  condition  that  the  radial  dis- 
placement of  an  electron  by  the  wave  be  as  large  as  the 
scale  length  kj1  characterizing  the  perpendicular 
eigenmode  structure.  This  is  judged  to  be  more  accurate 
than  the  criterion  for  displacements  as  large  as  the  beam 
radius,  which  requires  higher  wave  amplitudes.  Equation  (R9) 
indicates  that  the  displacements  will  be  largest  at  about 
0.8  times  the  beam  radius.  One  physical  consequence  of 
radial  displacements  as  large  as  k7*  a/(l  + G)  is  that 
wall  collisions  might  occur  unless  a sufficiently  wide 
vacuum  gap  was  available.  The  second  condition  in  Eq.  9rl3) 
is  that  finite  electron  gyroradius  effects  become  important 
because  of  the  perpendicular  kinetic  energy  excited  by  the 
wave.  The  third  condition  is  that  the  radial  motion  of 
the  electrons  becomes  nonlaminar  by  virtue  of  radial 
"crossings"  of  electron  trajectories.  This  may  be  easily 
seen  to  occur  when  | 3x1/3xx | exceeds  1,  which  occurs 
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first  for  particles  on  axis.  These  last  two  conditions 
in  Eq.  (R13 ) are  symptomatic  of  the  breakdown  of  the  fluid 
nature  of  the  wave  due  to  nonlinear  turbulence. 

Condition  (R14)  prescribes  when  the  perpendicular 
wave  electric  field  exceeds  the  self  radial  field  of  the 
beam,  which  tends  to  occur  first  on  axis.  This  requires  a 
wave  amplitude  larger  than  that  for  density  nonlinearity  by 
only  k2/k2.  Condition  (R14)  also  reveals  the  connection 
between  these  various  criteria  of  nonlinearity  and  the 

measurable  ratios  such  as  E /E  (a)  which  are  con- 

iinax  o-l 

veniently  extracted  from  computer  simulations  ;s  a 
diagnostic.  One  physical  consequence  of  the  nonlinearity 
suggested  by  Eq.  (R14)  is  the  possibility  of  electron  field 

emission  at  the  conductor  walls.  The  self-field  E is 

o 1 

typically  much  larger  than  the  threshold  amplitude  required 
for  field  emission,  but  is  in  the  wrong  direction.  One  might 
think  that  if  the  wave  fields  became  larger  than  the  self- 
field and  had  the  opposite  radial  sense,  then  field  emission 
might  result.  However,  this  is  a deception  of  the  linear 
theory  which  is  breaking  down,  since  Poisson’s  equation 
shows  that  even  for  nonlinear  waves,  the  electrostatic 
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portion  of  the  electric  field  must  always  point  toward  a 
region  of  net  negative  charge.  Hence,  the  only  possibility 
of  field  emission  is  due  to  the  inductive  or  electromagnetic 
portion  of  the  electric  field,  which  is  small  for  nonrela- 
tivistic  phase  velocity,  or  to  the  presence  of  the  ions 
whose  density  might  locally  exceed  the  nonlinearly  reduced 
electron  density.  To  quantify  this  criterion  for  field 
emission  would  require  a detailed  analysis  of  the  ion 
dynamics,  which  is  not  presented  herein.  However,  one 
suspects  that  the  required  ion  density  would  be  a signifi- 
cant fraction  of  the  electron  density,  since  the  electron 
wave  fields  are  geometrically  reduced  below  their  peak 
value  at  the  walls,  and  the  peak  value  does  not  exceed 
the  self-field  until  (n/r^  ) ^ k2a2/k^a  according  to  Eq.  (R14). 

Condition  (R15)  for  the  radial  electron  oscilla- 
tions in  the  wave  to  become  relativistic  will  be  stronger 

than  condition  (R13)  only  if  k exceeds  k . It  is 

z *- 

possible  that  this  might  be  modestly  true  in  the  growth 

section  and  at  the  beginning  of  the  accelerator  section, 

but  toward  the  end  of  the  accelerator  section  k will  be 

z 

below  k.  and  B,  <<  1. 
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The  significance  of  the  0^1  nonlinearity  is 

that  it  implies  a significant  fluctuation  y in  the 

electron  Y-factor  as  shown  in  condition  (R16) . This 

occurs  because  the  electrons  are  perturbed  by  the  wave  to 

move  with  or  against  the  radial  self -field  of  the  beam.  At 

the  beginning  of  the  accelerator  where  3 is  largest, 

E (a)/B  is  smallest  and  may  be  somewhat  less  than  1. 

The  net  product  of  these  terms  yields  an  expression,  also 

shown  in  Eq.  (R16),  which  indicates  that  Y / yq  is  generally 

somewhat  less  than  n/n  . 

o 

However,  condition  (R17)  indicates  that  the 

variation  y/y  is  magnified  by  k2c3/u,2  to  produce  a 
o p 

variation  in  the  fluid  electron  cyclotron  resonance 

|uu-k  c0  + n|  which  may  be  comparable  to  its  ambient 
z z 

value  (i.e.,  6uu)  . This  nonlinear  spreading  of  the 


cyclotron  resonance  occurs  due  to  the  y -dependence  of  Zlj 
it  is  sometimes  called  "resonance  broadening."  It  is  also 
symptomatic  of  the  breakdown  of  fluid  theory  and  the 
conversion  of  the  cyclotron  mode  to  a kinetic  regime  due 
to  nonlinear  turbulence.  Condition  (R17)  indicates  that 
this  nonlinear  resonance  broadening  effect  is  felt  at 
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wave  amplitudes  somewhat  below  those  for  density  non- 
linearity, and  hence  is  the  most  sensitive  nonlinear  effect 
yet  discussed.  Experience  would  suggest  that  this  resonance 
broadening  would  reduce  the  wave  growth  rate  when  the 
wave  amplitude  reaches  the  level  given  by  Eq.  (R17). 

Condition  (R18)  gives  the  corresponding  resonance 
broadening  which  occurs  due  to  changes  in  the  axial 
electron  velocity.  The  first  terms  of  Eq.  (R18) , which  are 
linear  in  the  wave  amplitude  n/nQ,  are  relatively  unim- 
portant as  usual  for  y » 1.  However,  the  last  term 

o 

of  Eq.  (R18) , which  is  quadratic  in  the  wave  amplitude,  may 

be  seen  to  be  more  (or  less)  important  than  the  y/y 

resonance  broadening,  depending  on  whether 

ksG3c2/2k2uu2a2  is  greater  (or  less)  than  1.  This 
z p 

quadratic  term  in  the  wave  amplitudes  reflects  the  per- 
turbation in  0 due  to  linear  perturbations  in  0 , 
according  to 

1 ~a 

0 - — 0 + linear  terms 

z 2 ± 

+ other  quadratic  terms  (R20) 
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which  follows  immediately  from  a perturbation  expansion  of 
the  definition  of  y(P) • For  the  parameters  which  have 
been  proposed  for  the  proof  of  principle  experiment  on 
Auto-Resonant  acceleration,  assuming  k ~ k at  the 

Z -L 

beginning  of  the  accelerator,  one  may  see  that  conditions 
(R17)  and  (R18)  "trigger"  at  roughly  the  same  wave 
ampl itude . 

With  the  exception  of  the  final  term  in  Eq.  (R18) , the 
effect  of  nonlinearities  has  only  been  explored  so  far 
through  the  retention  of  terms  which  are  linear  in  the 
wave  amplitude.  In  an  attempt  to  examine  nonlinear 
effects  which  are  quadratic  in  the  wave  amplitude,  a 
calculation  is  begun  in  Appendix  Q which  seeks  the  non- 
linear shift  in  the  frequency  of  the  electron  cyclotron 
wave  due  to  nonlinear  coupling  to  harmonics  of  the  linear 
eigenmode.  Although  this  calculation  is  presently  in- 
complete as  mentioned  in  Appendix  Q,  the  result  will 
scale  generally  as  indicated  in  condition  (R19) . The 
frequency  shift  will  be  quadratic  in  the  wave  amplitude, 
and  it  is  believed  that  the  coefficient  function  H is 
at  least  order  unity.  It  is  speculated  that  when  this 
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nonlinear  frequency  shift  exceeds  the  linear  growth  rate 

r of  the  cyclotron  wave,  the  growth  will  be  strongly 

modified.  For  growth  mechanisms  such  that  T « 6uu, 

which  is  the  regime  of  validity  of  Eq.  (R19) , one  sees  that 

this  nonlinear  effect  may  appear  much  sooner  than  the 

density  nonlinearities  and  even  sooner  than  resonance 

broadening  effects.  In  the  competition  between  Eq.  (RT7)  and 

(R19) , a low  value  of  T favors  (R19)  while  a high  value 

of  k a favors  Eq.  (R17).  Of  course,  if  k a is  too  large, 
z z 

the  cyclotron  eigenfunctions  should  be  reevaluated  due 
to  the  violation  of  assumption  (R5) . 

It  is  important  to  notice  that  the  mode  amplitude 
at  which  these  nonlinear  effects  appear  is  reduced  for 
higher  values  of  ka . Consequently,  one  would  expect 
the  higher-k  modes  to  saturate  in  growth  at  lower  ampli- 
tudes than  the  lower -k  modes.  When  exciting  or  growing 
the  cyclotron  wave,  care  should  be  taken  (as  in  designing 
computer  simulations)  that  the  lowest  k mode  is  not 
discriminated  against  by  any  frequency  or  wavenumber 
selection  which  may  be  imposed. 
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It  should  not  be  inferred  that  this  list  of 
electron  nonlinearities  is  complete,  or  that  there  are 
not  other  important  nonlinear  effects.  For  example,  the 
possibility  of  a nonlinear  three-wave  decay  instability 
was  previously  examined  and  found  to  be  disallowed  by  the 
coupling  selection  rules . There  are  also  a number  of  ion 
nonlinear  effects  which  could  produce  saturation  of 
electron-ion  instabilities,  for  example,  at  wave  ampli- 
tudes which  may  be  below  those  discussed  here  in  some 


cases . 


APPENDIX  S 


COMPUTER  SIMULATION  OF  THE  RESISTIVE  GROWTH  AND 
SATURATION  OF  THE  ELECTRON  CYCLOTRON  MODE 

Explanation  of  Code 

The  simulations  were  performed  with  a two  dimensional 
simulation  code  fully  described  in  Ref.  6.  The  code  is  a 
two  dimensional  Cartesian,  relativistic,  fully  electro- 
magnetic particle  code. 

The  code  has  been  modified  as  described  in  Ref.  7 
to  eliminate  the  numerical  Cherenkov  instability. 

The  primary  additional  modification  made  to  the  code 
as  previously  described  is  the  inclusion  of  a nonisotropic, 
homogeneous  resistive  medium  in  the  code.  This  is  implemented 
by  modifying  the  x component  of  the  current  as  shown 


6.  Drobot,  Adam  T. , "Theory  and  Experiments  in  the 
Simulation  of  Collisionless  Plasmas"  (unpublished  disserta- 
tion, Fusion  Research  Center,  Univ.  of  Texas  at  Austin,  1974). 

7.  Drummond,  W.  E.,  et  al.  "Turbulent  Heating  of 
Plasmas  by  High  Energy  Beam,"  Interim  Scientific  Report, 
AFOSR  Contract  No.  F44620-72-C-0071.  Arlington,  Virginia: 
Air  Force  Office  of  Scientific  Research,  Aug.  1975. 
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J - J „ + aE  (SI) 

X px  X 

where  J is  the  x component  of  the  particle  current 
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computed  in  the  normal  manner  by  bilinear  interpolation. 

The  aE  addition  to  the  current  will  simulate 
x 

the  action  of  a nonisotropic  resistive  medium  distributed 
throughout  the  volume  of  the  beam.  This  correction  is 
made  at  every  mesh  point  at  which  the  current  is 

defined.  The  purpose  of  the  correction  is  to  render  the 
negative  energy  electron  cyclotron  mode  unstable.  The 
linear  theory  and  dispersion  relations  for  this  case  are 
worked  out  in  Appendix  T. 

The  motivation  for  simulating  the  distributed 
resistivity  case  rather  than  having  the  resistivity 
concentrated  in  the  wall  is  computational  convenience. 


An  additional  motivation  is  an  enhanced  growth  rate.  Thi; 


will  become  evident  when  one  considers  the  complications 
involved  in  simulating  a resistive  liner. 

In  order  to  solve  the  fields  in  the  liner,  it  is 
necessary  to  either  implement  an  appropriate  jump  condi- 
tion at  the  boundary  or  else  solve  for  the  fields  across 
a transition  region  which  is  of  order  c/a  wide.  This 
would  imply  a very  tiny  Eulerian  cell  width  and  con- 
sequently a large  number  of  cells.  Implementing  jump 
conditions  at  the  boundary  would  involve  putting  in 
special  logic  to  determine  where  the  boundary  lay,  and 
then  inserting  a jump  in  the  appropriate  quantity. 

Either  one  of  these  would  require  a significant  amount 
of  reprograming. 

Another  reason  for  using  distributed  resistivity 
is  that  the  growth  rate  turns  out  to  be  higher  than  if 
the  resistivity  is  concentrated  in  the  walls.  The 
reason  for  this  is  the  E^  field  is  a monotonic  decreasing 
function  of  the  radial  coordinate  in  the  walls,  it  is 
therefore  evident  that  more  energy  will  be  dissipated 
the  further  in  the  resistive  liner  extends.  Because  of 
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the  negative  energy  character  wave,  this  results  in  a 
higher  growth  rate. 


It  may  also  be  expected  that  the  nonlinear  satura- 
tion mechanism  for  the  distributed  resistivity  model  will 
probably  give  valuable  insight  to  the  nonlinear  satura- 
tion of  the  resistive  liner  configuration. 

The  format  of  the  runs  was  to  allow  the  unstable 
wave  to  grow  from  noise.  We  could  then  determine  at  what 
level  saturation  occurred  and  if  any  other  unstable  waves 
would  grow.  Letting  the  wave  grow  from  noise  dictates 
long  computer  runs  since  the  wave  must  grow  by  four  to 
five  orders  of  magnitude  before  nonlinear  effects  might 
become  important  and  also  because  it  typically  takes  1500 
to  2000  time  steps  before  the  wave  becomes  discernible 
from  noise. 

The  first  simulation  run  has  the  following 
parameters . 

n At  = 0.3 
o 


I 


where 

fi  = e B / me  c 
o x 

4TTe  n3 
e 

u)  = 

po  m 

e 

e = electronic  charge 

= applied  guide  field 

m = electron  mass 

e 

c = velocity  of  light 

c At/Ay  = 0.3  50 

c At/ Ax  = 0.477 

At  = time  step 

Ax  = longitudinal  cell  size 

Ay  = transverse  cell  size 

Nx  = 40 

= number  of  active  cells  in  longitudinal 
direction 

Ny  = 13 

= number  of  active  cells  in  transverse 
direction. 
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Figure  SI  shows  that  the  most  unstable  mode 


corresponds  to  one  wavelength  in  the  machine  length.  The 
graph  represents  the  amplitude  of  the  longitudinal  E 
field  as  a function  of  longitudinal  position.  By  comparing 
such  graphs  produced  at  different  points  in  time,  it  is 
possible  to  calculate  a phase  velocity  of  approximately 
0.1c. 

Figure  S2  shows  E^  plotted  as  a function  of  time. 

The  growth  rate  associated  with  the  amplitude  of  E is 

x 

3.01  x lO-3!^  as  shown  in  Figure  S2.  This  corresponds 
to  a growth  rate  of  379  time  steps  per  e -folding. 

The  linear  theory  for  this  case  presented  in 
Appendix  T gives  a growth  rate  of 
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the  theoretical  e-folding  time  turns  out  to  be  429  time 
steps  per  e-folding.  This  constitutes  reasonable 
theoretical  agreement. 

Figure  S2  shows  that  the  amplitude  of  the  E^ 
field  takes  a dip  close  to  the  end  of  the  simulation. 

The  total  energy  (particle  plus  field)  increases  by  about 
4 parts  in  a thousand  during  this  same  interval . 
total  energy  had  been  constant  to  within  one  part  in  104 
until  this  point  in  the  run,  it  is  probable  that  this  dip 
is  due  to  loss  of  numerical  energy  conservation. 

The  maximum  field  amplitude  reaches  1/18  of 

the  radial  E field  amplitude  at  its  peak.  This  is  probably 
too  low  for  nonlinear  effects  to  become  important. 
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The  plasma  parameters  of  the  second  and  third 
simulations  are  the  same. 


Q a 
c 


1.85 


a = beam  radius 
c = velocity  of  light 

The  first  simulation  was  carried  out  on  a 15  by  42 
mesh  with  3174  particles.  The  beam  was  initially  filled 
half  way  to  the  wall.  The  simulation  was  run  a total  of 
10500  time  steps.  A plot  of  energy  versus  time  is 

shown  in  Figure  S3 . It  yields  an  amplitude  growth  rate 

T = 3.56  x 10-3  n 

The  approximate  growth  rate  from  the  linear  theory  is 

T = 1.86  x lo“3  ft 

o 

or  approximately  half  the  observed  growth  rate.  This  con 
stitutes  reasonable  agreement  since  some  of  the  approxima 
tions  made  in  the  theory  are  not  well  satisfied  in  the 
simulation. 
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ENERGY- ARBITRARY  UNITS 
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The  energy  plot  in  Figure  S3  shows  a clear 

saturation  where  the  growth  of  the  wave  stops,  the 
amplitude  declines  slightly  and  then  equilibrates.  The 
amplitude  of  field  at  this  point  is  approximately 

8.5%  of  the  D.C.  E field. 

y 

A brief  theoretical  survey  of  the  various  nonlinear 
effects  which  might  begin  to  cause  growth  saturation  has 
been  performed  for  this  computer  run.  In  examining  the 
behavior  of  the  linear  wave  eigenfunctions,  some  previous 
results  were  utilized  which  employed  a number  of  approxi- 
mations including  the  following: 

b - a 

« 1 (b,  a = wall,  beam  radius) 

d 

v ~a  « l 

uoa  <<  ka  ca  « ka  ca 

k ||  a s*  1 (S4) 

Although  these  approximations  are  poorly  satisfied  for  the 
parameters  of  this  run,  it  is  suspected  that  the  theoretical 
result  is  sufficiently  insensitive  to  the  precise 
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eigenmode  structure  that  the  results  of  the  calculation 
may  remain  valid  within  a factor  of  perhaps  3. 

One  of  the  first  nonlinear  effects  that  appears  to 
be  encountered  for  this  run  is  that  fluctuations  in  the 
electron  y -factor,  which  occur  as  the  electrons  are 
radially  perturbed  to  travel  through  the  DC  radial  electric 
field,  may  effectively  broaden  the  electron  cyclotron 
resonance  |uu-k||  cP  + q|  beyond  the  bounds  of  cold  fluid 
theory  when  the  wave  field  reaches  the  level 


(ktia)  (k±a) 


'xDC 


k4  a4 


0.04 


( S 5 ) 


In  addition,  the  perpendicular  component  of  the  electron 
velocity  may  become  relativistic  at  a field  amplitude  a 
few  times  larger  than  this . 

Although  there  is  no  certainty  that  these  nonlinear 
effects  will  cause  wave  growth  saturation,  it  is  possible 
that  saturation  will  result,  and  this  provides  an  ampli- 
tude level  which  may  be  compared  with  the  observed 
computer  calculation  results . 


243 


f 


it  *•  W ' 1 _ __ 


The  momentum  of  the  particles  drops  by  approximately 
9%  during  the  time  period  from  4800  time  steps  to  8500  time 
steps . This  indicates  that  the  wave  energy  is  being 
extracted  from  the  streaming  energy  of  the  electrons . 

A plot  of  total  energy  versus  time  is  shown  in  Figure 
S4 . The  total  loss  of  energy  is  approximately  5%  of  the 
initial  energy.  Some  of  this  energy  loss  may  be  accounted 
for  by  energy  dissipation  in  the  resistive  media.  The 
rest  of  the  loss  must  be  accounted  for  in  terms  of  trunca- 
tion errors  made  by  the  code. 

The  third  simulation  had  the  same  plasma  parameters 
as  the  second  run.  However,  the  conducting  walls  were 
brought  in  much  closer  to  the  beam.  The  initial  fill  had 
the  beam  radius  equal  to  0.928  times  the  wall  radius.  This 
run  corresponds  more  closely  to  the  linear  theory  since  the 
theory  assumed  the  beam  radius  approximately  equal  to  the 
wall  radius. 

Figure  S5  shows  energy  as  a function  of  time. 

The  amplitude  growth  rate  for  this  case  is 

» 

r = 4.66  x lO"3  Qq  (S6) 
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Figure  S4.  Total  Energy  vs.  Time 
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Figure  S5.  E Energy  vs.  Time-Run  2 
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The  corresponding  linear  growth  rate  is 

T = 2.03  x 10"3  n (S7) 

o 

The  energy  conservation  for  this  run  was  within  0.7%  for 
the  first  4800  time  steps  and  within  2.6%  for  the  first 
5600  time  steps.  The  amplitude  of  the  field  is  8.6% 

of  the  D.C.  transverse  field  at  4800  time  steps  and  22.4% 
of  the  D.C.  transverse  field  at  5600  time  steps. 

The  three  simulation  runs  tend  to  indicate  that  E 

x 

can  reach  approximately  15%  of  the  transverse  D.C.  electric 
field  before  saturation.  This  is  with  no  effort  being 
made  to  maximize  saturation  level.  It  must  also  be 
pointed  out  that  energy  conservation  became  a serious 
problem  towards  the  end  of  these  runs . The  validity  of 
these  results  can  be  trusted  only  as  long  as  the  energy 


conservation  remains  good. 


Extra  diagnostics  have  been  added  to  the 
code  which  involve  doing  longitudinal  fast  Fourier  trans- 
forms on  selected  field  and  current  components  at 
selected  transverse  locations.  These  show  both  phase  and 
magnitude  of  the  selected  Fourier  components  so  that  we 
can  deduce  phase  velocity  and  uu  of  the  various  modes  we 
wish  to  analyze.  This  package  involves  using  the  FFT 


routines  in  the  IMSL  Library. 


i 
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APPENDIX  T 

GROWTH  OF  THE  CYCLOTRON  WAVE  BY  AN 
ANISOTROPIC  RESISTIVE  BACKGROUND 


Introduction 

In  this  appendix  we  derive  the  dispersion  relation 
for  a relativistic  slab  of  electrons  propagating  inside 
a perfectly  conducting  waveguide  in  which  there  is 
a guide  magnetic  field  and  an  anisotropic  resistive 
plasma  which  obeys 

J = 0 = J 
px  py 

J = c E (Tl) 

pz  z 

where  J = (J  , J , J ) is  the  plasma  current,  E is 
p px  py  pz  z 

the  component  of  the  electric  field  in  the  + z (propaga- 
tion) direction,  and  is  the  (constant)  plasma  con- 

ductivity. From  this  dispersion  relation  we  will  obtain 
estimates  of  the  growth  rate  of  the  symmetric , lower 
Doppler  shifted  cyclotron  eigenmode. 
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Equilibrium  Solution 


We  choose  rectangular  Cartesian  coordinates  (x,  y, 

z)  with  the  beam  propagating  in  the  + z direction  with 

velocity  v =*  c , where  c is  the  velocity  of  light  in 
z 

vacuum.  The  beam  is  assumed  to  extend  from  x = +a  to 
x = - a in  x and  is  assumed  to  be  of  infinite  extent  in 
y and  z.  The  guide  magnetic  field  is  given  by  = 

A 

B z with  B constant.  We  shall  also  assume  that  the 
o o 

beam  is  extremely  relativistic,  i.e.  y » 1 where 
y = (1  - v2/c2 ) _1/2  . 

For  our  equilibrium  solution  we  take  the  "rigid 
rotor"  solution  to  the  Maxwell -Lorentz  equation,  which 
is  (c  = 1) 


E=4rrriex  B = 4tt  n e v (T2) 

x o y o z 


E=E=0  B = 0 (T3) 

y z x 


B = B 
z o 


cu4  (a3  - Xs  ) 

_E 

il2  y2 


B 

O 


(T4 ) 
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where  J is  the  beam  current  density  in  kiloamps , 

KA 

B is  the  guide  field  strength  in  kilogauss  and  a 
KG  cm 

is  the  beam  half  thickness  in  centimeters. 

Perturbed  Solutions 

We  consider  perturbation  of  the  form  f (x)  exp 
(st  + kz)  which  are  solutions  to  the  linearized 
relativistic  pressureless  fluid  equations  and  Maxwell 
equation  with  zero  order  fields  given  by  Eqs.  (T2)  through 
(T6).  Note  that  we  take  all  quantities  to  be  independent 
of  y . 

Neglecting  the  equilibrium  v^  and  the  diamagnetic 

term  in  B , the  fluid  equations  give 
z 

(s  + ik)  v = — (E  - B ) + Q v (T12) 

x my  x y y 


(s  + ik)  v = — (E  + B ) - Q v 
y my  y x x 


(T13 ) 


v = 0 (T14) 

z 

for  the  perturbed  velocity  field  v - (v  , v , v ) in 

— x y z 

terms  of  the  perturbed  fields  E and  5.  The  Maxwell 


equations  for  E and  B are 


~ , 4tt  lk  - 

7 • E = 4tt  p - <7  E 

— s z 


7 • B = 0 


7 • E — - sB  7 • B = sE  + 4tt  J + 4n  c E z (T15) 

z 


where 


J = e ri  v + e-n2v,  o = e ri 
o ' — ' 'o 


and  we  have  used  the  continuity  equation  for  the  resistive 
plasma 


so  + 7 • J = so  + ik  (j  E = 0 
P P P z 


(T16) 


We  now  define  two  functions  A+  and  A as 


)-7 


A 3 (l  + — ) (i  ± iE 
± \ s / \ x y 


Equations  (T15)  then  imply 


1 + A - 2(l  + - 2 5 

+ - \ s/xsx 


S E 
X z 


A -A  = 2i  E 11  + 


ik 


(T17) 


(T18 ) 


(T19) 
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VA+  + A-> 


2 (s  + ik)  (l  + J Ez 


(T20) 


Hi 


o that  if  A+  and  A are  known,  E may  be  determined 
and  B may  then  be  found  from  the  curl  E equation. 

Thus  the  solution  is  completely  determined  by  A+  and  A . 

Now  using  the  linearized  fluid  equations  (T12) 
through  (T14) , together  with  the  Maxwell  equation  (T15), 
we  find  that  A___  must  satisfy 


2ttc 

s 


52  A 


(T21) 


2ttj  _ 
s 


5i  A+  - ( 


2tt? 

s 


+ l)  ?2 


(T22) 


where 


>u2  (s  + ik) 

pa  = _ s2  - k2  - 2 (T23) 

(s  + ik)  ± iQ 

The  solution  to  Eqs.  (T21)  and  (T22)  is 


A_  = A1  sin  q+x  + B ^ cos  3+x  + A2  sin  q_x  + B2  cos  q_x 

(T24 ) 


254 


A+  = - A a_  sin  q+x  - Bl  a _ cos  q+x 


-A  a+  sin  q_x  - B>2  a+  cos  q_x 


(T25) 


where 


2b  ?2  a±  s - (b  + 1)  (§*  - §») 


± -y / (b  + 1)  2 (5*  - ?2)  + 4b2  ?2  ?2  (T26 ) 


2 q2  = (b  + 1)  (52  + §2) 


± W(b  + 1)2(?2  - ?2)2  + 4b2  *2  ?2  (T27) 


(T28 ) 


Equations  (T24)  through  (T28)  with  appropriate 
boundary  conditions  determine  the  solution  to  the 


linearized  perturbed  equations. 
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Dispersion  Relation 

We  assume  that  the  perfectly  conducting  waveguide 
is  located  at  x = ± a so  that  the  appropriate  boundary 
conditions  are 

E = E = 0 at  x = ± a 
z y 

which  from  Eqs.  (T19)  and  (T20)  gives 
A+  = A at  x = ± a 


(A+  +A^)=0  at  x = ± a 


These  four  conditions  give  four  linear  homogeneous 
equations  for  A^  , B , and  B2 . Setting  the  determi- 
nant of  the  coefficients  equal  to  zero  then  gives  the 
dispersion  relation  which  is 


a+  + 1 


a +1 


a. 


T cot  V - — i ■*- 


q a cot  q a 


X 


a+  + 1 


a, 


~ q a tan  q a 
1 + ^+ 


a +1 


a - 1 - 


q a tan  q a 


= 0 
(T29) 


•1 


'I 
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We  note  that  had  we  taken  B^  = B^  = 0 initially 
in  Eqs.  (T24)and  (T25)  the  dispersion  relation  would 
have  only  contained  the  first  factor.  Similarly,  setting 
= 0 gives  the  second  factor  as  the  dispersion 
relation . 

Growth  Rate  of  the  Symmetric  Mode 

For  the  symmetric  mode,  E (x)  = E (-x),  so  that 

z z 

we  require  B,  = B_  = 0 in  the  solutions  for  A . We 
12  ± 

thus  have  for  this  mode  the  dispersion  relation 
a+  + 1 a - 1 

q a cot  q a - — q a cot  q a = 0 (T30 

a-l+  + a - 1 - - 


For  the  lower  Doppler  shifted  cyclotron  mode  we 
have  (for  low  phase  velocity  waves,  yu/k  « 1 ) 


u)  =■  k - a 


(T31 


Writing  jl>  = -ju  + i T with  m and 


purely  real  we 


r " |k2+  ?2  + u)2  |s  *-R  ?- 


(T33 ) 


when  § and  ? are  the  real  and  imaginary  parts  of 
? respectively.  Thus  to  obtain  the  growth  rate  T we 
must  solve  Eq.  (T30)  for  C . 

Using  Eqs.  (T26)  through  (T28)  we  obtain 


+ 1+b  - (1+b)  ?2/§2  - 2b?  /? 

t — *T*  — 


a +1 
a ~ 


(2b + 1)  ?2 

[i  ““  b(?2  + ?2)  + e2 


(2b + i)  §2 

2 

l+  ~ (b  + 1)  - b?2/?2 


for  fjj  « k.  The  dispersion  relation  thus  becomes  (taking 


a = 1)  , 


f 
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q+  cot  q+ 


/,  4rrrr  \ 

= - l.1  * v ) q-  cot  9< 


(T38) 


with 


(T39) 


(2b  + 1)  §2 

b?3  + 


(T40) 


In  Eqs.  (T38)  through  (T40)  we  have  assumed  that 


S3/?3  « l 


(T41) 


Now  for  c = 0 (b  = 0)  the  dispersion  relation  Eq. 


(T38)  becomes 


q+  cot  q+  = - q cot  q 


(T42) 


and  for  ik  =-  q+  small  in  modulus,  this  implies 


q (ff  = 0)  = 2 


(T43 ) 


Expanding  q about  the  solution  (T43)  of  (T42)  we 


obtain 
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q 


where 


2 (4rrq)  2 

" a (Da  + (4rrq)  2 
o 


ia 


4rr  'ju  q 

o 

a)2  + (4iTq)  2 
o 


(T44) 


a 


cot 


cot  q^ 

q - q esc  2 q^ 


q=0 


0.32 


Using  Eq.  (T44),  Eq.  (T40)  may  be  solved  for  £2  for 
use  in  evaluating  the  growth  rate.  The  result  is 


4 (1  - 2 ib) 
1 + 4b2 


so  that  we  finally  obtain 


(T45) 


r =» 


4fl  uu2 

E 

+ 4 + u)a  Is 
P 


4rr  uu  q 

o 

uu2  + (4rrq)2 
o 


(T46) 


Note  that  Eq.  (T44)  implies  |52|  < k2  for  k £ 1 as  assumed 
in  obtaining  Eqs.  (T38)  through  (T40) . 

Now  the  second  factor  in  Eq.  (T46)  maximizes  at  a value 
of  0.5  for  ou  =» 


4rrq,  while  for  low  phase  velocity  waves 


APPENDIX  U 


A MODEL  OF  ION  PICKUP  AND  ACCELERATION 
BY  A TRAVELLING  WAVE 

The  design  which  has  been  formulated  for  the 
proof -of -principle  experiment  on  Auto-Resonant  accelera- 
tion contemplates  a puff -gas  injector  section  (Ref.  1) , 
wherein  the  subject  ions  will  be  created  by  primary  beam 
ionization  of  a neutral  gas,  then  trapped  or  picked  up 
by  the  preformed  travelling  cyclotron  wave  which  is 
supported  by  the  electron  beam,  and  finally  carried  out 
of  the  puff -gas  section  into  the  accelerator  section 
where  the  wave  and  ions  are  accelerated  above  the  initial 
wave  velocity.  Consequently  it  is  of  interest  to  determine 
the  degree  of  precision  with  which  the  wave  amplitude  and 
phase  velocity  must  be  controlled  in  order  to  achieve  an 
acceptable  pickup  efficiency. 


In  addition,  since  the  puff-gas  injector  section 


will  be  between  the  wave  growth  and  the  wave  acceleration 
sections,  it  is  desirable  to  generalize  the  analysis  to 
include  the  possibility  that  the  wave  is  still  growing,  or 
has  begun  accelerating,  when  the  ions  are  first  intro- 
duced. Such  a theoretical  framework  would  also  allow  the 
consideration  of  alternative  design  options.  For  example, 
one  might  envision  extending  the  wall  flare  of  the 
acceleration  section  back  into  the  puff-gas  section,  so 
that  an  axial  DC  electric  field  would  be  present  there  to 
help  eject  the  ions  out  into  the  accelerator  section. 

Then  one  could  conceivably  grow  and/or  accelerate  the 
wave  as  necessary  to  pick  up  the  ions  in  the  beginning  of 
the  accelerator  section. 

In  order  to  address  these  questions,  a simple 
one -dimens ional , nonrelativistic  model  of  the  ion-wave 
interaction  is  explored.  It  is  hypothesized  that  the 
ion  equation  of  motion  in  the  axial  direction  is 

t 

z = | E<fc)  sin  k [z  - J dr  v#  (t)  j (U1) 

o 


where  k 


mm 


= axial  wave  number,  assumed  constant 
e,M  = ion  charge,  mass 

E (t)  - amplitude  of  the  axial  wave  electric  field 

v#  (t)  = axial  wave  phase  velocity 

and  E (t)  and  v^  (t)  are  presumed  to  be  controlled 
externally,  independent  of  the  behavior  of  the  ions. 

One  seeks  the  ion  motion,  as  it  evolves  from  prescribed 
initial  conditions  subject  to  the  prescribed  functions 
E(t)  , v (t)  . 

It  is  convenient  to  make  a formal  transformation 
of  variables: 

t 

y h 2 - J dr  v#  (t)  + tt  (U2) 

o 

whence  one  obtains  the  transformed  equation  of  motion, 

y = - ^ E(t)  sin  ky  - v (t)  (U3) 

Here  y may  be  viewed  as  the  ion  coordinate  relative  to 
the  travelling  wave  frame  of  reference.  One  may  define  a 
Hamiltonian  function 
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H(y,y,t)  = | y2  _ cos  ky  + y v#  (t) 


which  generates  the  equations  of  motion  BH/dy  = y. 


3H/By  = - y so  that 


^H(y.  S',  t) 


e E (t) 


ky  + y v (t)  (U5) 


Therefore,  when  the  wave  amplitude  E and  .the  wave 
acceleration  A = v\  are  constant  in  time,  it  follows 
that  H = 0 so  that  H(y,  y)  is  a constant  of  the  motion. 
This  allows  an  immediate  determination  of  the  ion 
trajectory. 

Furthermore,  when  E and  A are  not  constant, 
but  vary  only  slowly  in  time  such  that  E and  v#  are 
small,  it  is  possible  to  define  a new  quantity 


J(H,t)  = £ dy  y[H,  y,  t] 


which  is  aaiabatically  invariant  and  thus  serves  as  an 
approximate  constant  of  the  motion  to  yield  the  ion 
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^ ***!*!*!'»'  v , | 


trajectories,  in  Eq.  (U6),  y(H,  y,  t)  is  obtained  by 
inverting  Eq.  (U4)  and  t held  fixed  while  <£  dy  is 
performed --between  turning  points  or  over  a wave  period. 

For  the  case  where  E and  A are  constant  in 
time,  it  is  instructive  to  examine  the  function 


eE 

v(y)  = - ^ ky  + Ay 


where  H(y,yj  = — y2  + V (y) . This  potential  function 
V(y),  which  is  a ramp-biased  sinusoid,  is  depicted  in 
Figure  Ul.  It  has  relative  maxima  a bit  to  the  right  of 
odd  multiples  of  n,  and  relative  minima  a bit  to  the 
left  of  even  multiples  of  rr . When  E and  A are 
constant,  H = -j  y3  + V(y)  is  a constant  of  the  motion, 
so  that  a particle  trajectory  is  characterized  by  a 
horizontal  line  on  Figure  ui  with  the  ordinate  H.  The 
particle  motion  is  confined  to  that  region  of  y-space 
where  H £ V (y) , with  turning  points  y^  where 
V(y  ) = H.  Given  the  initial  conditions  (y,  y)  after 
which  E,  A are  constant,  the  subsequent  motion  may  be 
inferred  immediately  from  inspection  of  H,  V on 


Figure  Ul . 


Figure  Ul.  Potential  Function  of  Accelerating  Wave 
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When  A > 0,  all  trajectories  will  have  either  a 
rightmost  turning  point  only  (i.e.,  untrapped  particles) 
or  else  two  turning  points  (i.e.,  trapped  particles). 

With  two  turning  points,  the  particles  are  confined  to 
a one -wavelength  region  of  y-space  and  hence  move  with 
the  wave  relative  to  the  laboratory  coordinate  system 
(i.e.,  z) . With  only  a rightmost  turning  point,  the 
particle  is  reflected  at  most  once,  after  which  it 
travels  continually  to  the  left  relative  to  the  wave, 
with  increasing  relative  speed.  in  the  laboratory  frame 
of  reference,  such  a particle  would  be  untrapped  and  the 
wave  would  be  accelerating  away  from  it.  When  A < 0, 
the  situation  is  symmetric,  with  untrapped  particles 
exiting  the  wave  in  the  forward  direction. 

It  may  be  observed  that  when 

I A | > — (U8) 

then  dv/dy  has  the  sign  of  A and  V(y)  is  monotonic 
increasing  (decreasing)  if  the  acceleration  is  positive 
(negative) . In  this  case,  there  exist  no  particles 
trapped  in  the  wave . 


268 


It  is  instructive  to  consider  the  following  series 
of  problems,  which  provide  some  insight  as  to  the  proper- 
ties of  wave  amplitude  and  acceleration  which  are  commensu- 
rate with  a high  ion  pickup  efficiency: 
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(a)  Initial  conditions — particles  are  distributed  uniformly 

in  z;  |y|  = |z  - v | = u for  all 
particles  ? 

dynamics — the  field  amplitude  is  grown 
suddenly  from  0 to  E,  while 
f = A = 0 . 

9 

(b)  I nitial  conditions — same  as  for  problem  (a)  ; 

dynamics — the  field  amplitude  is  grown  slowly 
(adiabatically)  from  0 to  E,  while 
v . = A = 0 . 

(c)  Initial  conditions — same  as  for  problem  (a) ; 

dynamics — same  as  for  problem  (a),  followed 
by  wave  acceleration  at  v.  = A, 

9 

while  E = constant. 

(d)  Initial  conditions — same  as  final  conditions  of 

problem  (b) ; 

dynamics --the  wave  is  accelerated  at 

v,  = A,  while  E = constant. 

9 


Problem  (a) 

This  problem  has  been  considered  previously  (Ref.  1, 
p.  308)  in  computing  the  pickup  efficiency  of  a wave 


- - 


2 


travelling  at  constant  amplitude  through  the  puff  gas 


section.  The  sudden  approximation,  which  requires 


fast  electron  velocity,  and  is  at  full  strength  when  an 


ion  is 


(U14) 


which  is  a curve  shifted  upward  from  V(y  ) by 


Inserting  Eq.  (U14)  into  Eq.  (U12)  allows  calculation  of 


particles  in  the  range  (Ull)  , those  are  trapped  for  which 


(U15) 


so  that  the  pickup  efficiency  is 


s . . = — Arc  cos 

pickup  tt  w 


kMua 

4eE 


(U16) 


For  this  sudden  growth  problem,  there  will  always  be  some 

untrapped  particles  near  y = ± n/k,  so  e <1. 

o pickup 

However,  the  larger  the  field  amplitude,  the  larger  the 

pickup  efficiency.  Marginal  pickup  (i.e.,  e . £ 0) 

pickup 

occurs  for 


— u2  <,  2 — = V(y  ) - V(v  ) 

2 Mk  v/peak'  ^valley' 


(U17) 


and  there  is  no  pickup  at  all  for  smaller  values  of  E/u2 
which  violate  Eq,  (U17) 

Problem  (b) 

This  problem  has  also  been  considered  previously 
(Ref.  8) . The  adiabatic  invariant  is  found  to  be 


9.  See  for  example,  J.  R.  Thompson,  Phys . Fluids 
14,  1532  (1971)  . 
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J(K.t) 


4 / eE  (t) 
tt  V Mk 


|kE2  (h-1  ) S(*  - 1) 


+ [e2(x)  - (1  - x2 ) K (h)  ]s  (1  - x)  } 

(U18) 


where 


X 


+ 


3 1 

eE  (t)  ./Mk  ) 


K,E2  = complete  elliptic  integrals  of  the  first, 
second  kind 


S (x)  = — + — sign  (x)  = unit  step  function 


and 


2tt  E (t) 
ku  E (t) 


(U19) 


is  required  to  validate  the  assumption  of  adiabat icity . 
Since  J is  an  adiabatic  invariant  for  all  particles, 
and  they  all  have  the  same  value  of  H = j u2  initially, 
it  follows  that  all  particles  will  have  the  same 
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H (t)  !”and  x.(t)j  at  future  times,  if  the  amplitude  growth 
is  halted  when  E(t)  = E,  with  the  acceleration  A 
remaining  zero,  the  final  value  of  H or  x.  may  be 
determined  from  the  J-invariance  equation. 


J(H,t) 


J(x,t) 


J 


initial 


u 


(U20) 


From  Figure  U1  with  v.  = A = 0,  or  from  Eq.  (U4)  one  may 
note  that  marginal  trapping  occurs  for 

x - 1 
eE 

H ~ nr  (U21) 

Mk 

Particles  remain  untrapped  for  larger  values  of  k,  H, 
and  the  particle  trajectories  in  (y,y)  phase  space- 
determined  from  Eq.  (U4)  with  H = constant  are  unbounded 
in  y.  Likewise,  for  - (eE/Mk) <H< (eE/Mk)  or  0 < x.  < 1, 
the  particles  are  trapped  and  their  trajectories  are 
closed  bounded  contours  in  phase  space.  From  Eqs.  (U18)  , (U20) 
and  (U21)  , one  finds  that  marginal  trapping  occurs  when 


274 


T 


— 


. -I.S'.  n.A 


eE  _ rrf_ 
Mk  ~ 16 


and  the  pickup  efficiency  is 


(U22 ) 


_ ST  eE  - — 1 
'pickup  l Mk  u3  16  J 


(U23) 


The  pickup  efficiencies  of  Eqs.  (U16)  and  (U23)  are  plotted  in 
Figure  U 2 to  illustrate  the  difference  between  sudden  and 
adiabatic  growth  in  the  field  amplitude  required  for 
pickup . 

Problem  (c) 

The  analysis  of  this  problem  can  proceed  exactly 

as  for  problem  (a) , except  for  the  minor  complication 

that  A ^ 0.  Equation  (U14)  again  describes  the  initial 

conditions  for  H(y  ) , with  V (y  ) given  by  Eq;  (U7)  and 

o o 

illustrated  on  Figure  Ul  . Condition  (U12)  defines  which 

particles  will  be  trapped  in  the  accelerating  wave,  and 

it  has  the  simple  graphical  interpretation  of  selecting 

those  particles  for  which  the  H(y  ) contour  falls  below 

o 
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the  leftmost  peak  of  V (y  ) which  occurs  at  y = y , 

o o peak 

given  by  Eq.  (U9) . without  loss  of  generality,  consideration 
may  be  restricted  to  the  coordinate  range  Eq.  (Ull) , since  the 
trapping  condition  (U12)  will  be  unaltered  for  any  other 
wavelength  range. 

Condition  (U17)  of  problem  (a)  for  marginal  pickup 
has  the  following  obvious  geometrical  generalization  when 
A ji  0: 


^ u 3 s V - V 

2 peak  valley 


where 


G ( a)  = J 1 - a2  + a Arc  sin  a-  j a 


(U24) 


(U25) 


and  the  equality  in  Eq.  (U24)  implies  e . » 0 while  the 

"pickup 

inequality  in  Eq.  (U24)  is  necessary  for  ? >0.  Since 

"pickup 

G(a)  falls  from  1 to  0 as  a = AM/eE  increases  from  0 
to  1,  it  is  obvious  that  larger  accelerations  reduce  the 
allowable  speed  of  particles  that  may  be  picked  up  with  a 
given  wave  amplitude,  and  condition  (U8)  continues  to  be 
necessary  for  any  pickup  whatever. 
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Furthermore,  since  G < 1,  Eq.  (U24)  is  more  difficult 
to  satisfy  than  Eq.  (U17) . This  means  that  if  no  particles  are 
trapped  prior  to  wave  acceleration,  none  can  be  trapped 
afterward  without  additional  growth  in  the  wave  amplitude. 
That  is,  if  the  particles  are  initially  travelling  faster 
than  the  wave,  and  the  wave  is  slowly  accelerated  to  a 
higher  velocity  than  the  particles  start  with,  then  no 
particles  can  be  picked  up  by  the  wave — no  matter  how 
slowly  the  acceleration  is  accomplished.  .The  wave  will 
accelerate  right  through  the  mean  particle  speed  without 
picking  up  any  particles . 

Given  that  condition  (U24)  is  satisfied,  then  the 
trapped  particle  coordinates  may  again  be  determined  by 
inserting  Eq.  (U14)  into  Eq.  (U12)  to  obtain 


ci  hy  - res  ky  = ^ 1 - a2  + a Arc  sir,  a - na  - ~ 


Mku‘ 


eE 


(U26) 


where  a = AM/eE.  Equation  (U26)  will  have  two  solutions 


for  y,  which  we  designate  y . and  v , where 

mm  max 


2v 


y < y , < y and  y < y < y + 

peak  mm  valley  valley  max  peak  k 


1 
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Those  particles  initially  in  the  range  gq#  (uil)  will  be 

trapped  for  which  y . £ y , £ y , so  that 

min  o trapped  max 

the  pickup  efficiency  is  given  by 


epickup 


k(y  - y . ) 
max  mm 


(U27) 


As  in  problem  (a) , because  of  the  sudded  growth  in  ampli- 
tude, there  will  always  be  some  untrapped  particles  near 

y y , no  matter  how  large  eE/Mku2  is , although 
o peak 

the  number  of  such  untrapped  particles  approaches  zero  as 
eE/Mku2  increases.  In  addition,  because  of  the  finite 
wave  acceleration,  there  will  always  be  some  untrapped 


particles  near  y s;  y ( equivalently  y ^ y , + — 

o peak  \ o peak  k 

in  range  (Ull)  \ , and  the  number  of  these  untrapped 
particles  remains  finite  even  when  Mku2/eE  approaches 
zero.  One  may  show  that  the  limiting  fraction  of  untrapped 
particles  obeys 


a (2tt  - 2tt  e 


lost 


sin  2tt  ) = 2 /l  — o.2^  sin2  t 

lost  V lost 


(U28) 


whose  bounding  solutions  are 
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a « 1 


'lost 


./a/" 


a =*■  1 


where 


e . , s l - e 
pickup  lost 


(U29) 


(U30) 


. and  the  equality  in  Eq.  (U30)  applies  when  Mku2/eE  approaches 
zero.  It  may  be  seen  that  eloat  is  monotonic  increasing 
with  a#  while  epickup  is  monotonic  decreasing  with  a, 
so  that  wave  acceleration  always  has  a negative  influence 
on  particle  trapping  for  this  problem. 

Problem  (d) 

For  this  problem,  the  wave  acceleration  begins  when 
the  adiabatic  growth  ceases . Suppose  that  this  initial 
phase  of  adiabatic  growth  results  in  the  particle's 
trapping  parameter  k being  reduced  to  some  value  h 

o 

when  the  wave  growth  is  terminated.  This  value  h = 

o 

x.o(eE/Mkus)  may  be  determined  from  Eqs.  (U18)  and  (U20) , and 


— — 'f.wp,wf  - 


the  particles  will  be  trapped  prior  to  wave  acceleration 

if  x < 1. 
o 

The  wave  acceleration  phase  of  this  problem  may  be 
examined  as  before  via  the  potential  function  V(y)  in  Eq. 
(U7)  . However,  instead  of  Eq.  (U14) , the  initial  conditions 
for  this  problem  are  that 


H - if  (2  »o  - X)  + A yc 


(U31) 


which  is  a linear  ramp  rather  than  a curve  paralleling 

V (y  ) . Furthermore,  if  x < 1 the  particles  do  not 
o o 

exist  throughout  the  entire  range  Eq.  (Jll) , but  are  restricted 
to  the  range 


h Arc  cos  (1  - 2x=)  < y < f;  Arc  cos  (1  - 2x2) 

K O O K O 


(U32) 


modulo  2n/k. 

For  the  case  in  which  x > 1,  the  possibility  of 
subsequent  trapping  may  be  examined  by  inserting  Eq.  (U31) into 
Eq.  (U12) . However,  even  for  the  most  favorable  point  y^  = 
ypeak'  °nS  ^n<*s  that  Eq.  (U12)  reduces  to 
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(U33) 


which  is  impossible  for  hq  > 1.  Therefore,  no  particles 
can  be  trapped  after  acceleration  if  they  were  not  trapped 
before  acceleration. 

For  the  case  in  which  * <1,  there  is  the 

o 

possibility  of  complete  trapping  due  to  the  restricted 

range  Eq.  (U32)  of  particles  confined  initially  below  the 

lip  of  the  potential  wells . Complete  trapping  will  occur 

H[y  = k_1  Arc  cos  (1  - 2 xs)  ]<  v (y  which  reduces 

o o peak 

to 


a Arc  cos  i 1 - 2 k2) 
o 


(1 


2 H*)  s 
o 


+ a Arc  sin  a - rra 


(U34) 


For  example,  if  hq  “ 2_1/a  , complete  trapping  requires 
a 0.217. 

When  k is  less  than  1 but  greater  yhan  the 
o 

value  for  which  Eq,  (U34)  is  marginally  obeyed,  then  some 


particles  will  be  spilled  by  the  wave  acceleration.  The 


pickup  efficiency  may  then  be  computed  as  before  by 


inserting  Eg.  (U31)  into  Eg.  (U12)  , and  solving  for  the  initial 
coordinates  of  marginal  trapping,  subject  to  the  range 
Eg.  (U32). 


In  order  to  illustrate  these  ideas  graphically,  as 
well  as  to  examine  situations  where  E(t)  and  ^ (t)  are 
arbitrarily  varying  in  time,  a computer  code  was  written 
to  solve  Eg.  (Ul)  numerically.  Examples  of  problems 
(a)  through  (d)  have  been  computed,  including  the  following 


Run  1 


0, 


particles 


- 0?  0 < 


Zparticles< 


2ir/k 


flit) 

Mku2 


. 0.04  kut 

0.01  e 


0 < t < 146/ku 


3 .438 


146  ^ 300 

» *: — < t < - — 


V#  (t) 
u 


The  phase  from  0 < t < 146/ku  is  the  adiabatic  growth  of 
problem  (b) . The  particles  first  become  trapped  at 
t » 103/ku,  and  when  t ■ 146/ku,  the  trapping  parameter 


0 < t < 


1+0.7  (kut  - 200) 


is  reduced  to  h =“  2-1/a  . The  phase  from  (200/ku)  < t < 
(300/ku)  is  the  wave  acceleration  of  problem  (d) . The 
acceleration  parameter  is  a = MA/eE  =-  0.204,  which  is 
small  enough  to  satisfy  (U34)  for  complete  pickup,  in 
the  computer  run,  each  of  20  test  particles  was  picked  up 
and  accelerated?  the  results  are  illustrated  in  Figure  U3 , 
where  the  mean  particle  velocity  is  compared  to  the  wave 
phase  velocity. 

Run  2 


0.1 


0 < t < 300/ku 


eE(t) 

Mku8 


v#  (t) 

— - — = -1+0.02  kut,  0 < t < 300/ku 


This  is  an  example  of  problem  (c) , where  a = 0.^  and 
eE/Mku2  is  so  small  that  Eq.  (U24)  is  violated  and  no 
particles  can  be  picked  up.  This  was  confirmed  in  the 
computer  run,  whose  results  are  illustrated  in  Figure  U4. 


Run  3 


@ t « 0, 


particles 


= 0, 


< z , < 

particles 


2rr/k 


g-E.(t) 

Mku2 


0.5 


, 0 < t < 300/ku 


v#  (t) 

— - — * -1+0.02  kut,  0 < t < 300/ku 

In  this  case,  eE/Mku3  is  large  enough  to  obey  Eq.  (U24) . 

From  Eqs.  (U26)  and  (U27) , one  finds  that  « . , “■  0.46 

pickup 

theoretically.  In  the  computer  run,  9 of  20  test 


particles  were  picked  up  and  accelerated;  the  results  are 


also  illustrated  in  Figure  U4  . 


In  conclusion,  this  one  dimensional  model  of  wave 
growth  and  acceleration  reveals  that  the  primary  require- 
ment for  particle  pickup  is  that  wave  growth  should  occur 
before  the  wave  is  accelerated  relative  to  the  particles . 
This  growth  must  proceed  beyond  the  point  of  marginal 
trapping,  and  then  the  wave  acceleration  must  be  gentle 
or  else  a large  fraction  of  the  particles  may  be  dumped 
from  the  potential  wells.  Complete  pickup  is  possible 
only  if  the  wave  growth  is  adiabatic  and  the  acceleration 
is  gentle. 
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APPENDIX  V 


VACUUM  REQUIREMENTS 


Operation  of  the  Auto-Resonant  Accelerator  depends 
initially  on  the  establishment  of  a rigid  rotor  equilib- 
rium in  which  the  net  radial  expansion  forces  due  to  the 
electron  beams  self -fields  and  centrifugal  effects  are 


balanced  by  a v0B  restoring  term.  Derivation  of  the 

U Z 


self -field  force  term  has  generally  assumed  that  there  is 
no  space  charge  neutralization  due  to  ionization  of  ambient 
gas  atoms.  If  a uniform  fractional  neutralization, 
n+/nk  = f,  is  assumed  then  self-field  force  term 


becomes 


m (Er  - W 
o 


(JU 


po 


, 1 'N  r 

U3  " f)  2 


(VI) 


It  is  clear  that  for  neutralization  to  be  an  unimportant 
effect  it  is  necessary  that 


f <<  rr 


(V2) 


The  rate  of  neutralization  build-up  due  to  ionization  of 
neutral  ambient  ions  by  primary  beam  electrons  is  given  by 


L. 
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I 


— ' 


' 


df 

dt 


(n  - n)  a v =“  n a c 
a + z a 


(V3 ) 


where 


n^  = ambient  gas  density 

a = ionization  cross-section  =“  10~18cm2 

(assuming  oxygen  and  nitrogen  are  the 
dominant  atomic  species) . 


Also 


n 


3-5  X 1016P 


(V4) 


where 


Using 


P = ambient  pressure  in  torr. 

cL 


y =“  7 


then  using  Eqs . (V2)  and  (V4)  there  obtains 


r 

- 


t 'I 

6*.  • 

: : 


J. 

1 


B. 

1 

B „ 


E . 
1 


vi/a 


(Wl) 


where 


E^  = initial  ion  energy  =-  0.25  MeV 


E = final  ion  energy  =30  MeV 


The  energy  fluence  incident  on  the  beam  stop  is  then 


O = J-  E T 
f e 


where 


E^  = electron  kinetic  energy  =3.0  MeV 


(W2) 


= beam  pulse  length  =“  3 00  ns 


Since 


then 


J - 9.55  kA/cm2 


cp  =*■  785  j/cm2 
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Using  Spencers  (Ref.  9 ) calculation  for  the  energy  depo- 


sition profile  produced  by  electrons  in  low  z materials. 


the  ratio  of  peak  dose  to  incident  fluence  is 


J/gm 


E dx 
o 


j/cm2 


Then 


D =-  712  J/gm  = 170  cal/gm 


For  most  elements  this  amount  of  deposited  energy  is 


sufficient  to  melt  the  material.  Further,  the  short 


deposition  time  of  the  energy  will  produce  thermo- 


mechanical stresses  which  could  produce  spalling  of  the 


material . 


Based  upon  its  low  Gruneisen  coefficient  and  high 


dose  level  to  produce  melting,  graphite  is  chosen  as  the 


best  beam  dump  material.  Specifically,  ATJ  graphite  is 


recommended.  The  thickness  required  to  completely  stop 


the  electron  beam  is  (Ref.  9) 


9.  Spencer,  L.  B.,  "Energy  Dissipation  by  Fast  Electrons," 
N.B.S.  Monograph  No.  1,  1959. 
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t =-  1.7  gms/cm3 


Using 


p = 2.0  (±  10%)  gms/cm3 


a thickness  of  0.85  cm  (0.354  inch)  is  sufficient. 


Radiation  Production 


The  0°  bremsstrahlung  radiation  yield  from  a gold 
convertor  bombarded  by  3.0  MeV  electrons  is  (Ref. 10) 


Y = 2.5  x 104  Rads/coulomb  @ 1 meter 


Since  the  yield  scales  linearly  with  atomic  charge  (Z) , 
the  value  for  carbon  is  a factor  of  13  lower.  The  total 
charge  per  pulse  incident  on  the  beam  dump  is 


Q = I T 

o 


Assuming  the  machine  fires  once  every  5 minutes  (At  = 

1/12  hour) , the  dose  rate  at  1 meter  (assuming  no  shield- 
ing) would  be 


Martin,  T.  H.,  "Determination  of  Bremsstrahlung 
Production  Efficiencies  from  Data  Obtained  on  PHERMEX 
at  27  MeV,"  AFWL  TR  73-65,  1973. 
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D = ~r  y 1 t/ At 
13  O 


2.25  X 105  ^ 
Hr 


To  reduce  this  to  a safe  level  of  roughly  2 mR/Hr,  an 
attenuation  of  10s  is  required.  For  3.0  MeV  bremsstrah- 
lung,  attenuation  is  dominated  by  Compton  scattering  so 
that  for  all  materials  the  10th  value  thickness  is 
60  gms/cm2.  Thus,  a total  material  thickness  of 
300  gms/cm2  is  required.  Since  concrete,  with  a density 
of  2.4,  is  the  cheapest  approach,  a total  thickness  of 
1.25  meters  («  49  inch)  is  required. 

Although  the  angular  yield  of  bremsstrahlung  is 
primarily  forward,  the  large  angle  component  is  still 
lto  10%  of  the  forward  yield.  Thus  nearly  4rr  shielding 
is  called  for.  Fortunately,  this  is  only  required 


immediately  around  the  beam  dump.  It  represents  a cube 
of  roughly  8 feet  on  a side.  The  floor  loading  pro- 


duced is  roughly  8.3  psi  (1,200  pounds/sq.  ft.).  Use  of 
a steel  pallet  to  distribute  this  load  over  a larger  area 


may  be  required. 
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Assuming  that  a 30  amp,  30  MeV  proton  beam  is 
collectively  accelerated  once  every  5 minutes,  the 
protons  striking  the  graphite  beam  stop  will  undergo 
nuclear  interactions  which  generate  high  energy  (>  1 MeV) 
neutrons . To  a very  good  approximation  the  neutron 
yield  is  isotropic.  The  neutron  yield  for  30  MeV  protons 
on  graphite  has  been  measured  to  be  (Ref.  11) 

Y ~ 4rr  x 1010  neutrons/|jC  of  protons 

Since 

Q ~ 6 y,C 

then  at  1 meter  the  neutron  flux  is 

YQ  _1_ 

^ 4tt  R2  At 

7 x 107  neutrons/cm2  • Hr  (W6) 

Given  that  (Ref.  11) 

2.88  x 104  neutrons/cm2  • Hr  = 1 

Hr 

11.  Burrill,  E.  A.,  "Neutron  Production  and  Protection, 
HVEC  Report,  High  Voltage  Engineering  Corp.,  Burlington 
Mass . 
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